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Preface 


Sliding mode control (SMC) has been studied since the 1950s and widely used in 
practical applications, such as robot manipulators, aircraft, underwater vehicles, 
spacecraft, flexible space structures, electrical motors, automotive engines, power 
electronics systems, and fuel cell power systems, due to its insensitivity to 
parameter variations and robustness against external disturbances. Among these 
applications, many research endeavors have been focused during the past decades, 
on control and observation problems in fuel cell power systems and power elec- 
tronics systems, using sliding mode technique. One of the most distinguished 
properties of SMC is that it utilizes a discontinuous control action, which switches 
between two distinctively different system structures such that a new type of system 
motion, called sliding mode, exists in a specified manifold. This particular char- 
acteristic of the motion in the manifold is its insensitivity to parameter variations, 
and complete rejection of external disturbances. 

Polymer electrolyte membrane fuel cells (PEMFCs) have emerged as the most 
prominent technology for energizing future’s automotive world. They are clean, 
quiet and efficient, and have been widely studied in automotive applications over 
the past two decades due to their relatively small size, lightweight and easy man- 
ufacture. While there are still major issues concerning cost, liability, and durability 
to be addressed before they become a widely used alternative to internal com- 
bustion engines (ICEs), fuel cells are expected to lead the world towards 
fossil-fuel-independent hydrogen economy in terms of energy and electro-mobility. 
One hindrance of fuel cells in general, as a kind of independent electrical power 
sources, is that their dynamic response is slow. Therefore, a fuel-cell-based power 
system requires additional storage elements with fast response time in order to 
handle rapid load variations. 

This book is organized as follows. In Chap. 2, a brief review of sliding mode 
control for both linear and nonlinear systems are presented. 

In Chap. 3, we focus on sliding mode observer (SMO) design for both linear 
uncertain systems and nonlinear systems. Two kinds of SMO design based on 
Utkin’s and Lyapunov’s methods for uncertain linear systems are recalled. Then, 
the SMO designs are extended to three forms of nonlinear systems (companion 
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form, triangular input form and algebraical observable form). The applications of 
SMO based fault detection and isolation (FDI) have been presented and fault 
reconstruction via traditional SMOs is discussed. Then, our theoretical contribution 
in Second-Order SMO-based FDI is introduced. A novel adaptive SOSM Observer 
is developed and its application in FDI is presented. In the end, two illustrative 
examples are shown to both SMO designs and their applications in FDI. 

In Chap. 4, the model of a PEMFC system which includes a stack voltage model 
and an air feed system model is discussed. The air feed system is modeled as a 
four-state model which considers the dynamics of oxygen partial pressure, nitrogen 
partial pressure, compressor speed, and supply manifold pressure. Then, a real-time 
PEMFC emulator is designed using experimental data obtained from a 33 kW 
PEMEFC unit containing 90 cells in series. Finally, the proposed air feed system 
model is validated experimentally through the hardware-in-loop (HIL) test bench 
which consists of a physical air-feed system, based on a commercial twin screw 
compressor and a real-time PEMFC emulator. 

In Chap. 5, we focus on the problem of robust control in PEMFC system, in 
order to maximize the fuel cell net power and avoid the oxygen starvation by 
regulating the oxygen excess ratio to its desired value during fast load variations. 
The oxygen excess ratio is estimated via an extended state observer from the 
measurements of the compressor flow rate, the load current and supply manifold 
pressure. An HIL test bench which consists of a commercial twin screw air com- 
pressor and a real-time fuel cell emulation system, is used to validate the perfor- 
mance of the proposed extended state observer (ESO)-based second-order sliding 
mode (SOSM) controller. The experimental results show that the proposed con- 
troller is robust and has a good transient performance in the presence of load 
variations and parametric uncertainties. 

In Chap. 6, the estimation problem of the PEMFC system is presented. First, the 
design of an SOSM observer is presented for the PEMFC air-feed system, in order 
to estimate the hydrogen partial pressure in the anode channel of the PEMFC, using 
the measurements of stack voltage, stack current, anode pressure and anode inlet 
pressure. The robustness of this observer against parametric uncertainties and load 
variations is studied, and the finite time convergence property is proved via 
Lyapunov analysis. Then, an algebraical observer is designed for the partial pres- 
sures of oxygen and nitrogen in the cathode of the PEMFC. The states of the 
PEMEC air-feed system are presented in terms of a static diffeomorphism involving 
the system outputs (compressor flow rate and supply manifold pressure) and their 
time derivatives, respectively. The implementation of the algebraical observer on 
the HIL test bench is described. The effectiveness and robustness of the observer are 
validated experimentally. 

In Chap. 7, we focus on the fault diagnosis problem of the PEMFC system, 
considering a fault scenario of sudden air leak in the air supply manifold. The 
design of an SOSM observer is presented for the PEMFC air-feed system, con- 
sidering state estimation, parameter identification and fault reconstruction problems 
simultaneously. An adaptive-gain SOSM observer is developed for observing the 
system states, where the adaptive law estimates the uncertain parameters. The 
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oxygen starvation phenomenon is monitored through an estimated performance 
variable (oxygen excess ratio). Satisfactory experimental results are obtained to 
show the effectiveness of the proposed observer. The effect of parameter variations 
and measurement noise are considered during the observer designs. 

In Chap. 8, we concentrate on the control problem on the power side of the 
PEMEC system. Fuel cell power system requires power conditioning circuits with 
precise power control algorithms behind them such that the output power is com- 
patible with the constraints of the power bus. These circuits include DC/DC con- 
verters which are used to control the voltages of the fuel cell and storage elements 
to the desired value. Furthermore, robust control design for DC/DC power con- 
verters is discussed, using sliding mode techniques. The main focus is on the 
necessary modification and improvement of conventional sliding mode methods for 
their applications to fuel cell power systems. 

In Chap. 9, we consider the control and observation problems of three-phase 
AC/DC power converters. The controller design is based on the system model and 
has a cascaded structure which consists of two control loops. The outer loop reg- 
ulates the DC-link capacitor voltage and the power factor providing the current 
references for the inner control loop. The current control loop tracks the actual 
currents to their desired values. To design the proposed controller, the load con- 
nected to the DC-link capacitor is considered as a disturbance, which directly 
affects the performance of the whole system. Theoretical analysis is given to show 
the closed-loop behavior of the proposed controller and experimental results are 
presented to validate the control algorithm under a real power converter prototype. 

This book is a research monograph whose intended audience is graduate and 
postgraduate students, academics, scientists and engineers who are working in the 
field. 
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Chapter 1 M®) 
General Introduction ectics 


This chapter introduces in detail the fuel-cell-based industrial power systems and 
several additional types of equipment required to make the fuel cell work at the 
optimal operating point. It also stresses some system failures or mechanical faults of 
the vulnerable power generation systems that can cause the shutdown or the perma- 
nent damage of the fuel cell. This monograph is focused on the polymer electrolyte 
membrane fuel cell air-feed system and power electronics systems. The motivation 
behind concentrating on these systems and the importance of their observation are 
discussed, and the outline of the monograph is described in this chapter. 


1.1 Fuel-Cell-Based Industrial Power Systems 


Polymer electrolyte membrane fuel cells (PEMFCs) have emerged as the most promi- 
nent technology for energizing future’s automotive world. They are clean, quiet and 
efficient, and have been widely studied in automotive applications over the past two 
decades due to their relatively small size, light weight and easy manufacture [1, 5, 
12, 13]. While there are still major issues concerning cost, liability and durability to 
be addressed before they become a widely used alternative to internal combustion 
engines (ICE), fuel cells are expected to lead the world towards fossil-fuel indepen- 
dent hydrogen economy in terms of energy and electro-mobility. 

One hindrance of fuel cells in general, as a kind of independent electrical power 
sources, is that their dynamic response is slow. Therefore a fuel-cell-based power 
system requires additional storage elements with fast response time in order to handle 
rapid load variations. The most common elements are rechargeable batteries and 
super capacitors [15]. Hybrid electrical vehicles (HEV), such as Audi Q5, carry 
high-power batteries as well that can share the load with fuel cells. Obviously, their 
integration in the power system introduces additional converters in order to control 
their charging and discharging on the power bus. For example, let us consider a 
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Fig. 1.1 PEMFC power system 


typical fuel-cell based hybrid automotive power system, shown in Fig. 1.1. The power 
electronics topology has numerous interconnected components, i.e. 


e Power elements (fuel cell stack and battery) 

A boost-type unidirectional DC/DC converter, for boosting fuel cell voltage to 
meet the power bus requirements 

A boost-type bidirectional DC/DC converter for connecting the battery to the bus 
A three-phase AC/DC rectifier, for occasional battery charging through external 
power source 

A three-phase DC/AC inverter for vehicle propulsion and traction motor(s) 

A multi-cell converter for DC loads (power windows, windshield wipers etc.). 


It is clear that such a system requires sophisticated bus control algorithms in the 
background. In addition, automotive fuel cell applications have more rigorous oper- 
ating requirements than stationary applications [3], as the risks of mechanical faults, 
leakage and explosion are higher in mobile systems. Technological malfunction can 
also turn the fuel-cell car into a potential hazard of life. These applications need 
precise control of performance, in order to guarantee reliability, health and safety 
of both, the fuel cell and the user. Along with control, health monitoring and safety 
systems are essential for the application of fuel cells in automobiles. These monitor- 
ing and control systems need to be precise, yet cost-effective and computationally 
simple, so that they can be easily implemented on commercial automotive embedded 
system platforms. Both control and health-monitoring systems require precise mea- 
surements of different physical quantities in the fuel cell. However, it is not always 
possible to use sensors for measurements, either due to prohibitive costs of the sens- 
ing technology or because the quantity is not directly measurable, especially in the 
conditions of humidified gas streams inside the fuel cell stack. In these cases, state 
observers serve as a replacement for physical sensors, for obtaining the unavailable 
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quantities, are of great interest. Observers form the heart of a general control problem 
and serve various purposes, such as identification, monitoring and control the system 
(Fig. 1.2). This book is dedicated to the problems associated with state observation 
in fuel cell systems. 


1.2. Motivation and Features of the Book 


We have seen that fuel-cell-based power generation systems are complex. Several 
additional types of equipment are required to make the fuel cell work at the optimal 
operating point. They are therefore vulnerable to system failures or mechanical faults 
that can cause the shutdown or the permanent damage of the fuel cell. Thus, reliable 
state observation and fault detection and isolation (FDI) schemes for such systems 
are necessary. As the entire fuel cell system cannot be covered in one study, the 
work carried out for this book is focused on the PEMFC air-feed system and power 
electronics systems. The motivation behind concentrating on these systems and the 
importance of their observation are discussed in the following subsections. 


1.2.1 Sliding Mode Control, Observation and FDI 


Sliding mode control (SMC) is a special kind of nonlinear control which has proven 
to be an effective robust control strategy for incompletely modeled or nonlinear sys- 
tems since its first appearance in the 1950s [4, 8, 17]. One of the most distinguished 
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properties of SMC is that it utilizes a discontinuous control action which switches 
between two distinctively different system structures such that a new type of system 
motion, called sliding mode, exists in a specified manifold. This particular charac- 
teristic of the motion in the manifold is its insensitivity to parameter variations and 
complete rejection of external disturbances. SMC has been applied primarily to the 
control of variable structure systems (VSS), its analysis and design are well pre- 
sented in books, survey and tutorial papers [7, 8, 18, 19, 21], both from theoretical 
and implementation perspective. 

SMC suffers from the so-called chattering phenomenon, which is undesirable 
because it often causes control inaccuracy, high heat loss in electric circuitry, and 
high wear of moving mechanical parts [21]. The chattering in SMC systems is usually 
caused by (1) the unmodeled dynamics with small time constants, which are often 
neglected in the ideal model; and (2) utilization of digital controllers with finite 
sampling rate, which causes so-called discretization chattering. Theoretically, the 
ideal sliding mode implies infinite switching frequency. Since the conventional SMC 
action is constant within a sampling interval, switching frequency cannot exceed 
that of half the sampling frequency, which leads to chattering as well. In practical 
applications, the high switching frequency is not possible or undesirable due to 
limitations of switching devices, such as switching losses, time delay, response time 
constant, presence of dead zone, hysteresis effect and saturation of device switching 
frequency [16]. In addition, the chattering action may excite the unmodeled high- 
order dynamics, which could damage actuators, systems and even leads to unforeseen 
instability [9]. Various approaches have been developed to address these problems, 
more details can be found in [6, 8, 18, 19, 21]. 

The idea of using a dynamical system to generate estimates of the system states 
was proposed by Luenberger in 1964 for linear systems [10]. FDI is usually achieved 
by generating residual signals, obtained from the difference between the actual sys- 
tem outputs and their estimated values calculated from dynamic models. The basic 
configuration of observer-based FDI is shown in Fig. 1.3. 

Sliding mode techniques known for their insensitivity to parametric uncertainty 
and external disturbance, have been intensively studied and developed for observation 
and FDI problems, existing in the fuel cell power system. In particular, higher-order 
sliding mode (HOSM) approaches are considered as a successful technique due to 
the following advantages [11]: 


e Possible to work with reduced order observation error dynamics 
Possible to estimate the system states in finite time 

Possible to generate continuous output injection signals 
Possible to offer chattering attenuation 

Robustness with respect to parametric uncertainties. 
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Fig. 1.3. Basic configuration of observer-based FDI [14] 


1.2.2. PEMFC Power System 


Fuel cells produce electricity through hydrogen and oxygen reaction. In PEMFCs, 
the anode and cathode sides are fed by hydrogen and oxygen, respectively. In fuel 
cell automobiles, hydrogen is stored in pressurized cylinders whereas air is used as 
an oxygen source, pumped into the cathode by a compressor. The air-feed system 
introduces an interesting challenge in the overall PEMFC system performance. As the 
PEMEFC system works as an autonomous power plant in automobiles, the compressor 
motor is also powered by the PEMFC. Therefore, the net power of the system is the 
difference of the power produced by the fuel cell and that consumed by the air-feed 
system (consumption by the other auxiliary systems is negligible). Experimental 
studies have shown that the air-feed system can consume up to 30% of the fuel- 
cell power under high load conditions [20]. Therefore, it needs to be operated at its 
optimal point, at which it supplies just sufficient oxygen necessary for the hydrogen 
and oxygen reaction. 

Unfortunately, such type of control is not possible without the knowledge of exact 
oxygen partial pressure. Conventional sensors can only give the total air pressure 
inside the cathode, which contains the partial pressures of other mixture gases such 
as nitrogen, carbon dioxide, etc. Imprecise knowledge of the oxygen quantity in the 
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cathode can lead to serious problems, such as oxygen starvation during load transi- 
tions and hot-spots on the membrane surface [12]. To solve this problem, observers 
can be adopted. Observers can serve two important roles in the air-feed system. 
First, they can provide a precise estimate of the oxygen partial pressure. Second, 
they can detect any immediate variations in the nominal values of pressures through- 
out the air-feed system in order to identify anomalous behavior, thereby detecting 
and identifying faults and failures. 


1.2.3. Power Electronics System 


Fuel cell power systems require power conditioning circuits with precise power 
control algorithms behind them in order for the output power to be compatible with 
the constraints of the power bus. Power converters are key components in managing 
the energy flow through such hybrid systems. Control of power electronic systems 
requires the knowledge of several states, whereas practical systems are equipped with 
a limited number of voltage sensors due to cost concerns. An observer can augment 
the number of available-for-control states by using the output voltage measurements 
to observe the unknown current values. Along with the objective for control, they 
can also be used to identify system faults by detecting abnormal currents. 


1.3 Outline of the Book 


This book is organized as follows. In Chap. 2, a brief review of SMC for both linear 
and nonlinear systems are presented. The main contents of this monograph are shown 
in Fig. 1.4. 


Chapter 3 firstly provides SMO designs for both linear uncertain systems and 
nonlinear systems. Two kinds of SMO design based on Utkin’s and Lyapunov’s 
methods for uncertain linear systems are recalled. Then, the SMO designs are 
extended to three forms of nonlinear systems (companion form, triangular input 
form and algebraical observable form). The applications of SMO based FDI are 
presented and fault reconstruction via traditional SMOs is discussed. Then, our 
theoretical contribution in Second Order SMO based FDI is introduced. A novel 
adaptive SOSM Observer is developed and its application in FDI is presented. In 
the end, two illustrative examples are shown concerning both SMO designs and 
their applications in FDI. 

Chapter 4 describes the model of a PEMFC system which includes a stack voltage 
model and an air feed system model. The air feed system is modeled as a four- 
state model which considers the dynamics of oxygen partial pressure, nitrogen 
partial pressure, compressor speed, and supply manifold pressure. Then, a real- 
time PEMFC emulator is designed using experimental data obtained from a 33kW 
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PEMEFC unit containing 90 cells in series. Finally, the proposed air feed system 
model is validated experimentally through the HIL test bench which consists of 
a physical air-feed system, based on a commercial twin screw compressor and a 
real-time PEMFC emulator. 

Chapter 5 focuses on the problem of robust control in PEMFC system, in order 
to maximize the fuel cell net power and avoid the oxygen starvation by regulating 
the oxygen excess ratio to its desired value during fast load variations. The oxygen 
excess ratio is estimated via an extended state observer from the measurements 
of the compressor flow rate, the load current and supply manifold pressure. An 
HIL test bench which consists of a commercial twin screw air compressor and a 
real-time fuel cell emulation system, is used to validate the performance of the 
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proposed ESO-based SOSM controller. The experimental results show that the 
proposed controller is robust and has a good transient performance in the presence 
of load variations and parametric uncertainties. 

Chapter 6 studied the estimation problem of the PEMFC system. First, the design 
of an SOSM observer is presented for the PEMFC air-feed system, in order to 
estimate the hydrogen partial pressure in the anode channel of the PEMFC, using 
the measurements of stack voltage, stack current, anode pressure, and anode inlet 
pressure. The robustness of this observer against parametric uncertainties and load 
variations is studied, and the finite time convergence property is proved via Lya- 
punov analysis. Then, an algebraical observer is designed for the partial pressures 
of oxygen and nitrogen in the cathode of the PEMFC. The states of the PEMFC 
air-feed system are presented in terms of a static diffeomorphism involving the 
system outputs (compressor flow rate and supply manifold pressure) and their time 
derivatives, respectively. The implementation of the algebraical observer on the 
HIL test bench is described. The effectiveness and robustness of the observer are 
validated experimentally. 

Chapter 7 studied the fault diagnosis problem of the PEMFC system, considering 
a fault scenario of sudden air leak in the air supply manifold. The design of an 
SOSM observer is presented for the PEMFC air-feed system, considering state 
estimation, parameter identification and fault reconstruction problems simultane- 
ously. An adaptive-gain SOSM observer is developed for observing the system 
states, where the adaptive law estimates the uncertain parameters. The oxygen 
starvation phenomenon is monitored through an estimated performance variable 
(oxygen excess ratio). Satisfactory experimental results are obtained to show the 
effectiveness of the proposed observer. The effect of parameter variations and 
measurement noise are considered during the observer designs. 

Chapter 8 focuses on the power side of the PEMFC system. Fuel cell power system 
requires power conditioning circuits with precise power control algorithms behind 
them such that the output power is compatible with the constraints of the power bus. 
These circuits include DC/DC converters which are used to control the voltages of 
the fuel cell and storage elements to the desired value. Robust control design for 
DC/DC power converters is discussed, using sliding mode techniques. The main 
focus is on the necessary modification and improvement of conventional sliding 
mode methods for their applications to fuel cell power systems. 

Chapter 9 considers the control and observation problems of three-phase AC/DC 
power converters. The controller design is based on the system model and has a 
cascaded structure which consists of two control loops. The outer loop regulates 
the dc-link capacitor voltage and the power factor providing the current references 
for the inner control loop. The current control loop tracks the actual currents to their 
desired values. To design the proposed controller, the load connected to the dc-link 
capacitor is considered as a disturbance, which directly affects the performance of 
the whole system. Theoretical analysis is given to show the closed-loop behavior 
of the proposed controller and experimental results are presented to validate the 
control algorithm under a real power converter prototype. 
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Chapter 2 ®) 
Basic Theory of Sliding Mode Control creek 


This chapter introduces some fundamentals of SMC including SMC design methods 
and main approaches to alleviate or limit the chattering. Moreover, some second- 
order SMC algorithms and definitions are introduced. It also presents the basics for 
the theoretical results developed in the subsequent chapters. 


2.1 General Development of Sliding Mode Control 


SMC is a special kind of nonlinear control which has proven to be an effective 
robust control strategy for incompletely modeled or nonlinear systems since its first 
appearance in the 1950s [15, 21, 35]. One of the most distinguished properties of 
SMC is that it utilizes a discontinuous control action which switches between two 
distinctively different system structures such that a new type of system motion, called 
sliding mode, exists in a specified manifold. This particular characteristic of the 
motion in the manifold provides insensitivity to parameter variations and complete 
rejection of external disturbances. The analysis and design of SMC have been well 
presented in the literature [19, 21, 36, 38], from both theoretical and implementation 
perspectives. 

SMC, as one of the robust control strategies, is famous for its strong robust- 
ness to system uncertainties in a sliding motion. However, the uncertainties should 
satisfy the so-called ‘matching’ condition, that is, the uncertainties act implicitly 
within channels of the control input. If a system has mismatched uncertainties in 
the state matrix or/and the input matrix, the conventional SMC approaches are not 
directly applicable. Therefore, in the past two decades, many researchers investi- 
gated the SMC of uncertain systems with mismatched uncertainties/disturbances, 
see for example, [7, 8] and references therein. To mention a few, in [23], the 
SMC of uncertain second-order single input systems with mismatched uncertainties, 
was considered; In [8, 34], the authors investigated the SMC design for uncertain 
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systems, in which the uncertainties are mismatched and exist only in state matrix. 
The related approaches were then developed in [9, 10] to deal with a more compli- 
cated case that the mismatched uncertainties are involved in not only the state matrix 
but the input matrix. In addition, the integral SMC techniques were extensively 
used to deal with uncertain systems with mismatched uncertainties, see for example, 
[5, 6], and some other SMC approaches to deal with uncertain systems can be found 
in [13, 32]. 

In the past two decades, SMC has successfully been applied to a wide variety of 
practical systems such as robot manipulators, aircraft, underwater vehicles, space- 
craft, flexible space structures, electrical motors, power systems, and automotive 
engines [11]. In general, SMC suffers from the so-called chattering phenomenon, 
which is undesirable because it often causes control inaccuracy, high heat loss in 
electric circuitry, and high wear of moving mechanical parts [38]. In addition, the 
chattering action may excite the unmodeled high-order dynamics, which could dam- 
age actuators, systems and even leads to unforeseen instability [22]. The chattering 
in SMC systems is usually caused mainly by the following reasons: 


1. Utilization of digital controllers with finite sampling rate, which causes the so- 
called discretization chattering. Theoretically, the ideal sliding mode implies infi- 
nite switching frequency. Since the conventional SMC action is constant within 
a sampling interval, switching frequency cannot exceed that of half the sampling 
frequency, which leads to chattering. 

2. The unmodeled dynamics with small time constants, which are often neglected 
in the ideal model. 


In this section, we will present some preliminary background and fundamental 
theory of SMC at first, which is helpful to some readers who have no or limited knowl- 
edge on SMC, and then an overview of recent development of SMC methodologies 
and chattering reduction methods will be given. 


2.2 Fundamental Theory of SMC 


Let us consider the following nonlinear system: 
X(t) = fx, t) + g(x, tute), (2.1) 
where x € R” is the state variable vector, u(t) € R” is the control input, f(-, -) and 


g(-, -) are continuous functions in x, u and ¢t vector fields [12, 36]. 
The sliding mode controller 


u(t) = [wi (t) u(t) «+» Un (t)] (2.2) 
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+ if 5; 
ince de Oe OS oa Bi ay (2.3) 
u; (t), if s;(x) <0, 


where u; (t) # u,; (t), and s(x) € R” called sliding manifold is the switching vec- 
tor function s(x) = [ si(x) $2(X) +++ Sm(X) le It undergoes discontinuities on the 
surface s(x) = 0. 

Note that the SMC law (2.3) is designed to ensure that the sliding surface (s(x) = 
0) is reached and then motion on the sliding surface is maintained. This means that 
the so-called ‘reachability condition’ should be satisfied by manipulating the control 
law u(t). The necessary and sufficient condition for the system (2.1) to satisfy the 
reachability condition is expressed as 


s(x)8(x) <0. (2.4) 


Condition (2.4) guarantees that the trajectory of the system states always points 
towards the sliding surface. A more strict reachability condition called “7-condition’ 
is given as follows 


5(x)S(x) < —nls@, DI, (2.5) 


where 77 is a positive scalar. Condition (2.5) ensures that the sliding surface is reached 
in finite time despite the presence of uncertainty. The system state trajectories under 
the reaching condition (2.4) or (2.5) are called the reaching phases [12, 19]. 


2.2.1 SMC Design Methods 


Several SMC design methods have been proposed in literature which mainly consist 
of two steps [12, 36]: 


Step 1: Design a sliding manifold s(x) which provides desired performance in the 
sliding mode, such as stability, disturbance rejection capability and tracking; 
Step 2: Design a discontinuous feedback control u(t) which will force the system 
states to reach the sliding manifold in finite time, thus the desired performance is 

attained and maintained. 


For ease of implementation, the sliding variable s;(x), i = 1,2,...,m is chosen 
as a linear combination of the state variables, expressed as 


5;(x) = > ajx;(t), (2.6) 
j=l 
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where a; denotes the sliding coefficients and x;(t) € x(t). The main objective of 
the sliding mode controller is to drive the system state trajectories onto the specified 
sliding surface in a finite time and maintained there for all subsequent time. Typical 
SMC strategies will be introduced in the following part. 


Equivalent-Control-Based Design 


For the system (2.1), assuming that the term $ 9 g(x, t) is non-singular, the control 
law u(t) is designed as follows, 


U(t) = Ueg(t) + uN (t), (2.7) 
where ueq(t) represents a continuous component and u(t) represents a discontinu- 
ous component. 


The equivalent control ueq(t) is derived from the so-called equivalent control 
method, i.e., in the case when s(x) = s(x) = 0. Thus, ueq(t) is calculated as 


a To 
Weg(t) = — (Foe. ») at t). (2.8) 


Substituting the above equivalent control (2.8) into the original system (2.1), it fol- 
lows that the motion of sliding mode is determined by 


a 19 
i() = c =9 (Shae ») =| fxd, (2.9) 


where (2.9) is considered as the equation of the sliding mode in the manifold s(x) = 0. 
The high frequency switching action uv y(t) can be designed as 


= 
uny(t) = -3 (> g(x, ») sign(s(x)), G>0 (2.10) 


such that the derivative of the Lyapunov function V = $sT(x)s(x) is negative, that 


is 
V Ss OG) =5'@) ote OO atx, t)uy(t) < —B\ls@)Il. (2.11) 


Remark 2.1 The physical meaning of the equivalent control can be interpreted as 
the low-frequency component of the discontinuous control law u(t), because the 
high-frequency uy (t) can be filtered out by a low pass filter of the system 


Tz2+z=u(t), T<1, (2.12) 


which means Z ~ Ugg. 


2.2 Fundamental Theory of SMC 15 


Reaching Law Approach 


The reaching law specifies the dynamics of a switching function, which can be 
described by the following differential equation: 


S(x) = —TYsign(s(x)) — Kg(s(x)), (2.13) 
where 


Y = diag{é1, €2,.--,€m}, 6; > 9, 
K = diag{k,, ko, ...,km}, ki > 0, 


sign(s, (x)) gi(si(x)) 

sign(s2(x)) g2(s2(x)) 
sign(s(x)) = : » g(s(x)) = ; 

sign (Sm (x)) Im (Sm (x)) 


gi(0) = 0, s;x)gi(s;(x)) > 0, i =1,...,m. 


Equation (2.13) is a general form of reaching law, and some special cases are 


1. The constant rate reaching law: 
S(x) = —Ysign(s(x)). 
2. The constant plus proportional rate reaching law: 
S(x) = —Ysign(s(x)) — Ks(t). 
3. The power rate reaching law: 
S)(x) = —e;|s;(x)|“sign(s;(x)), O<a <1. 


The reaching law approach not only guarantees the reaching condition but also spec- 
ifies the dynamics of the motion during the reaching phase [19]. 


2.2.2 Chattering Problem 


Chattering problem is one of the main obstacles for applying SMC to real appli- 
cations. The chattering in SMC systems is usually caused by (1) the unmodeled 
dynamics with small time constants, which are often neglected in the ideal model; 
and (2) utilization of digital controllers with finite sampling rate, which causes so 
called ‘discretization chattering’. From control engineers’ point of view, chattering is 
undesirable because it often causes control inaccuracy, high heat loss in electric cir- 
cuitry, and high wear of moving mechanical parts. In addition, the chattering action 
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Table 2.1 Main approaches to alleviate or limit the chattering problems 


Approach 


Boundary layer 
approach 


Operation principle 
To insert a boundary layer near the sliding surface so that a continuous control 
action replaces the discontinuous one when the system is inside the boundary 
layer [33]. For this purpose, the discontinuous component of the controller: 
uy (t) = —Kssign(s(x)) 
is often replaced by the saturation control: 

wok 5(x) 
nO ESTO 
for some, preferably small, 6 > 0. 
The boundary layer approach has been utilized extensively to the practical 
applications. However, this method has some disadvantages such as: 
e It may give a chattering-free system but a finite steady-state error must exist 
e The boundary layer thickness has the trade-off relation between control 
performance of SMC and chattering mitigation 
e Within the boundary layer, the characteristics of robustness and the accuracy 
of the system are no longer assured 


Reaching law 
approach 


Since the amplitude of chattering depends on the magnitude of control, the 
intuitive way of chattering reduction is to decrease the amplitude of the 
discontinuous control [17]. This technique affects the robustness property of the 
controller and degrades the transient response of the system. Therefore, exists a 
trade-off between the chattering reduction and the system performance. A 
compromised approach is to decrease the amplitude of the discontinuous 
control, uy (t), when the system state trajectories are near to sliding surface (to 
reduce the chattering), and to increase the amplitude when the system states are 
far from the sliding surface 


Dynamic SMC 
approach 


The main idea of the dynamic SMC approach is to insert an integrator (or any 
other strictly proper low-pass filter) between the SMC and the controlled plant 
[4]. The concept is illustrated here: 


Dynamic SMC 


Augmented System 


The time derivative of the control input, w, is treated as the new control input 
for the augmented system. Since the low-pass integrator filters out the high 
frequency chattering in w, the control input to the real plant, u, becomes 
continuous which offers a possibility to reduce chattering. Such a method can 
eliminate chattering and ensure zero steady-state error, however, it should be 
noted that the system order is increased by one and the transient responses must 
be degraded [37] 


may excite the unmodeled high-order dynamics, which probably leads to unfore- 
seen instability [24]. Therefore, various methods have been proposed in literature to 
reduce or soften the chattering action [1, 4]. Among others, the main approaches to 
avoid or limit the chattering problems are shown in the Table 2.1. 
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2.3 SOSM Control 


Apart from the above-mentioned chattering elimination/reduction approaches, 
HOSM technique is an effective method for chattering attenuation [3, 14, 25]. In 
this technique, the discontinuous control is applied on a higher time derivative of the 
sliding variable, such that both the sliding variable and its higher time derivatives 
converge to the origin. It is of great interest to see that the discontinuous control does 
not act upon the system input directly, and the real control input single is continu- 
ous, therefore, chattering effect is alleviated while the properties of robustness and 
finite time convergence of classical SMC are retained. Among various HOSM con- 
trol techniques, SOSM control technique is the most popular and effective method 
which forces the sliding variable and its first order time derivative to zero in finite 
time [3, 26-28]. 


2.3.1 Definitions 


Consider a discontinuous differential equation in the sense of Filippov [16] 
x = v(x), (2.14) 


where x € X C R’ is the state vector, v is a locally bounded measurable (Lebesgue) 
vector function. The equation can be replaced by an equivalent differential inclusion 


x € V(x). (2.15) 


If the vector-field v is continuous almost everywhere, V(x) is the convex closure 
of the set of all possible limits of V(y) as y > x, while {y} are continuity points 
of v. Solutions of the equation are defined as absolutely continuous functions x(t), 
satisfying the differential inclusion almost everywhere [31]. Let a constraint function 
given by 


s(t, x(t)) = 0, (2.16) 


where s : R” — Risa sufficiently smooth function. 


Definition 2.1 ((31]) Suppose that 


e Successive total time derivatives s and 5 are continuous functions of the system 
state variables. In other words, the discontinuity does not appear in the s; 
e The set s = s = 0 is non-empty and consists of Filippov’s trajectories [31]. 


Then, the motion on the set s = 5 = 0 is called a second sliding mode with respect 
to the constraint function s (Fig. 2.1). 
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s=s=0 


Fig. 2.1 SOSM trajectory 


2.3.2 SOSM Dynamics 


Consider the following nonlinear system [25] 


X(t) = f(t, x(t), u), (2.17) 
s(t) =s(t,x), 
where x € X C R” is the state vector, u € U C R is the bounded control input, t 
is the independent variable time and f is a sufficiently smooth uncertain vector 


function. The control objective is to force the sliding variable and its time derivative 
s(t) and s(t) to zero in finite time, i.e. 


s(t) = s(t) =0. (2.18) 


Assume that the control task is fulfilled by its zero dynamics with respect to the 
sliding variable s(t, x). By differentiating the sliding variable s twice, 


0 0 
S(t) rr Ap te ee 


(2.19) 


5(0) OS Rie Sera See. 
Ot Ox Ou 
————— a 


y(t,x) Vt.x) 


Depending on the relative degree [20] of the nonlinear SISO system (2.17), two 
cases are considered 
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- 2 _ * 0 . 7 
e Case a: relative degree r = 1, i.e., ays # 0; , 
; : — at bs ite 
e Case b: relative degree r = 2, ie., a3 = 0, am #0. 


Relative Degree 1 


In this case, the control problem can be solved by the classical first order SMC, 
nevertheless second order SMC can be used in order to avoid chattering. Shortly 
speaking, the time derivative of the control u(t) may be considered as the actual 
control variable. A discontinuous control u(t) steers the sliding variable s and its 
time derivative s to zero, so that the plant control u is continuous and the chattering 
is avoided [4, 25]. 

The second time derivative 5 (2.19) is described by the following equation 


S= p(t, x) + y(t, x)u(t), (2.20) 


where y(t, x) and y(t, x) are some bounded functions. The following conditions are 
assumed [25]: 


1. The control values belong to the set U = {u : |u| < Uy}, where Uy = Constant 
> 1. 

2. There exists wu; € (0, 1) such that for any continuous function u(t) with |u(t)| > 
uy, there is t,, such that s(t)u(t) > O for eacht > f,. 

3. There exist positive constants sg, K,, and Ky such that if |s(¢, x)| < so, then 


0< Km < y(t,x) < Ky, Vue U, x EX, (2.21) 


and the inequality |u| > uo entails su > 0. 
4. There exists constant C such that within the region |s(t, x)| < so the following 
inequality holds, 


lp(t,x)| <C, Vue U, x EX. (2.22) 


Condition 2 means that there exists a proper control u(t) forcing the sliding variable 
into a set for any initial value of state, given that the boundedness conditions on the 
sliding dynamics defined by Conditions 3 and 4 are satisfied. It follows from (2.20), 
(2.21) and (2.22) that all solutions satisfy the differential inclusion 


8 €[-C,C]+[Kn, Kylu(t). (2.23) 


Relative Degree 2 


In this case, the control problem statement can be derived by considering the variable 
u as a State variable and u as the actual control. Suppose the system dynamics (2.17) 
is affine in the control law, i.e., 


ft, x@), u) = a(t, x) + b(t, x)u, (2.24) 
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where a : R"*+! — R" and db: R"*+! — R" are sufficiently smooth uncertain vector 
functions. Equation (2.19) can be rewritten as 


s(t) = og x)+ aes x)a(t,x)+ OG. x)b(t, x)u (2.25) 
Ot Ox Ox 


0 O 
= ape) + a ee): 


. oe a @) 
S(t) = Bs *) + Be a ee) 


2 2 
+ lax st, x) +a'@, aa a(t, +5 a ah 23 Fm e at, »| 


x Jaw. x)+ b(t, uc] 


= p(t, x) + y(t, x)u(t). 


It follows from (2.21), (2.22) and (2.25) that all solutions satisfy the following dif- 
ferential inclusion 


€ [-C, C]+ [Kn, Kylu(t). (2.26) 


2.3.3 SOSM Controllers 


In this part, the most well known SOSM controllers are introduced, i.e., the super- 
twisting controller, the twisting controller and the sub-optimal controller. These con- 
trollers are insensitive to some model perturbations and external disturbances. Given 
that the expression for the sliding manifold is known, it is possible to design the 
constant parameters of the controllers [31]. 


Super-Twisting Control Algorithm 


The super-twisting algorithm (STA) is a unique absolutely continuous sliding mode 
algorithm ensuring all the main properties of first order sliding mode control for 
system with Lipschitz continuous matched uncertainties/disturbances with bounded 
gradients [18]. The STA was developed to control systems with relative degree one 
in order to avoid chattering in variable structure control (VSC). The trajectories on 
the second sliding manifold are shown in Fig. 2.2. 

Consider the system (2.23), the control algorithm is defined as follows [25] 
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Fig. 2.2, STA phase 


trajectory s 
S 
u(t) = uj(t)+uo(t), 
: —U, if |u| > 1, 
uy(t) = : ve F 
—asign(s), otherwise, (2.27) 


—Alsol’sign(s), if |s| > so, 
uz(t) = 
—Al\s|Psign(s), otherwise, 


where a, A are positive constants and p € (0, 1). The sufficient conditions for the 
finite time convergence to the sliding manifold are 


C 4C Ky(a+C) 
Oe See, 
— Kun _ K2 Kn(a—C) 


(2.28) 

The STA does not need the evaluation of the sign of the time derivative of the 
sliding variable. For the choice p = 1, the origin is an exponentially stable equilib- 
rium point. The choice p = 0.5 assures that the maximum real second order sliding 
is achieved. For 0 < p < 0.5 the convergence to the origin is even faster. The choice 
0 < p < 1 assures the finite time convergence to the origin [25, 29]. 


Twisting Control Algorithm 


This algorithm is characterized by a twisting around the origin, shown in Fig. 2.3. 
The finite time convergence to the origin of the plane is due to the switching of the 
control amplitude between two different values. The control amplitude switch at each 
axis crossing which requires the sign of the time derivative of the sliding variable s. 

In case the relative degree r = 1. Consider the system (2.23), the twisting algo- 
rithm is defined by the following control law [25] 
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Fig. 2.3. Twisting algorithm 
phase trajectory S 


—-u, if |u| > 1 
u(t) = }—a,sign(s), if ss <0, |u| <1 (2.29) 


—aysign(s), if ss > 0, |u| < 1 


where @y > Q, > 0 and the sufficient conditions for the finite time convergence to 
the sliding manifold are 


4Ku C 
An > —, An > K.’ KinQm os KuQm + 2C. (2.30) 


0 m 


In case the relative degree r = 2. Consider the system (2.26), the twisting algo- 
rithm is defined by the following control law [27] 


u = —r\sign(s) — resign(s), 7) > r2 > 0, (2.31) 
the sufficient conditions for the finite time convergence to the sliding manifold are 
(11 +12)Km > (1 —12)Ku +2C, (11 —12)Km > C. (2.32) 


A particular case of the controller with prescribed convergence law [25, 31] is 
given by 
; : 1, x 
u = —asign (5 + A\s|#sign(s)) ,a>0, \>Oand ak, —C > 3" (2.33) 
Controller (2.33) is close to a terminal sliding mode controller [39]. 
Sub-Optimal Control Algorithm 


The SOSM controller was developed as a sub-optimal feedback implementation 
of the classical time optimal control for a double integrator [30]. This algorithm 
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Fig. 2.4 Sub-optimal 
algorithm phase trajectory S 


ensures the finite time convergence of s and s to zero, confining the trajectories 
within limit parabolic arcs (including the origin). Both twisting and jumping (in 
which s and s do not change sign) behaviors are possible (see Fig. 2.4). Unlike most 
SOSMC algorithms, sub-optimal control does not require continuous estimate of s, 
only depends on upon the instances when the value of s is zero. 

Let the relative degree r = 2. Consider the system (2.26), the control algorithm 
is defined by the following control law [2, 25] 


u(t) = —a(t) Vysign(s — =) (2.34) 


ates a*, if (s 7 55*) (s* —s) > 0, 
1, otherwise, 


where s* is the latter singular value of the function s(t) which corresponds to the zero 
value of s and a* is a positive constant. The sufficient conditions for the finite time 
convergence to the sliding manifold are 


3Kin 
ae (0,1) N10, : (2.35) 
Kyu 
C 4C 
Vu > max : , 
a* Kin 3Kin = a* Ky 


2.4 Conclusion 


In this chapter, we recalled the basic concepts of classical SMC and SOSM control 
for uncertain nonlinear systems. Three most common SOSM controllers, i.e., the 
super-twisting controller, the twisting controller and the sub-optimal controller, are 
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introduced. This control method has proven to be an effective robust control strategy 
for uncertain nonlinear systems due to its attractive characteristics such as finite 
time convergence, insensitivity to parameter variations, and complete rejection of 
external disturbances. SOSM controllers can not only reduce the chattering effect 
but also improve the accuracy of tracking performance, make this technique suitable 
for implementation of practical systems. 

In the subsequent chapters, we shall concentrate on the applications of sliding 
mode technique to the observer design and FDI. Then, SMC and FDI of a complex 
nonlinear system, i.e., PEMFC will be presented with the objective of increasing its 
efficiency and service life. 
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Chapter 3 M®) 
Sliding Mode Observer cro 
and Its Applications 


SMOs have found wide application in the areas of fault detection, fault reconstruction 
and health monitoring in recent years. Their well-known advantages are robustness 
and insensitivity to external disturbance. Higher-order SMOs have better perfor- 
mance as compared to classical sliding mode based observers because their output 
is continuous and does not require filtering. However, insofar as we are aware, their 
application in FDI has remained unstudied. In this chapter, we shall develop the 
theoretical background of SMOs and SMO-based FDI. A bibliographical study of 
existing approaches in these fields will be followed by a brief presentation of some 
established first order and second order SMO algorithms. Then, first order SMO based 
FDI methods will be demonstrated. Finally, our contribution in bridging the gap of 
second order SMO and adaptive second order SMO based FDI will be presented, 
followed by two illustrating examples. 


3.1 State of the Art 


System observation is essential for obtaining unmeasurable states for precise control 
applications. In general, they are also used for cutting costs by replacing some phys- 
ical sensors for observable states. Yet, modeling inaccuracy and parametric uncer- 
tainty in complex physical systems hinder correct state observation and induce errors. 
Sliding mode technique is known for its insensitivity to external disturbances, high 
accuracy and finite time convergence. These properties make it an excellent choice 
for observation of higher order nonlinear systems such as fuel cell systems. 

The early works were based on the assumption that the system under consideration 
is linear and that a sufficiently accurate mathematical model of the system is available. 
When the system under consideration is subject to unknown disturbances, the fault 
signal and unknown disturbance are very likely to produce a similar residual signal. 
This problem is known in the literature as robust FDI [6], which usually involves 
two steps: the first step is to decouple the faults of interest from uncertainties and the 
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second step is to generate residual signals and detect faults by decision logics. Several 
practical techniques for these steps have been proposed in contemporary literature, 
for example geometric approaches [22], H,.-optimization technique [23], observer 
based approaches (e.g. adaptive observers [35, 38], high gain observers (HGO) 
[2, 33], unknown input observers (UIO) [3, 13, 26]). 

Edwards et al. [6] proposed a fault reconstruction approach based on equivalent 
output error injection. In this method, the resulting residual signal can approximate 
the actuator fault to any required accuracy. Based on the work of [6], Tan and Edwards 
[30] proposed a sensor fault reconstruction method for well-modeled linear systems 
through the linear matrix inequality (LMI) technique. This approach is of less prac- 
tical interest, as there is no explicit consideration of disturbance or uncertainty. To 
overcome this, the same authors [31] proposed an FDI scheme for a class of linear 
systems with uncertainty, using LMI for minimizing the Lz gain between the uncer- 
tainty and the fault reconstruction signal. Linear uncertain system models can cover a 
small class of nonlinear systems by representing nonlinear parts as unknown inputs. 

Jiang et al. [14] proposed an SMO based fault estimation approach for a class 
of nonlinear systems with uncertainties. Yan and Edwards [36] proposed a precise 
fault reconstruction scheme, based on equivalent output error injection, for a class 
of nonlinear systems with uncertainty. A sufficient condition based on LMI was 
presented for the existence and stability of a robust SMO, based on strong structural 
condition of the distribution associated with uncertainties. In their later work [37], 
this structural constraint was relaxed and the fault distribution vector and the structure 
matrix of the uncertainty are allowed to be functions of the system’s output and input. 
All these works require that the bounds of the uncertainties and/or faults are known. 


3.2 First-Order SMO 


In this section, several design methods of traditional SMO will be recalled. 


3.2.1 SMO Design Based on Utkin’s Method 


Consider initially the following linear uncertain system [32]: 


x(t) = Ax(t) + Bu(t) + Gd(x, u, t), 


3.1 
y(t) = Cx(0), — 


where x € R” is the state, u € R” is the control input, y € R? is the measurable 
output. The matrices A, B and C are of appropriate dimensions. It is assumed that 
d(x, u, t) is unknown, but bounded 


Id@uHll<sp, Vt>0, (3.2) 
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where ||-|| represents the Euclidean norm. Gd(x, u, t) represents the system uncer- 
tainties, with G is a full rank matrix in R”“?. The matrices B and C are assumed 
to be of full rank and the pair (A, C) is observable. Furthermore, without loss of 
generality, the output distribution matrix C can be written as C = [¢ 1 Co ]; where 
C, € R?*®-?), C2 € R?*? and det(C2) 0. The objective is to estimate the states 
x(t) only from the measurements of input u(t) and output y(t). 

Two cases are considered: d(x, u,t) = 0 and d(x, u, t) € 0. For the first case, a 
coordinate transformation is introduced in order to facilitate the observer design 


x1 (t) _ Ino 0 o- 
| = C, a = Tx, (3.3) 
where T is non-singular. With the transformation (3.3), the system (3.1) can be 


written as 
X(t) = Ayixi(t) + Ajzy(t) + Byu(t), 


yt) = Aaixi(t) + Axy(t) + Bou(t), “ 
where 
TAT = be | , TB = B CT 7 (3.5) 
The proposed observer has the following form 
£1(t) = Au£i(t) + Ay (t) + Biu(t) + Lv, ae 


$(t) = An ky (t) + Any(t) + Bou(t) + v, 


where <j (t), (¢) represent the estimates for x(t), y(t), L ¢ R“~?*? is the constant 
gain matrix and the discontinuous terms vy is defined 


y; = osign (y(t) — Ji(t)), i=1,---,p (3.7) 


where 1; is the ith component of v. 
Denote the errors e;(t) = x)(t) — X1(t) and e,(t) = y(t) — $(t). Then, the error 
dynamical system is obtained 


é\(t) = Ayer (t) — Ly, (3.8) 
éy(t) = Arie (t) — v. (3.9) 
Thus, for a large enough scalar o, an ideal sliding motion is induced in finite time on 


the surface 
S=feeR’ : e = Ce = 0}. (3.10) 
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During the sliding motion e, = é, = 0, Eq. (3.9) is written as 
Veq = Ariei(t), (3.11) 


where 1,, represents the equivalent output error injection term which is generated by 
a low pass filter. Substituting (3.11) into Eq. (3.8), it follows that the reduced order 
sliding motion is governed by 


é1(t) = (Ar) — LAai) e1(t). (3.12) 


Since the pair (A, C) is observable, then the pair (Aj;, A21) is also observable. 
Therefore, there exists L such that the matrix A,; — LA2; is stable. Consequently, 
X(t) converges to x, (t) asymptotically. 

The main disadvantage of the above formulation is the requirement of large value 
o, to ensure sliding mode for a broad range of initial state estimation errors, especially 
when the underlying system in unstable. A trade-off between the requirement of a 
large o and its subsequent reduction to prevent excessive chattering (whilst still 
ensuring sliding mode) is usually taken into account. Slotine et al. [27] proposed a 
method which include a linear output error injection term 


S(t) = An ki(t) + Aooy(t) + Bou(t) + Grey (t) + v, (3.13) 


where the linear gain G; should be chosen to enhance the size of the so-called sliding 
patch, i.e., the domain of the state estimation error in which sliding occurs. Under 
certain conditions, the properties of global convergence of state estimation error and 
robustness can be achieved. 

For the second case d(x, u,t) 4 0. Using the transformation (3.3), the linear 
uncertain system (3.1) can be transformed into following canonical form, 


X(t) = Ay x(t) + Ayay(t) + Byu(t) + Gi d(x, u, t), 


y(t) = Arixi(t) + Any(t) + Bou(t) + Grd (x, u, t). oi 
Under the SMO (3.6), the equivalent control signal will be 
Veq = Az1e\(t) + God(x, u, t). (3.15) 
The error dynamics of e; will become 
€i(t) = (Ay — LAai) e1(t) + (G1 — LG) d(x, u, t). (3.16) 


It is clear that e;(t) will not approach zero if d(x, u,t) is nonzero. It should be 
noted that even if G, is zero, the equivalent control signal will still introduce the 
uncertainties into its observer error dynamics. A direct approach is to select the 
gain L such that G; — LG2 = 0. However, it may be very difficult to satisfy this 
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condition. Another more reasonable approach is to force the estimation error to be 
below an acceptable threshold. However, it requires that ||d(x, u, t)|| is small enough, 
as discussed in [27, 35]. 


3.2.2 SMO Design Based on Lyapunov Method 


The SMO proposed by Wallcott and Zak [34] attempts to provide exponentially con- 
vergent estimate of the state described in (3.1) in the presence of matched uncertainty. 

Recall that the pair (A, C) is assumed to be observable, thus, there exists a matrix 
K such that Ag = A — KC is Hurwitz. Therefore, for every real Symmetrical Posi- 
tive Definite (SPD) matrix Q € R”™", there exists areal SPD matrix P as the unique 
solution to the following Lyapunov equation 


Aj P+ PAo =—Q. (3.17) 
It is also assumed that the structural constraint 
= (FC); (3.18) 


is satisfied for some F' € R4*?. 
The observer proposed by [34] has the form 


R(t) = AR(t) + Bu(t) + K(y — C8) +p, (3.19) 
where 


FCe(t) : 
P(t, Y, WEG] |FCe|’ if nr #0, (3.20) 
0, otherwise. 


Denote e(t) = x(t) — x(t), then, the following error dynamical system is obtained: 

e(t) = Ape(t) —v + Gd(x, u, t). (3.21) 
It can directly proof the stability by using V(e) = e! Pe asaLyapunov function can- 
didate, it is shown that V(e) < —cV, for some positive value c, thus e(t) converges 
to zero exponentially. Furthermore, an ideal sliding motion takes place on 


Sr = {e(t) ER" : FCe(t) = O}, (3.22) 


in finite time. 
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Remark 3.1. It should be noted that in case of p > m, the sliding on Sr is not the 
same as sliding on S in (3.10). Therefore, the structure of observer (3.19) is different 
from that in (3.13). 


The main difficulty in designing the above observer is the computation the matrices 
P and F such that (3.17) and (3.18) are satisfied. In [34], a symbolic manipulation 
tool was used to solve a sequence of constraints that ensure that the principal minors 
of both P and the right hand side of (3.17) are positive and negative, respectively. It 
is convenient for low order systems, but not for high order systems. The conditions of 
the structural requirements (3.17) and (3.18) to be solvable were given by Edwards 
et al. [6] as follows: 


e rank (CG) = m; 
e any invariant zeros of (A, G, C) lie in the left half plane. 


For a square system (p = m), the above two conditions fundamentally require the 
triple (A, G, C) to be relative degree one and minimum phase. A key development in 
[6] is that there in no requirement for the pair (A, C) to be observable. The SMOs can 
be designed as long as the triple (A, G, C) satisfy the above two conditions. Details 
of the constructive design algorithms can be found in [6, 29]. Floquet and Barbot 
[9] and Floquet et al. [10] show that the relative degree condition can be relaxed 
if a classical SMO is combined with the sliding mode robust exact differentiators 
[17]. Additional independent output signals can be generated from the available 
measurements. 


3.2.3. SMO Design Based on Slotine’s Method 


Let us consider a nonlinear system in companion form [28] 


x1 = x2, 
x2 = x3, 
Le f(x, t), (3.23) 


where x? = [x1 Xqee+ | is the state vector, x; is a single measurable output and 
Ff (x, t) is a nonlinear uncertain function. 
An SMO is designed as follows 


Xx; = —aye; + X2 — kysign(e)), 


Xp = —aye; + X3 — kysign(e)), 


ky = —ane; + f — kysign(es), (3.24) 
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where e; = x; — x1, f isanestimate of f (x, t) and the constants a;,i € {1,2,--- ,n} 
are chosen to ensure asymptotic convergence for a classical Luenberger observer 
when k; = 0. The corresponding error dynamics are given by 


é; = —ayje; + e2 — ky sign(ey), 
é7 = —aze; + €3 — kysign(es), 
€, = —Aney + f= k, sign(e}), (3.25) 


where f =" f — f(x, t) is assumed to be bounded and the following condition holds 
ky = Vfl (3.26) 
The sliding condition ¢ (e1)” < Ois satisfied in the region 


e2<kit+aje, ife, >0, 
eS = ecarep, EU. (3.27) 


Therefore, the sliding mode is attained on e; = é; = 0, it follows from Eq. (3.25) 
that 


e2 — kisign(e;) = 0. (3.28) 


Substituting (3.28) into (3.25), it follows that the reduced order sliding motion is 
governed by 


a ky 

2 = &3 ki e2, 

: x ky 

ar a, (3.29) 
ky 


Mrna-| || =0. (3.30) 
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Assuming that kj, i € {2,--- ,n} are proportional with k; and the poles determining 
the dynamics of sliding patch are critically damped, i.e., are real and equal to some 
constant values —y < 0, then Slotine et al. [28] show the precision of the state 
estimation error 


< (2y)'ki, i=0,---,n-2. (3.31) 


The effect of measurement noise on SMOs was also discussed in [28], the system 
does not attain a sliding mode in the presence of noise, but remains within a region 
of the sliding patch which is determined by the bound of the noise. Moreover, it was 
demonstrated that the average dynamics can be modified by the choice of k; which 
in turn can tailor the effect of the noise on the state estimates. 


3.3. SOSM Observer 


It should be noted that the traditional first-order SMOs require low pass filters to 
obtain equivalent output injections. However, the approximation of the equivalent 
injections by low pass filters will typically introduce some delays that lead to inac- 
curate estimates or even to instability for high order systems [9]. To overcome this 
problem, continuous SOSM algorithms are used to replace the discontinuous first- 
order sliding mode, such that continuous equivalent output injection signals are 
obtained. In the following, three kinds of SOSM algorithms will be introduced. 


3.3.1 Super-Twisting Algorithm 


The STA is described by the differential equation 


= —)[¥,|? sign(%1) + %, 
X_ = —asign(X)) + d(x), (3.32) 


2 a eT : ; : 
where x = [x1 | € R’ are state variables, \, a are gains to be designed and the 
function #(x) is considered as a perturbation term, which is bounded 


|p(x)| < ®, (3.33) 


where ® is a positive constant which is assumed to be known. The solutions of (3.32) 
are all trajectories in the sense of Filippov and Arscott [8]. 
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Frequency Characteristics of STA 


It has been shown in [16] that the frequency characteristics of STA can be used to 
evaluate the performance of STA observer and determine the choice of its parameters. 
Let x; = Asin(wt) and use the describing function method, we have 


1 Qn 
= i [1 |? sign(%)) sin(wr)dwr 
T JO 


2 Tv 
=— / |A sin(wT) [2 sign (A sin(wT)) sin(wT)dwt (3.34) 
T Jo 
2. 4. , 3 
= —A? |sin(wT)|? dwr, 
T 0 
i ere nee oe a , 4 
— sign(x,) sin(wr)dwrt = — sign (A sin(wT)) sin(wT)dwr = —. (3.35) 
T JO T Jo T 
It is easy to get 2 fY |sin(wr)|dwr = +, 2 ff |sin(wr)|? dwr = 1 and 
1<A=2 f* |sinwr)|? dur < 4. (3.36) 


Thus, the describing functions of the nonlinear terms (|- [2 sign(-), sign(-)) in (3.32) 
are N,(A) = 4. and N>(A) = 4, respectively. 
Az 


The linearization of system (3.32) is 


Ge AAAS, 


2= ye + d(x). (3.37) 
TA 


The nature frequency and its damping coefficient of system (3.37) are 


— 27a — Aryv= 
Wy = mad’ S= AT: (3.38) 


From (3.37) and (3.38), we can find that, when the observation error x, is large, that 
is A is large far from the origin, the gains in (3.37) are small. On the other hand, 
when the observation error x, is small, that is A is small near the origin, the gains are 
large. Therefore, we can conclude that the behavior of STA is strong near the origin, 
contrarily, is weak far from the origin. The frequency characteristic of STA is shown 
in Fig. 3.1. 


Convergence of STA 


A strong Lyapunov function was proposed to analyze the stability and finite time 
convergence of the STA [19] 


Viet’ Pe, (3.39) 
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Bode Diagram 


Magnitude (dB) 


Phase (deg) 


Frequency (rad/sec) 


Fig. 3.1 Frequency characteristics of STA 


welll whe Yak ig 4 =r]. 
where ¢7 = [(, Go] = [ [41/? sign(%1), %] and P = 3 ae é is a SPD 


matrix. The derivative of the vector ¢ is given by 


ba | eee), Me [— Pz 
C= Zac Bog, a= [723], 
B=[01]', G = [10]¢ = ce. (3.40) 


It should be noted that VC) in (3.39) is continuous but not differentiable on the 
set S = {(X1, x2) € R?|x1 = 0}. Thus, the classical versions of Lyapunov’s theo- 
rem [1, 21] can not be applied since they require a continuously differentiable of 
at least a locally Lipschitz continuous Lyapunov function. However, by means of 
Zubov’ theorem [41] which requires only continuous Lyapunov functions, it is still 
possible to show the convergence property [19]. 


Proposition 3.1 [4, 19]. Suppose that the condition (3.33) is satisfied and there 
exists a SPD matrix P. such that the Algebraic Riccati Equation (ARE) 
A’P+PA+@°C'C+PBB'P = -Q <0 (3.41) 


is satisfied. Then, all trajectories of system (3.32) converge to the origin in finite time. 
Moreover, the quadratic form V (¢) = ¢' P¢ is a strong, robust Lyapunov function 
for the system (3.32) and the finite convergence time T is estimated 
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Ze he 
T= eV io) (3.42) 


where Xo is the initial error and y(P) is 


nin Q)ain(P) 


ay 


(3.43) 


Proof of this Proposition 3.1 can be found in [4, 19]. 


3.3.2 Modified Super-Twisting Algorithm 


It is easy to see that the standard (constant gains) STA always requires the condition 
(3.33). It does not allow to compensate uncertainties/disturbances growing in time or 
together with the state variables due to its homogeneous nature. This means that the 
standard STA can not ensure the sliding motions even for systems for which the lin- 
ear part is not exactly known. Thus, it is very important to design non-homogeneous 
extension of the standard STA allowing exact compensation of the smooth uncer- 
tainties/disturbances bounded together with their derivatives by the known functions, 
which could grow together with the state [11]. 

As shown in [20] that linear growing perturbations are included by means of the 
addition of linear terms to the nonlinear SOSM terms (SOSML). The behavior of 
the STA near the origin is significantly improved compared with the linear case. 
Conversely, the additional linear term improves the behavior of the STA when the 
states are far from the origin. In other words, the linear terms can deal with a bounded 
perturbation with linear growth in time while the nonlinear terms of STA can deal 
with a strong perturbation near the origin. Therefore, the SOSML inherits the best 
properties of both the linear and the nonlinear terms. 

The SOSML algorithm is described by the following differential equation 


= —\ 1% |? sign() — ky +X, 
%) = —asign(%1) — ke) + $(%), (3.44) 


where A, a, ky, kq are positive gains to be determined and the perturbation term $(x) 
is bounded by 


|D(x)| < 51 + 62 |i, (3.45) 


where 6; and 62 are some positive constant and are assumed to be known. The 
solutions of (3.44) are all trajectories in the sense of Filippov [8]. 
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Bode Diagram 


30 


Magnitude (dB) 


Phase (deg) 


-90 + a : ts 4s 1 3 e — : 
10 10" 10 
Frequency (rad/sec) 


Fig. 3.2 Frequency characteristics of SOSML 


Frequency Characteristics of SOSML 


Let x; = A; sin(w ft), the linearization of system (3.44) is given as 
; a 
¥ = —(AAA,?* +k) x1 + %, 


ipa -(= a kn) + $(&). (3.46) 


The nature frequency and its damping coefficient of system (3.46) are 


1 
4 NAA? +k 
eed ee. eee ee (3.47) 
TA, Gj Bee ae, 


TA, 


Comparing with STA, the behavior of SOSML is both strong when the observation 
error X, is near the origin and far from the origin. Better dynamic characteristics are 
obtained which can be seen from the frequency characteristic of SOSML shown in 
Fig. 3.2. 


Convergence of SOSML 


A strong Lyapunov function was proposed to analyze the stability and finite time 
convergence of the SOSML 


- : oe Coes rer eee —— 
Vi(&) = 20 %1| + kak? + 58 + 5 Ale sign() 4 Kady = Xa)". 
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The Eq. (3.44) can be written as a quadratic form 

Vi) = CPE, (3.48) 
where 


f= [62x B= Ea sign(X1), X1, %| : 


4datr rk, —A 
P= 5 Aky Qkaks —k)y 
—X —k, 2 
The matrix P is SPD under the condition 


4aky > (8a + 9”)k}. (3.49) 


The derivative of the vector ¢ is given by 


. 1 [73 %3 -3 0 0 0 
f= h Spee -O Het |S 0 BG), (3.50) 
Il} a 0 0 O80 1 
— eS 
Aj Ad 


Proposition 3.2. Suppose conditions (3.45) and (3.49) hold. Then, all trajectories 
of system (3.44) converge to the origin in finite time if 


Amin (£21) > O1y/ »? =F ky - 4, 


Amin (22) > b2Amax(A®P), (3.51) 


where the matrices {2,, S23 and A® are defined in (3.55) and (3.56), respectively. 
Moreover, the quadratic form V,(x) = &' P\€ is a strong, robust Lyapunov function 
for the system (3.44) and the finite convergence time T, is estimated 


2-1). 
T, = —V/ Go), (3.52) 
V1 


where Xq is the initial error and 7, is some positive constant. 


Proof. The time derivative of Lyapunov function candidate (3.48) along the trajec- 
tories of the system (3.50) is calculated as, 


0 
1 
= — ET (ATP, + PAE + ET(AZP, + PiAdE+2E7P,| 0 |. (53) 


‘= 
‘TE $(%) 
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It follows from £)é) = ce and £)63 = Sst that 


1&1 1&1 
1 
Vi = “aAie me — ETQE+ HigiE + K€" gf [0 1 OF €, (3.54) 
1 
where 
A= [-A —ky 2] > 
d + 2a 0 —X 
= . O 2ky + SkK —3ky |, 
—r —3k) 1 (3.55) 


a+2 0 0 
{27 = ky 0 ka +k -ky 
0 —k)y 1 


It is easy to show that §2;, (22 are positive definite under the condition (3.49). There- 
fore, Eq. (3.54) can be rewritten as, 


: 1 
Vi S —~Amin(QIEM? — Amin(22) Ell? + 41 U.gr ll Ell + F2Amax(A®) EI? , 


fi 
where 
4 0 0 
Ad =|03—k +10 (3.56) 
0 0 1 
Further, using the fact that 
1 
VV? 
lls léll < =4—., 
Pt) (3.57) 


llgull = 2 +45 +4, 


it follows that 


vi = (Amin > diy »? ag ke ae s) él ~ (Amin (S22) — b2Amax(A®)) él’. 


Under the conditions that cy = Amin($21) — 51,/A? + kX + 4 and cp = Amin(@2) — 
62Amax(A®) are positive, it follows 


Nic 


Vi <-1nVeZ - 2M, (3.58) 
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where 
a c2 


= 8 ———— (3.59) 
Auax (P1) Amax (Pi) 


n= 


By the comparison lemma [15], it follows that V(x(t)); and x(t) converge to zero 
1 


units of time. 


Bos 
in finite time and reaches that value at most after 7; = ee 


3.3.3. Adaptive-Gain Super-Twisting Algorithm 


The main disadvantage of STA and SOSML is that it requires the knowledge of the 

bound of uncertainty, see conditions (3.33) and (3.45). However, in many practical 

cases, this boundary can not be easily obtained. Therefore, an adaptive-gain SOSML 

is proposed, which handles the uncertainty with the unknown boundary. An adaptive 

law of the gains of the proposed algorithm is derived via the so-called “time scaling” 

approach [25], which are adapted dynamically according to the observation error. 
The adaptive-gain SOSML algorithm is described as follows 


¥ = —A() Li? sign(&1) — k\(NE + ho, 
X2 = —a(t)sign(x1) — ka(t)x1 + o(), (3.60) 


the perturbation term ¢(x) is bounded by 
IO@)| < o1 +02 1x1, (3.61) 


where a and a2 are some positive constant and are assumed to be unknown. 
The adaptive gains A(t), a(t), k,(t) and k,(t) are formulated as 


At) = Aovi(t), a(t) = aol(t), 
kG =i Oy hay = kal, (3.62) 


where Xo, Qo, ky, and k,, are positive constants to be defined and /(t) is a positive, 
time-varying, scalar function. 
The adaptive law of the time-varying function /(f) is given by: 


. k, if |x 
gai = ae (3.63) 
0, otherwise, 


where k is a positive constant. 
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Theorem 3.1 Consider system (3.60), (3.62), and (3.63). Suppose that the condition 
(3.61) hold with some unknown constants 0, and 02. The trajectories of the system 
(3.60) converge to zero in finite time if the following condition is satisfied 


Ack, > 8kx, a0 + IADK},- (3.64) 


Proof. A new state vector is introduced to represent the system (3.60) in a more 
convenient form for Lyapunov analysis. 


G1 13 (t)|¥1 |? sign(%1) 
¢=|6]= l(t), ; (3.65) 
G3 Xo 


Thus, the system in (3.60) can be rewritten as 


36 =% @ 2 eo 0 T06 
C= Te] 0 —W0/CHIM| 0 —kyI/C+] aie]. 6.66) 
1 
—ayo 0 O 0 —ka, 0 o(x) 
EE et — FO 
Aj Ad 


Then, the following Lyapunov function candidate is introduced for the system (3.66): 


1 1 2 
VQ) = 2a06f + kao + 5G +5 Aoi tkoG—G),  — G.67) 
which can be rewritten as a quadratic form 


1 4ao + MA ok ro — Xo 
VO=C PC, P = 5 Nokyg KR + kag —kry | - (3.68) 
—o —k, 2 


As (3.67) is a continuous Lyapunov function and the matrix P is positive definite. 
Taking the derivative of (3.67) along the trajectories of (3.66), 


: 1 
(a rat CINTA + OB + 7 1D an (3.69) 


where g) = [—Ao —ky, 2, Q= [1 © 0], and 
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> Dy + 2a0 0 —Xo 
oie - 0 2kay + 5kX, —3krg |. 
= —3k), 1 


3.70 
ao +242 0 0 — 


{22 = ky, 0) Ka + ki. —k), 
0 —k, 1 


It is easy to verify that §2; and {22 are positive definite matrices under the condition 


(3.64). 
Using the fact Amin(P)IClI? < V <Amax(P)IGI? and |Ilgil] = 2 +k +4. 


Equation (3.69) can be rewritten as 


. Amin 2 L Amin Q 5 
V<-1y) ( Dy! I(t) ( BD jig GT 
Ninax(P) Saal E Sete) gai 
Or a i(t) 
—C A®P —— AQ, 
= i> - 21(t) 
where 
AQ = (ao + AH)CT + 2Aok 61 G2 + Wak}, G> — A01G3 — kr 263) 
= UOC, 
4 0 0 
A® = |04-k,+10 
2 ’ 
i ; (3.72) 
4ay + 5 + Aoky + 0 0 
OS 0 Qkagk?, + Aokry + 2 0 
0 0 a) 
In view of (3.72), Eq. (3.71) becomes 
¥<— (pain — 2B) y 
~ Dakax (P) Yanin(P) (3.73) 
(1) Amin(22) 02 Amx(A®) _ I) Huu) 
Yinax(P) T(t) Amin(P) 2) Amin P) )” * 
For simplicity, we define 
Amin(Q) 71 V a aa ks 7 4 Amin (£22) 
1 = a rr | 2 = — oi oo 3 = a > Shige a 
Max(P) Min(P) Amax(P) (3,74) 
— _ Xmax(A®) _ Amax(Q) 
Va = 02 oan tS =" a 


Amin(P) , DAnii(P) 
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where 71, 72,73, Y4 and 7s are all positive constants. Thus, Eq. (3.73) can be simplified 
as 


: i 
Vs—-€On - 2) V? = (tn - i ay) (3.75) 


Because /(t) > Osuch that the terms /(t)y; — y2 and/(t)y3 — iw _ Ons are positive 
in finite time, it follows from (3.75) that 


V <—cV2 —@V, (3.76) 


where c, and c2 are positive constants. By the comparison principle [15], it follows 
that V (¢) and therefore ¢ converge to zero in finite time. Thus, Theorem 3.1 is proven. 


3.4 SMO-Based FDI 


Let us now explore the utilization of SMOs for system fault detection and recon- 
struction. Most observer-based FDI schemes generate residuals by comparing the 
measurement and its corresponding estimate provided by observers. Wrong esti- 
mates will be produced when faults occur, thus, a nonzero residual would raise an 
alarm. In this section, it is shown that SMO can not only detect the faults but also 
reconstruct the fault signal (its shape and magnitude). Fault reconstruction is of great 
interest in active fault tolerant control which can be employed in controller design [7]. 
Consider a nonlinear system [36] 


x = Ax+G(x,u)+ EW(t,x,u)+ Df, u,b), 
y=Cx, (3.77) 


where x € R”,u € R” and y € R? are the state variables, inputs and outputs, respec- 
tively. A ¢ R’*”, E € R’*’, De RR"? and C € R?*”" (n > p > q) are constant 
matrices. The matrices C and D are assumed to be of full rank. The known non- 
linear term G(x, u) is Lipschitz with respect to x uniformly for u ¢ U, where U 
is an admissible control set. The bounded unknown function f(y, u,t) € R4 rep- 
resents the actuator fault which needs to be estimated and the uncertain nonlinear 
term W(t, x, wu) represents the modeling uncertainties and disturbances affecting the 
system. 
First, some assumptions will be imposed on the system (3.77). 


Assumption 3.1. rank(C[E, D]) = rank({E,D]) = q <p. 
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Assumption 3.2. The invariant zeros of the matrix triple (A, [E, D], C) lie in the 
left half plane. 


Assumption 3.3. The function f(t, 1) and its time derivative are unknown but 
bounded: 


IfOwOl<oy.un, IFO. dl < aa.ut, (78) 
where a(y, uv, t) and ag(y, u, t) are two known functions. 


Assumption 3.4. The nonlinear term W(t, x, w) and its time derivative are unknown 
but bounded: ; 
IMG x, wll< 6, IWGx, wll < Ga, (3.79) 


where (3 and (3, are known positive scalars. 


Under Assumption 3.1, there exists a coordinate system in which the triple 
(A, [E, D], C) has the following structure 


Aj A2 O(n—p)xr O(n—p)xq 
fle va E D> , [9px (n=p) C2| , (3.80) 


where A, € R“~?)*"—P), Cy € R?*? is nonsingular and 


_ | Op-@xr _ | %-axa@ 
EF — Eo ? Dz 4 Dy ? (3.81) 


with E>) € R@*" and Do € R4*? of full rank. 
Under Assumption 3.2, there exists a matrix L € R“-P)*P with the form 


L= [Ly Oe—sixa|, (3.82) 


with L; € R®-?)*@-® such that A, + LA; is Hurwitz. 
Without loss of generality, the system (3.77) has the form 


Xy = Ayx, + Agx2 + Gi (x, uv), 
Xo = Azxy + Agx2 + Go(x, u) + EnV (t,x, u) + Dof(y,u,t), 
y = C2x, (3.83) 


where x := col(x;, X2), x} € R"-?, x. € R?, Gi(x,u) and Go(x, uv) are the first 
n — p and the last p components of G(x, u), respectively. The matrices Ey and D2 
are given in Eq. (3.81). 
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3.4.1 First-Order SMO-Based FDI 


Consider the system (3.83), there exists a coordinate transformation z = Tx 


ge 
| 5 ale (3.84) 


where L is defined in (3.82). Thus, in the new coordinate z, the system (3.83) has 
the following form 


Zl _ Fiz + FyZ2 + (feo L] G(T'z, u), 
Zo = Agzi + F322 + Go(T'z, u) + ExW(t, T~'z, u) + Do f(y, u, 2), 
y = C222, (3.85) 


where z := col(z, Z2), z} € R" ”, z € R? and 
Fi =A,;4+ LA3, Fo = Ap+ LA,—F\L, F3 = Aq — LAs, (3.86) 


with the matrix F, is Hurwitz. 
Consider the following dynamical observer for the system (3.85) 


4 = Fit PoCy'y + [hp L] G(T, w), 
By = Agi + Fado — K(y — Co8o) + Ga(T!2, u) + vy, 9, 2, D), 
y = Cra, (3.87) 


where Z := col(Z;, Cy 'y) and ¥ is the output of the observer system. The gain matrix 
K is chosen such that 


F := ()F3Cy'!+ CK, (3.88) 


is asymmetric negative definite matrix given that C2 is nonsingular. The output error 
injection term v(y, 9, Z, u, t) is defined by 
ee 


v := ky, Z,u,t)Cy ly—3 


, ify—y # 0, (3.89) 


where k(y, Z, u, t) is a positive scalar function to be determined later. 
Let e; = z;—2Z, andey = y—¥Y = Cy(Z2 — 22). Then, the error dynamical 
system is described by 


é) = Fie, + [In—p L] (G(T 'z, u) — G(T '2, w)) (3.90) 
éy = CrA3e1 + Fey + Co (G2(T'z, uw) — G2(T~'2, u)) 
+ CoE .W(t, T7'z, u) + CrDo f(y, u, t) — Cov. Gon 
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Remark 3.2. The gain matrix K is introduced to guarantee that the following matrix 
F, O 
C,A3 F 
(3.86) and (3.88) are both Hurwitz. 


is stable. It can be directly obtained since the matrices F; and F in Eqs. 


Proposition 3.3 [36]. Suppose that Assumptions (3.1), (3.2), and (3.3) hold and the 
following matrix inequality 


stats pac i gat 2 
A’P’+PA+-—PP +e Gl_p+oP <0 (3.92) 
E 
is solvable for P where 
P:=P[I,pL], A = [Ai As)" (3.93) 
for P > 0, € and @ are some positive constants. yg is the Lipschitz constant for 


G(x, u) with respect to x. Then, the error dynamical system (3.90) is asymptotically 
stable, i.e. 


llev@ll <M le(Olen*, (3.94) 
where e, (0) is the initial error and M := ae 


Remark 3.3. The inequality (3.92) can be transformed into a standard LMI feasibility 
problem: for a given scalar 9 > 0, find matrices P, Y and a scalar € such that: 


PA, +A{P+YA3+AlY¥'+o0P+ey¥_ P Y 


P —el,_p» 0 <0, (3.95) 
id 0 -el, 
where Y := PL with P > 0. Moreover, a convex eigenvalue optimization problem 


can be posed which is to maximize yg by determining the values of P, Y and ¢ [36]. 


Proposition 3.4 [36]. Suppose that Assumptions (3.1), (3.2), and (3.3) hold and the 
gain k(-) is chosen to satisfy 


K(y, Z,u,t) = (l|C2Aall + IlCall ye) @@ + [C2 all 8 


3.96 
+ ||CyDall a(y, u, t) +, 78) 


where 1 is a positive constant and @(t) has the following dynamics 


: 1% F 
w(t) = 5 owl); w(0) = M |le:(0)|| . (3.97) 
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Then, the system (3.91) is driven to the sliding surface 
S= {(e1, ey ley = o} ; (3.98) 


in finite time and remains on it thereafter. 


It follows from Propositions 3.3 and 3.4 that system (3.87) is an SMO of system 
(3.85), where 9 is the observer output which will be used in the FDI. The proposed 
SMO will be analyzed in order to reconstruct or estimate the fault signal f(y, u, ft) 
in the presence of the uncertainty W(t, x, u). 


Fault Reconstruction Via First-Order SMO 


During the sliding motion, e, = é, = 0. Since C2 is nonsingular, it follows from Eq. 
(3.91) that 


A3e, + (G2(T~!z, u) — G2(T~'2, u)) 


+ E,V(t, T~!z, U) + Dorf, u, t) Veg = 0, (3.99) 


where vy, is the equivalent output error injection signal, obtained from a low pass 
filter. 
Since lim,_, 5. e1(t) = 0 and 
dim || G2(T'z, w) — G2(T~'2, u) | < 76 lle ll > 0. 
Thus, Eq. (3.99) can be written as 
Dz f(y, Ut) = Veq — EqW(t, T'z,u) +0), (3.100) 


where lim,_, 50 d)(t) = 0. In the case when W(t, T~!z, wu) = 0, the estimate of fault 
signal is 


FO) = Dp veq, (3.101) 


where Dz is the pseudo-inverse of Do, i.e. a Dy» = I, and 
lim (fo ~ f(yu, 1)) ei 
too 


In the case when W(t, T~'z, uw) 4 0, multiply both sides of Eq. (3.100) by Df, 
it follows 


f(t) = D} vq — Dj EyW(t, T~'z,u) + Did (0). (3.102) 
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In view of Assumption 3.4 and Eq. (3.102), it follows 


|fw =F; u,1)| < ||DIE,| 6+ |DTao], (3.103) 


where lim,-_, 5, D3 d) (t) = 0. 

The objective here is to choose an appropriate matrix DJ such that the effect of the 
uncertainty W(t, T~'z, u) is minimized, i.e. min | D} £5). Let D = {x e€ RI? | 
XD, = I, }. The set D can be parameterized as 


= {(DED2)' DI +pD" |we Reve] (3.104) 


where D3’ € R'?-”*? spans the null-space of D2 which implies that D!Y D2 = 0. 
Thus, for any Dy € D, it follows 


Dy Ey = (D1 Dy) DEE, + DY Ep. (3.105) 
The objective is transformed into the following optimization problem 


min = {| (DTD) ' DIE, + uDN E, || (3.106) 


eR (Pa) 
This can be easily solved using LMI optimization approach [40]. 


Remark 3.4. Inthe case when W(t, T~!z, uw) = 0, detection is inherent since precise 
reconstruction is achieved. However, when precise reconstruction is unavailable, 
detection is more difficult since the presence of uncertainty will make the equivalent 
output error injection be nonzero. Thus, it is difficult to distinguish the fault from the 
uncertainty. Provided the size of the bound | Dy E, | Gis relatively small compared 
to the size of the fault, a reasonable solution is to set appropriate thresholds, and a 
level of detection can still be achieved. 


3.4.2 SOSM Observer-Based FDI 


In the above section, a standard first-order sliding mode algorithm was employed 
to force e, and é, to zero exactly in finite time. Then, the fault reconstruction is 
achieved by analyzing the output error injection which was generated by low pass 
filters. Therefore, it is an approximate method [6, 32]. However, there are two main 
drawbacks in the above method. The first one is that the approximation of the equiv- 
alent injections by low pass filters will typically introduce some delays, that lead to 
inaccurate estimates or to instability for high-order systems. The second one is that a 
priori information of the bound of the fault signal is not always available in practice. 
For the first problem, the solution is to replace the discontinuous injection signal by 
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a continuous one. The solution of the second problem is to design an adaptive law for 
the gains of the SOSM algorithm. The gains of the algorithm stop increasing when 
the observation error converges to zero exactly. 

In this section, state estimation, parameter identification and fault reconstruction 
are studied for a class of nonlinear systems with uncertain parameters, simultane- 
ously. Adaptive-gain SOSM algorithms are employed handle the uncertainty with the 
unknown boundary by dynamically adapting their parameters. The approach involves 
a simple adaptive update law and the proposed adaptive-gain SOSM observer. The 
adaptive law is derived via the so-called “time scaling” approach [12], which are 
adapted dynamically according to the observation error. The uncertain parameters 
are estimated and then injected into an adaptive-gain SOSM observer, which main- 
tains a sliding motion even in the presence of fault signals. Finally, once the sliding 
motion is achieved, the equivalent output error injection can be obtained directly and 
the fault signals are reconstructed based on this information. 


3.4.2.1 Adaptive-Gain SOSML Design for FDI 
Consider the following nonlinear system 


xX = Ax+ g(x, u)+ o(y, u)O+ w(y, u) f(t), 
y=Cx, (3.107) 


A, Az 
A3 Aq 
the control input which is assumed to be known, y € VY C R? is the output vector. The 
function g(x, uv) € R” is Lipschitz continuous, é(y, u) € R’*? and w(y, uv) € R"*" 
are assumed to be some smooth and bounded functions with p > q + r. The unknown 
parameter vector 6 € R? is assumed to be constant and f(t) € R” is a smooth fault 
signal vector, which satisfies 


where the matrix A = |. € R" is the system state vector, u(t) € U C R” is 


IfOl <p, IFO < pr, (3.108) 


where /;, p2 are some positive constants that might be known or unknown. 
Assume that (A, C) is an observable pair, and there exists a linear coordinate 


I i 
P ? x= [zt zl’, with z} € R? andz. € 


transformation z = Tx = 
oa In—p 


R”~?”, such that 


Ai A12 

Ax: A22 
Hurwitz stable. 

ecT!= (Zn 0], where J, € R?*? is an identity matrix. 


e TAT'= | where the matrix Ax. = Ay — HA> € R“-?)*@-P) is 
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Assumption 3.5. There exists a function w;(y, u) such that 


Tu(y,u) = ies li (3.109) 


where wi (y, uw) € R?*’. 


Remark 3.5. Assumption 3.5 is a structural constraint on the fault distribution w(-, -). 
It means that the faults only affect on the system output channel. It should be noted 
that there are no such structural constraints on the uncertain parameters distribution 


System (3.107) is described by the following equations in the new coordinate 
system, 


2= TAT '24+Tg(T'z, u) + TH, uO + Tw(y, u) f (0), 
y = CT !z. (3.110) 


By reordering the state variables, the system (3.110) can be rewritten as 


y =Anyt+ Agiz2+ gi(Z2, y,4u) + O1(y, WO + a1 (y, uw) fC), 
Z2 = Agrza + Ary + g2(Z2, y, u) + b2(y, u)O, 
y=u, G.111) 


where 


_ | a@Q,u) = | ga, y, uw) 
roc.) = [90%], Tg(T zu) = ok (3.112) 


o\(,:): R? x R” > R’, ¢2(,-): R? x R" > R”?, giG,-,-): R? x R"? x 
R” — R? and g(.,-,-): R? x R™? x R” — RR’. 

We now consider the problem of an adaptive SOSM observer for system (3.111), 
in which the uncertain parameter is estimated with the help of an adaptive law. Then, 
an SOSM observer with gain adaptation is developed using the estimated parameter. 
Finally, based on the adaptive SOSM observer, a fault reconstruction method which 
can be implemented online is proposed. The basic assumption on the system (3.111) 
is as follows: 


Assumption 3.6. There exists a nonsingular matrix T € R®*?, such that 
T _ P| (y, u) 04 xr 
T [¢10, 4) wi(y, w)] = in (3.113) 


where ®)(y,u) € R?%47, &o(y, uv) € R"*” are both nonsingular matrices and 
bounded in (y,u)€ VY xU. 
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Remark 3.6. The main limitation in Assumption 3.6 is that the matrix [1 (y, 4), 
wily, u)| must be block-diagonalizable by elementary row transformations [37]. For 
the sake of simplicity, the case of only one fault signal and one uncertain parameter 
is considered (q =r = 1). 


Let z, = T y, where T is defined in Assumption 3.6. Then, the system in (3.111) 
can be described by 


: = = = Pi (y, u) 0 
Zy = TAyyt+TArz2+ 7TH, 22,4) + 0 O+ Pe | f(t), 


22 = Aarz2 + Ary + g2(Z2, y, u) + b2(y, u)O, 


va Ty, (3.114) 
where 
= Ay = Ay, = W, (y, 22> u) 
T-A = = ¥ T-A — 7 ; T. : : —_ A . 
ul | 12 a gi(y, 22, U) Wess 200) 


T Joss 
Let us define z, = kare Z| , where zy, € IR’, zy, € R’. Then, in view of (3.114), 
we can obtain 


Zy, = Any + Ariz2 + Wy, (y, 22, u) + Di(y, WO, 
Zy = Ay + Azz. + Wo (y, 22, u) + r(y, u) f(t), 
Z2 = Ary + ArzZ2 + go(y, 22, U) + Gr(y, u)O, 
y=T"[zy, z] - (3.115) 


Regarding the above system (3.115), the adaptive SOSM observer is represented by 
the following dynamical system 


2y, = Any + Aoi2o + Wo, (, 22,4) + O1(y, WO + pley,), 
2y, = Any + An22 + We, (y, 22, u) + U(ey,), 
2p = Any + Annda + gry, 22, u) + dry, u)4, (3.116) 


where ju(-) is the SOSM algorithm 


Ls) = A(t) |s|2 sign(s) + ait | sign(s)dt +k)(t)s + kat f sdT, 
0 0 


and the adaptive gains X(t), a(t), k\(t) and k,,(t) are determined later. 

The observation errors are defined as e,, = zy, — Zy,, Cy, = Zy, — 2), €2 = Z2—- 
2, 0=0- 6. The estimate of 9, denoted by 6, is given by the following adaptive 
law 
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6 = —K(y, u) (Any + Agi 2. + Wo (y, 22, u) + Pi(y, WO — in) » Bt) 


where K (y, u) is a positive design matrix which will be determined later. 


Remark 3.7. It can be seen that the adaptive law (3.117) depends upon Z,,. A real- 
time robust exact differentiator proposed in [17] can be used to estimate the time 
derivative of z,, in finite time. The differentiator has the following form 


: I al 
Zo = —AoLg |Z0 — Zy, |? Sign(zo — Zy,) + 21, 
Z1 = —aoLosign(Zo — Zy,), (3.118) 


where Zo and z, are the real-time estimations of z,, and Z,, respectively. The param- 
eters of the differentiator Ay = 1, a9 = 1.1 are suggested in [17]. Lo is the only 
parameter needs to be tuned according to the condition |Z,,| < Lo. 


Subtracting (3.116) from (3.115), the error dynamical equation is described by 


by = Ane, + Holy, 22, 22, U) + d2(y, WO, (3.119) 
8 = —K(y,u) (Arier + Woy(y, 22, 22,0) + Pi(y, 8), GB.120) 
by, = —p(ey,) + Anier + Oi (y, u)O + Wo, (y, 22, 22, U), (3.121) 


éy, = —p(ey,) + Aner + D2(y, uw) f(t) + Wy (y, 22, 22,4), (3.122) 


where 


92(ys 22, 22, U) = ga(y, 22, 4) — gal, 22, 4) 
Wy, Z2, 22, U) = Way, 22, U) — Ws, (y, 22, u) 
Wa, Z2, 22, U) = Wa, Z2,U) = W,,(y, 22, U) 


Some assumptions are imposed upon the error dynamical systems (3.119)— 
(3.122). 


Assumption 3.7. The known nonlinear terms go(y,2Z2,u), Wg,(y, Z2,u) and 
W,(y, Z2, u) are Lipschitz continuous with respect to z2 


Ilg2(y, Z2, 4) — gay, 22, wll < yallz2 — Zall, (3.123) 
| Wo, (y, 225 u) _ Wi, (y, Zo, u)|| < Vo I|Z2 ~ Zoll, (3.124) 
|| Woo (Ys 22, U) — Wo (y, 22, WII < Yoollz2 — Zell, (3.125) 


where Y, ,Y, and 2 are the known Lipschitz constants of W,, (y, Z2, 4), Wg, (y, Z2, U) 
and go(y, Z2, u), respectively [39]. 
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Assumption 3.8. Assume that the Hurwitz matrix A 9 satisfies the following Riccati 
equation 


Ax Pi + PiAr. + PL Pi + 2+ e)In—p = 0, (3.126) 


which has a SPD solution P; for some ¢ > 0 [24]. 


Assumption 3.9. Assume that the positive design matrix K (y, u) satisfies the fol- 
lowing equation 


K(y, u)®\(y, uv) + Of (y, u)K"(y, w) — 95, K(y, u) K"(y, u) — lg = 0, (3.127) 


for some ¢€ > 0. 


Assumption 3.10. Itis assumed that ||®,(y, )||, || ®2(y, w)|| are bounded in (y, uw) € 
Yxu. 


Now, we will first consider the stability of the error systems (3.119) and (3.120). 


Theorem 3.2. Consider systems (3.119) and (3.120) satisfying Assumptions (3.7)- 
(3.10). Then, the error systems (3.119, 3.120) are exponentially stable, if for any 
(y,u) € Y x U, the following matrix 


01= EIn—p Pida(y, uw) — AN KT, w) 


o3(y,u)P; — K(y, u)An él, 5 G28) 


is positive definite. 


Proof. A candidate Lyapunov function is chosen as 
Vi (e2, 0) = e} Pier + 678, (3.129) 


and the time derivative of V; along the solution of the system (3.119) and (3.120) is 
given by 
Vi = €3 (As Py + Pi Anz)er — 2e3 A5, K'(y, u)0 + 2€3 Pido(y, u)O 
+2¢3 Pi Go(y, 22, 22, u) — 26°K (y, u) We, (y, 225 22, ) 
—0" (K(y, wWP1(y,u) + OT (y, WK" (yw) 6 
< €3 (Ax Pi + PiAn + 3 Pi Pi + 2In—p)er 
+2e; (Pida(y, u) — A}, K'(y, u)) 6— "6 


= [et BT] 0, Ht 


Hence, the conclusion follows from the assumption that Q, given in (3.128) is pos- 
itive definitive in (y,u) € Y x U. 
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Remark 3.8. Theorem 3.2 shows that lim e2(t) = 0 and lim A(t) = 0. Conse- 
t>oo t> ow 


quently, the errors e2, 6 and its derivatives é7, 0 are bounded. Under Assumptions 
(3.7) and (3.10), the time derivatives of the nonlinear terms in the error dynamics 
(3.121) and (3.122) are bounded: 


a A v Aa 
| .F (Ane. + @y(y, u) f(t) + We (y, 22, 22, »)| 


A 
o 


(3.130) 


A 
> 
M 


d j= as x 
| ai (Aaie2 + Bi (y, u)O + Wo, (y, 22, 22; u)) | = 


where x; and yz are some unknown positive constants. 


According to the results of Theorem 3.1, the finite time convergence of the systems 
(3.121) and (3.122) can be obtained directly. Theorems (3.1) and (3.2) have shown 
that systems (3.116) and (3.117) are an asymptotic state observer and an uncertain 
parameter observer for the system (3.115), respectively. In the next part, we will 
develop the fault reconstruction approach based on those two observers. 


3.4.2.2 Fault Reconstruction via SOSM Observer 

The fault signal f(t) will be reconstructed based on the proposed observer by using 
an equivalent output error injection which can be obtained once the sliding surface 
is reached and maintained on it thereafter. 


It follows from Theorem 3.1 that e,, and é,, in (3.122) are driven to zero in finite 
time. Thus, the equivalent output error injection can be obtained directly 


[u(ey,) = Arrer + Go(y, u) f(t) + Wo (¥, 22, 22, u). (3.131) 


From Assumption 3.6, ®2(y, uv) is a bounded nonsingular matrix in (y,u) € Y x U 
and lim;-. oo Wy, (y, 22, 22, u) = O, then the estimate of f(t) can be constructed as 


f@ = OF! (y, Wuey,). (3.132) 


Theorem 3.3. Suppose that conditions of Theorems (3.1) and (3.2) are satisfied, 
then f (t) defined in (3.132) is a precise reconstruction of the fault signal f (t) since 


tim If) — f@Il = 0. (3.133) 
Proof. It follows from (3.131) and (3.132) that 


IFO) — FOI = Oy! y, w(Are2 + Wo.) 


x _ 7 (3.134) 
< 12 (y, w)Aaallllezll + Yello Cy, w)lllleall. 
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From Theorem 3.2, lim |le2|| = 0, it follows that 
too 


tim |f@ — fll = 0. (3.135) 


Hence, Theorem 3.3 is proven. a 


3.5 Illustrative Examples 


Let us now see some applicative examples of SOSM observer design and its appli- 
cation to FDI. 

In this example, we shall briefly outline the application of observers described in 
this chapter to a class of nonlinear systems. Consider the following nonlinear system 
[18], 

dy = xy — 37 — me +393 +00), 
Xp = —X, — 3x7 + Die + 2x,xX2 + 2x1x3 — 2x, w(t), (3.136) 


P 2, 2 
X3 = Xp +X — XQ — XB — 3, 


y'=[y1 y2] = [a1 23], 


where x? = [x1 x2 x3] are the states, w(t) = 2 - sin(5t) is an unknown disturbance, 
and y are the measurable outputs. Define the following coordinate transformation, 


xX] 
z= @(x) = |xi+x|. (3.137) 


X3 


The system (3.136) can be transformed into the following system 


Z= 21-2423 4+ w(t), 

Zo = —z1 + 42123 — Szt, 

23= 2) -2—23+20—Zz, 

y' = [z1 zs]. (3.138) 


The parameters for the simulation are shown in Table 3.1. 


Table 3.1 Parameters for the numerical simulation 


Original states Values Estimate states Values 
x1 (0) 0 41 (0) 0.3 
x2 (0) 0 X2 (0) 0.39 


x3(0) 0 X3(0) 0.3 
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3.5.1 Adaptive-Gain SOSM Observer Design 


The adaptive-gain SOSM observer is designed as follows 


Z=—ytyt+vy,), 
2. = —y1 + 4yiy2 — Sy}, + Biv (Za — 20), 


yi + 2yp — 2 — yg FYE), (3.139) 


NX 
w 
II 


where ey, = y) — 2, ey, = y2 — 23, the error injection term v(-) is obtained from 
the adaptive gain SOSM algorithm (3.60), (3.62), and (3.63) and 


i. at hey 26, 


%=-vV(ey,), E, = 3.140 
2 (y.) : 0, otherwise, ( ) 
where € is a small positive constant. 
Denote e2 = z2 — 22, the error dynamical system is given as follows: 
éy, = —V(ey,) —Z2+u(t), 
ey, = “UC y,) — 22, 
éy = —E\v3(Z2 — 20). (3.141) 


With appropriate gains of the adaptive gain SOSM algorithm, the error dynamical 
system (3.141) converges to zero in finite time. Moreover, the estimate of the distur- 
bance is obtained as W(t) = v(ey,) — v(ey,). The state observation and disturbance 
estimation is shown in Fig.3.3. It can be seen from Fig.3.4 that the adaptive law 
(3.63) is effective under disturbance effect. 


20 2.5 


— v(t) 


== = a(t) 


State estimate 


0 05 1 15 2 
Time (s) 
(a) Estimate of x2 


Fig. 3.3. Estimates of state and disturbance 
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Adaptive gain 


1 
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Time (s) 


Fig. 3.4 Adaptive law L(t) of the SOSM algorithm (3.63) 


3.5.2. SOSM Observer Based FDI 


An example is given in order to illustrate the design of adaptive-gain SOSM observer- 
based fault reconstruction scheme. Consider a pendulum system [5] 


xX, =x. + f(t), 

; 1 g GD V, 

Xp = -u — — sin(x,) — —x, 

2 7 a 1 7 2 

oe (3.142) 


where x1, x2 are the angle of oscillation and angular velocity, respectively. M = 1.1 
kg is the pendulum mass, g = 9.815 m/s? is the gravitational force, L = 1 m is the 
pendulum length, J = ML? = 1.1 kg- m* is the arm inertia, V, = 1.8 kg- m/s? is 
the pendulum viscous friction coefficient and f(t) is considered as a bounded fault 
signal. For simulation purposes, it was taken as 

f@) = 0.5 sin(2t) + 0.5 cos(5r). (3.143) 
The controller u is known which is used to drive x; to the desired value x{ = sin(t) 


u = —30sign(s) — 25sign(s), (3.144) 


where s = x; — xf. 
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The observer based on STA or SOSM algorithm is designed as 


Ky) = Xp + wley), 


Li ! date je 3 
= —u — = sin(x) — —Xo, 
J OL eG 
(3.145) 


=> 
ll 
=> 


where ey = y — ¥. Let en = x) — Xo, the error dynamics are given as follows 


ey = —M(ey) + er + f(t), 
; Vs 
& = — Te. (3.146) 


During the sliding motion e, = é, = 0, a fault reconstruction signal a (t) = wey) 
is introduced given that lim;-,. e2(t) = 0. 
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Fig. 3.5 The performance of observer (3.145) 
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Fig. 3.6 Fault reconstruction and its error for the pendulum system 
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The initial values x, (0) = —1, x2(0) = 3 are set for the pendulum and x, (0) = 
X2(0) = Oare set for the two observers. The state estimations and fault reconstruction 
are shown in Figs. 3.5 and 3.6, respectively. 


3.6 Conclusion 


In this chapter, we have discussed the SMO design and its application to fault recon- 
struction. Traditional first order SMO design for linear uncertain systems has been 
introduced first and their drawbacks have been discussed. Then, SOSM algorithms 
have been developed to generate continuous output error injection signals, which 
eliminate the use of any low pass filters and their applications to fault reconstruction 
have also been presented. Finally, our contributions in adaptive-gain SOSM observer 
and its use in FDI have been introduced. An adaptive-gain SOSM algorithm has 
been proposed, which does not require the a priori knowledge of the uncertainty’s 
boundary. Additionally, adaptive-gain SOSM observer based fault reconstruction has 
been considered for a class of nonlinear systems with uncertain parameters. State 
estimation, parameter identification and fault reconstruction are performed, simul- 
taneously. Finally, two illustrative examples have been provided to show both the 
observer designs and their applications to fault reconstruction. 

In the subsequent chapters, we shall concentrate on the applicative aspect of the 
above discussed sliding mode algorithms. Namely, a complex nonlinear system, 1.e., 
PEMFC will be presented and the observation and FDI of its subsystems will be 
discussed. 
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Chapter 4 ®) 
Introduction of Proton Exchange crest 
Membrane Fuel Cell Systems 


This chapter introduces some fuel cells classified into five different categories based 
on the electrolyte chemistry. The power generation mechanism of the PEMFC is 
presented. The PEMFC air-feed system model including supply manifold model, 
compressor model, cathode flow model, etc. is detailed. For an understandable SMC 
scenario design to the PEMFC, an experimental validation performed through the 
HIL test bench is provided, based on a commercial twin screw compressor and a 
real-time PEMFC emulator. 


4.1 Introduction of PEMFCs 


Fuel cells are electrochemical devices that convert the chemical energy of a gaseous 
fuel directly into electricity. They are under intensive development in the past few 
years as they are regarded as an efficient carbon free electricity production technol- 
ogy. The basic principle of the fuel cell was first discovered by Grove [3] in 1839, 
through an experiment that demonstrated that the reaction of hydrogen and oxygen 
produces electrical current. Modern fuel cells consist of an electrolyte sandwiched 
between two electrodes as shown in Fig. 4.1. 

Fuel cells can be classified into five different categories based on the electrolyte 
chemistry: PEMFC, solid oxide fuel cell (SOFC), molten carbonate fuel cell (MCFC), 
phosphoric acid fuel cell (PAFC) and aqueous alkaline fuel cell (AFC). Among these 
kinds of fuel cells, PEMFCs are suitable for automobile applications due to high 
energy density, low working temperature, simple structure and long life as well as 
low corrosion [11]. Table4.1 provides a summary of various fuel cell types and 
corresponding characteristics. 

The fuel cell under consideration throughout this monograph is the PEMFC. The 
electrolyte membrane of a PEMFC has a special property that allows only positive 
protons to pass through while blocking electrons. Hydrogen molecules are split into 
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Fig. 4.1 Fuel cell reaction 


Table 4.1 Summary of typical fuel cell characteristics (Table 16, [5]) 


Electrolyte Operating Anticipated Comments 

material temperature applications 

PEMFC 80°C Stationary and | Minimum contamination and material 

automobile problem 

AFC Approx 100 °C | Space program | Susceptible to contamination, very 
expensive 

PAFC Approx 100 °C | Stationary Higher temperature and longer warm up 
time makes unsuitable for vehicles 

SOFC 1000 °C Stationary Very high temperature create material 
problems, steam generation could 
increase efficiency by cogeneration 

MCFC 600 °C Stationary Same as SOFC 


protons and free electrons at the anode. These protons flow to the cathode through 
the electrolyte and react with the supplied oxygen and return electrons to produce 
water. During this process, the electrons pass through an external load circuit and 
provide electricity. There are two important subsystems required for proper operation 
of PEMFCs. First is the air-feed system that supplies oxygen to the cathode and 
therefore indirectly regulates the net power output of the fuel cell. Then, there is the 
power converter system that forms a link between the fuel cell output and the power 
bus of the power system. This chapter is dedicated to the description and modeling 
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of the air-feed system. We will start by a short description of the power generation 


in the PEMFC and highlight the importance of air-feed system in it. 


4.2 PEMEFC Stack Voltage 


In order to understand the importance of the air-feed system, let us look at the 
power generation mechanism of the PEMFC. A typical PEMFC polarization curve 
is shown in Fig. 4.2. The cell voltage is modeled from its static characteristic, which 
is a function of stack current, cathode pressure, reactant partial pressures, fuel cell 
temperature and membrane humidity [11] 


Ute = E — vVact — Vohm — Ucone: 


The open circuit voltage E can be calculated as 


E = 1.229 — 0.85 - 107? (Tye — 298.15) 


+ 4.3085 - 10-°Ty¢ E (pu) + sin (ro) 


(4.1) 


(4.2) 


where Ty, is the temperature of the FC (in Kelvin), py, and po, are the par- 
tial pressures of hydrogen and oxygen, respectively (in bar). It can be seen that 
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Fig. 4.2. Typical fuel cell voltage 
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the air-feed system indirectly controls the output voltage and hence the power, as it 
controls the oxygen pressure in the cathode. The current density i is defined as 


Isp 


i=, 
Afe 


(4.3) 


where /,,(A) is the stack current and A fe(cm’) is the active area. The activation loss 
Vacr, Ohmic loss Voym and concentration loss Ueone are expressed as follows: 


1. Vacr = aln (4) is due to the difference between the velocity of the reactions 


in the anode and cathode [6], a and ig are constants which can be determined 
empirically. It should be noted that this equation is only valid for i > ig. There- 
fore, a similar function that is valid for the entire range of i is preferred: 
Vact = Vo + Va (i — el 1) vo (volts) is the voltage drop at zero current den- 
sity, and v, (volts) and b; are constants that depend on the temperature and 
the oxygen partial pressure [1]. The values of vo, vg and b; can be determined 
from a nonlinear regression of experimental data. The activation loss is shown in 
Fig. 4.3(a). 
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Fig. 4.3. Voltage drops due to different types of losses in FC. a Activation losses; b Ohmic losses; 
c Concentration losses; d Total losses 
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2. Uohm = 1Ronm 1s due to the electrical resistance of the electrodes, and the resis- 
tance to the flow of ions through the electrolyte [6]. Ronm (Q - cm?) represents 
the fuel cell internal electrical resistance. The ohmic loss is shown in Fig. 4.3(b). 


. \ ba 
3: -Vigae= (d3 - ) results from the drop in concentration of the reactants due 


Umax 


to the consumption in the reaction. b3, b4 and ig, are constants that depend on 
the temperature and the reactant partial pressures. i,,¢, is the current density that 
generates the abrupt voltage drop. The concentration loss is shown in Fig. 4.3(c). 


The stack voltage V;,; is calculated as the sum of the individual cell voltages, 
Vor = NVfe, (4.4) 
where n is the number of cells. Thus, the fuel cell power is calculated as 


Pst = Ts Vt : (4.5) 


4.3 PEMEFC Air-Feed System Model 


Itis clear from the above discussion that the power generation inside the PEMFC core 
depends upon several variables and factors, such as oxygen and hydrogen pressures, 
temperature and various physical parameters of the fuel cell. All these must be con- 
trolled rigorously for proper and safe operation of the fuel cell. Hence an operational 
PEMEFC system contains the fuel cell core and a certain number of auxiliary systems 
for control (Fig. 4.4). 

The aim of the air-feed system is to regulate the oxygen quantity in the cathode. It 
usually consists of an electromechanical air compressor that maintains the required 
oxygen pressure and mass-flow in the cathode of PEMFC. As mentioned in the 
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Fig. 4.4 Fuel cell system scheme 
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introduction, this system can consume up to 30% of the fuel cell power and requires 
precise control in order to optimize the net power output of the PEMFC system. 
This system has been modeled under the following assumptions [11, 12]: 


All gases obey the ideal gas law; 

The temperature of the air inside the cathode is equal to the stack temperature, also 
is equal to the temperature of the coolant exiting the stack. The stack temperature 
is well controlled; 

The input reactant flows are humidified in a consistent and rapid way and the high 
pressure compressed hydrogen is available; 

The water inside the cathode is only in vapor phase and any extra water in liquid 
phase is removed from the channels; 

The flooding of the gas diffusion layer is neglected; 

The spatial variations are neglected, it is assumed that the flow channel and the 
gas diffusion layer are lumped into one volume; 

The anode pressure is well controlled to follow the cathode pressure. 


4.3.1 Supply Manifold Model 


For the supply manifold, the inlet and outlet mass flows are the compressor flow 
W-.p and supply manifold flow Wm our, respectively. The supply manifold model is 
described by the following equation 


dDsm _ RaTep,out 
dt Vs 


(Wep W nein) , (4.6) 


where V,,, is the supply manifold volume and T,, ,,; is the temperature of the air 
leaving the compressor which is calculated as follows: 


q-1 
T, 2) a 
Tep,out = Tym + an ( = ) 7 ’ (4.7) 
Nep Patm 


where 7),p is the compressor efficiency (its maximum value is 80%). 

The relationship between the flow and the pressure drop can be simplified as a 
linear nozzle equation since the pressure difference between the supply manifold 
Psm and the cathode pq is small: 


Wsm,out = Ksm,out (Psm ~ Pca) ’ (4.8) 


where Ksn,our is the supply manifold outlet flow constant. 
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4.3.2 Compressor Model 


The air feed system consists of a twin screw compressor and a permanent magnet 
synchronous (PMSM) motor [8]. In this part, the model of the air compressor is 
discussed, it is used to provide oxygen to the cathode side of the fuel cell system. 
The inputs to the model consist of inlet air pressure pep jn, inlet air temperature 
T-p,in and supply manifold pressure p;,,. The inlet air is typically atmospheric and 
its pressure and temperature are assumed to be Pgrm = 1 atm and Torn, = 25°C, 
respectively. 

The dynamic state in the model, i.e., compressor speed wp is given by the fol- 
lowing equation 

dWep 1 


dt = qe. (Te Tep fwep) ’ (4.9) 


where J,, is the motor and compressor inertia (kg - m7), Tem is the compressor motor 
torque input (N- m), 7;p is the torque of the compressor (N - m) and f is the friction 
coefficient. 
Tom = "Nem k; I, > 
q-1 


Colin on \ oT (4.10) 
to = a (2 ) 1 We, 


NepWep Patm 


where 1. is the motor mechanical efficiency, k, is the motor constant, J, is the motor 
quadratic current, C, is the specific heat capacity of air (1004 J- kg"!. K7), vis 
the ratio of the specific heats of air (y = 1.4) and W,, is the compressor mass flow 
rate. 

The mass flow rate of the twin screw compressor W,, depends on its angular speed 
Wep, Which is independent of the supply manifold pressure p,,, [13]. The following 
relation is given as 


1 
Wep = 5, he Vepr/tr Pawep: (4.11) 


Table 4.2 Parameters used in compressor model 


Symbols Parameters Values 

Pa Air density 1.23kg/m? 

k; Motor constant 0.31 N m/A 

f Motor friction 0.00136 V/(rad/s) 

Jep Compressor inertia 671.9 x 1075 kg m? 
Nep Compressor efficiency 80% 

Nem Motor mechanical efficiency 98% 

Vepr/tr Compressor volume per turn 5 x 1074 m3 /tr 

Cp Constant pressure specific heat of air 1004 J/kg K) 
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where 7),—, is the volumetric efficiency, V.p,/:- is the compressed volume per turn 
and pq is the air density. The parameters used in the compressor model are given in 
Table 4.2. 

The compressor motor power P,, is calculated as follows 


Pop = TemWep+ (4.12) 


4.3.3 Cathode Flow Model 


The thermodynamic properties and mass conservation are used to model the behavior 
of the air inside the cathode [11]. The dynamics of the oxygen, nitrogen and vapor 
partial pressures are described by the following equations 


d RT. 
=e = aor Le (Wo,,in _ Wo), out = W ox react) , 
dt Mo, Vea (4.13) 
d RT . 
BM: = esi oe (Wy, in Wny,out) ’ 
dt Mn, Vea 


where Wo,,in is the mass flow rate of oxygen entering the cathode, Wo, our is the 
mass flow rate of oxygen leaving the cathode, Wo, act iS the rate of oxygen reacted, 
Wy,.in 18 the mass flow rate of nitrogen entering the cathode, Wy, ou; is the mass 
flow rate of nitrogen leaving the cathode, Mo, is the molar mass of oxygen, My, is 
the molar mass of nitrogen and V,, is the cathode volume. 

The inlet mass flow rate of oxygen Wo, in and nitrogen Wo, our are calculated 
from the inlet cathode flow Wea in 


XO y,atm 
Wo,.in —_ Wea,in» 
1 + Watm (4 14) 
1- XO,,atm : 
Wyy,,in ———— Wea.in, 
1 + Watm 
where x9, atm 18 the oxygen mass fraction of the inlet air 
YO>,atm Mo, 
XO,,atm = (4.15) 


Yor,atmMo, ate (1 = YOr,atm) My, , 


with the oxygen molar ratio yg, atm = 0.21 for inlet air and the humidity ratio of 
inlet air wy; is defined as 


M, Po 
YOo,atmMo, F (1 = VOu-atm) My, Patm — Pv 


(4.16) 


Watm = 


where M, is the molar mass of vapor and p, is the vapor pressure. 
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The vapor partial pressure p, is determined 


Pv = Patm Psat (Tye) , (4.17) 


where dam is the relative humidify at ambient conditions (its value is set to 0.5), 
which can be used to describe the relation between the vapor partial pressure and 
the saturation pressure and Psa; (Tye) is the vapor saturation pressure. The relative 
humidify ¢ is defined as the ratio of the mole fraction of the water vapor in the 
mixture (air and water vapor) to the mole fraction of vapor in a saturated mixture at 
the same temperature and pressure. The value of relative humidify equals 1 means 
that the mixture is saturated or fully humidified. If there is more water content in the 
mixture, the extra amount of water condensed into a liquid form will be removed 
from the channels. 

The saturation pressure Psq; depends on the temperature and is calculated from 
the equation given in [9] 


logig (Psat) = — 1.69 x 107° T#, + 3.85 x 107" T}, 
— 3.39 x 10“*TF, + 0.143Ty, — 20.92, (4.18) 
where the saturation pressure ps; is in kPa and the temperature Ty, is in Kelvin. 


The cathode inlet flow rate Wain is assumed to be the same as the supply manifold 
outlet flow rate Wem our (4.8) 


Wea,in = Kea,in (Psin — Pca) ’ (4.19) 


where the cathode pressure p-q is assumed to be spatially invariant, which is the sum 
of oxygen, nitrogen and vapor partial pressures, 


Pca = POr oa PN, + Psat (Tye) . (4.20) 


The following equations are used to calculate the outlet mass flow rate of oxygen 
Wo,,our and nitrogen Wy, our in (4.13) 


X0O>,ca 
Wo,,out _ ——— Wea, outs 
1+ Wea,out 
i (4.21) 
— X0O),ca 
Wyy,, out = ———— Wea. outs 
1+ Wea,out 


where x9, cas Wca,our are the oxygen mass fraction, humidity ratio inside the cathode, 
respectively. 
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7 _ Yo ,caMo, 

On,ca = ’ 

— YOr,caMo, + (1 ~ Yor,ca) My, (4 22) 
M, Psat , 

Wea,out = 


Yoo,caMo, + (i — Voucs) My, PO, a PNy 


Unlike the inlet flow, the oxygen mole fraction of the cathode outlet flow yo, ca 
is not constant since oxygen is reacted, which is calculated as 


Po 
YO3.ca = —————.. (4.23) 
Po, “F PN 
Using (4.22) and (4.23), Eq. (4.21) can be rewritten as 
Mo, Po, 
Wo,,out = OnP Or Wea,outs 
Mo, Po, + My, Pn, + My Psat 
(4.24) 
My, Pn, 


Mo, POr + My, Pn + M, Psat ca,out 


The total cathode outlet flow rate Weg our is calculated by the nozzle flow equation 
[10], 


1 


1 y-1 2 al 
rat (a)'| 2 : (42) | » tortie > (Gir) 
Wea,out = fe en) 


7 i 
CpAr pea 3 (_2_\ 7-H Pain ~ (_2_\71 
Jit; \7H for ia = (yet) 


The rate of oxygen consumption Wo, +eacr is a function of the stack current J,,, which 
is calculated using electro-chemistry principles 


—, (4.26) 


where F is the Faraday number. Oxygen excess ratio in the cathode \o, is defined as 
the ratio between the amount of oxygen supplied by the compressor and the oxygen 
reacted in the FC stack. This value is considered as a performance variable of the 
system since it determines the safety of the fuel cell. 


Wo in 
Aor 2 29 


= aE (4.27) 
Wo,,react 


4.3.4 Dynamic Model with Four States 


In view of Eqs. (4.6), (4.9), and (4.13), the nonlinear model of the fuel cell is com- 
pleted with the following four states [10] 
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x =[Po, Pv; Wep Psm|] - (4.28) 


The four state dynamic model is written as follows 


. b3x] Wea,ou 
1 = Dy (x4 — x1 — x2 — ba) — — byl, 
b4x, + bsx2 + be 
7 b3x2 Wea,ou 
Xq = bg (x4 — x1 — X2 — bp) — 
bax, + b5x2 + be 
b bia (4.29) 
X3 => box3 i ( = ) 1 Wep + by3u, 
x3 bi 


bin 
. Xx 
X4 = dig [ + Dis (5) = | x [Wop — dis (x4 — x1 — X2 — bp). 
ul 


The stack current /;, is traditionally considered as a disturbance and the control input 
u is the motor’s quadratic current. The outputs of the system are 


y= [y1 y2 ya] = [Peon Wep Ve] , (4.30) 


where the stack voltage V;,, the supply manifold pressure ps, and the compressor 
air flow rate W,, are given in Eqs. (4.4), (4.6), and (4.11), respectively. The system 
performance variables are defined as 


z= [z: z2]’ =[Prer Ao.) (4.31) 


where Py; is the fuel cell net power and Xo, is the oxygen excess ratio. 

The fuel cell net power P,,.; is the difference between the power produced by the 
stack P;, and the power consumed by the compressor. Thus, the net power can be 
expressed as 

Pret = Pst — Pops (4.32) 


where the stack power P,, and the compressor power P,, are given in Egs. (4.5) 
and (4.12). The oxygen excess ratio 9, is defined as the ratio between the oxygen 
entering the cathode Wo, i, and the oxygen reacting in the fuel cell stack Wo, react 


Wo in big (Psm —_ Pca) 
a (4.33) 
x Wo), react Diol st 


Due to the reasons of safety and high efficiency, it is typical to operate the stacks with 
this value equals 2 during step changes of current demand [11]. It should be noted 
that positive deviations of \g9, above 2 imply lower efficiency, since excess oxygen 
supplied into the cathode will cause power waste and negative deviations increase 
the probability of the starvation phenomenon. The parameters b;, i € {1,--- , 19} 
are defined in Table 4.3. 
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Table 4.3 Parameter definition 
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b1 = TEV bs My Vea T+waim 
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Te 


1 XOp.atm 
b= 3g Nv—ce Vepr/tr Pa big= Kea,in Taeain 
nMo. 
bio = 4F* 


4.4 Experimental Validation 


In this section, the experimental validation is performed through the HIL test bench 
shown in Fig. 4.5. It consists of a physical air-feed system, based on a commercial twin 
screw compressor and a real-time PEMFC emulator [7, 8]. The complete architecture 


Twin screw compressor 


Fig. 4.5 Test bench [8] 
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Fig. 4.6 Scheme of HIL system used in the experiments 


of the experimental system is shown in Fig.4.6. The test bench is controlled by 
National Instruments CompactRIO real-time controller and data acquisition system, 
using a sampling frequency of 1 kHz. 

The objective of the air-feed system is to provide sufficient quantity of oxygen 
to the PEMFC cathode, keeping the oxygen excess ratio (Ag,) at its optimum value. 
The twin screw compressor of the PEMFC air-feed system has a flow rate margin 0- 
0.1 kg/s at a maximum velocity of 12000 RPM, it consists of two helical rotors which 
are coupled directly to its motor. Air intake is at the opposite side of the mechanical 
transmission and the output pressure is regulated by a servo valve. The compressor 
is driven by a PMSM through an inverter controlled by 3-phase currents I, Ip, Ic. 
The system is controlled by a robust sub-optimal SOSM controller presented in [8]. 


4.4.1 PEMFC Emulator 


The PEMFC emulator (Fig. 4.6) has been designed using experimental data obtained 
from a 33 kW PEMFC unit containing 90 cells in series. It has been developed on the 
FPGA platform of CompactRIO running at a clock frequency of 40 MHz and an over 
all loop frequency of 10 kHz. The nominal value of the PEMFC system parameters 
obtained from the PEMFC unit (as modeled in the emulator) are given in Table 4.4. 
The emulator parameters can be varied through software, in order to perform stability 
and robustness tests. The advantage of this emulator is that it provides an economical 
and safer alternative for real-time experiments on auxiliary fuel cell system, that does 
not require excessive hydrogen consumption. 

The PEMFC emulator runs a detailed 9th order dynamic PEMFC model, com- 
prising of the following state vector 


rT 
x= [Po PH» PN; Psm; "sm, Pv,ans Pv,ca,s Prm; Tst] ’ (4.34) 
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Table 4.4 Emulated PEMFC system parameters [8] 


Symbols Parameters Values 

n Number of cells in fuel cell stack 90 

R Universal gas constant 8.314 J/(mol K) 

Ra Air gas constant 286.9 J/(kg K) 

Patm Atmospheric pressure 1.01325 x 10° Pa 
Tain Atmospheric temperature 298.15 K 

Afe Active area 800 cm? 

F Faraday constant 96485 C/mol 

Ma Air molar mass 28.9644 x 1073 kg/mol 
Mo, Oxygen molar mass 32 x 107-3 kg/mol 
My, Nitrogen molar mass 28 x 10-3 kg/mol 
Mn, Hydrogen molar mass 2 x 1073 kg/mol 

My Vapor molar mass 18.02 x 1073 kg/mol 
Cp Discharge of the nozzle 0.0038 

Ar Operating area of the nozzle 0.00138 m? 

7 Ratio of specific heats of air 1.4 

Jep Compressor inertia 6719mg- m2 

f Motor friction 1.36 mV/(rad/s) 

kt Motor constant 0.31 N m/A 

Cp Specific heat of air 1004 J/(Kg - K) 

Nem Motor mechanical efficiency 98% 

Vea Cathode volume 0.0015 m3 

Vom Supply manifold volume 0.003 m? 

Van Anode volume 0.0005 m3 

Vepr/tr Compressor volume per turn 5-104 m3 /tr 

Kea,in Cathode inlet orifice constant 0.3629 x 1079 kg/(Pa s) 
Kan,in Anode inlet orifice constant 0.21 x 1075 kg/(Pa s) 
Pa Air density 1.23kg/m3 

X0p,ca,in Oxygen mass fraction 0.23 


where py, is hydrogen pressure in the anode, ms, is mass of air in supply manifold, 
Pv.an 18 Vapor partial pressure in the anode, Py.cq is vapor partial pressure in the 
cathode, p,m is pressure of return manifold and 7;; is the stack temperature. This 
dynamic model is based on Pukrushpan et al. [11], with the added temperature model 
described by a lumped thermal model [11] 


ae Ossi = W.Cp, (Tsi _ 


: TAs 
Osou = Ist + Vact + Ts: Rohm ’ 


AF 


Tein) > 


(4.35) 
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where m,, is the heat mass of the stack, C,,, and Cp, are the specific heat, W, is the 
coolant flow rate considered as a control variable, T,.;, is the coolant temperature 
at the stack inlet and Deon is the internal energy source. The latter is calculated 
as a function of the stack current, temperature, electrical resistance of stack layers 
Ronm, Faraday’s number F and the entropy change As. The physical parameters 
were obtained through extensive experimentation. 

The PEMFC works under several safety constraints, the most important of which 
are the anode and cathode pressure difference (to be kept minimum) and stack tem- 
perature (to be regulated). These constraints require additional control in real fuel 
cells, which have been replicated in the emulator as well. An internal controller is 
implemented in the emulator to ensure that the anode pressure follows the cathode 
pressure at all times. The temperature is also controlled and its value can be set 
through software. This added characteristic is a key feature of the emulator, as com- 
pared with other works [2, 4, 11, 14], since the heat produced in real fuel cells (due to 
heat produced by irreversible energy occurring in the chemical reactions and Joules 
losses) is a major factor that causes large parametric variations in the system. 


4.4.2 Validation of 4th Order Model with PEMFC Emulator 


The experimental stack voltage, the predicted stack voltage by the 4-th order model 
and the predicted stack voltage by the PEMFC emulator, as well as the experi- 
mental data obtained from the 33 kW PEMEFC stack are shown and compared in 
Fig. 4.7. From this figure, it can be concluded that the 4-th order model replicates the 
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Fig. 4.7 Stack voltage response 
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Stack Current (A) 


Fig. 4.8 Stack current under load variation 


characteristics of the Fuel Cell with sufficient precision with a relative error less than 
2.5% (+1.6 volts). 

To demonstrate the FC model characteristics, the stack current shown in Fig. 4.8 
was varied between 100 and 450 A during the tests and a static feed-forward controller 
is used to control the compressor voltage so that the oxygen excess ratio maintains 
near 2. 

It can be seen from the Fig.4.9(b) that a step increase in the load current (i.e. t 
= 20 s) causes a drop in the oxygen excess ratio. The performance variables (Pye, 
Xo,) in Fig.4.9 and the states variables (p0,, Py, Wcps Psm) in Figs. 4.10 and 4.11 
show that the four states model matches well with the emulator outputs. 


4.5 Conclusion 


In this chapter, we introduced different types of fuel cells, i.e. PEMFC, AFC, PAFC, 
SOFC and MCEC. The advantages and disadvantages of these fuel cells were sum- 
marized. Then, we focus on the type of PEMFC, which is suitable for automobile 
applications due to its high energy density, low working temperature, simple struc- 
ture and long life. PEMFC is supplied with hydrogen and air at the anode and the 
cathode, respectively. 

A typical PEMFC system consists of four subsystems: the air feed subsystem, the 
hydrogen supply subsystem, the humidifier subsystem and the cooling subsystem. 
The detail models of the air feed subsystem including air compressor and fuel cell 
stack were given. The studied air feed system consists of a33 kW PEMFC containing 
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90 cells in series and a twin screw compressor. Under several practical assumptions, 
this model is reduced into a four states model while preserving the dynamic behavior. 
Finally, experimental validation was performed through the HIL test bench which 
consists of a physical air-feed system, based on a commercial twin screw compressor 
and a real-time PEMFC emulator. In the following chapter, disturbance-observer 
based SOSM control will be designed for the PEMFC system based on the proposed 
model. 
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Chapter 5 ®) 
Sliding Mode Control of PEMFC Systems sx" 


In this chapter, a model-based robust control is proposed for PEMFC system, based 
on SOSM algorithm. The control objective is to maximize the fuel cell net power and 
avoid the oxygen starvation by regulating the oxygen excess ratio to its desired value 
during fast load variations. The oxygen excess ratio is estimated via an ESO from 
the measurements of the compressor flow rate, the load current and supply manifold 
pressure. A HIL test bench which consists of acommercial twin screw air compressor 
and a real-time fuel cell emulation system, is used to validate the performance of 
the proposed ESO-based SOSM controller. The experimental results show that the 
proposed controller is robust and has a good transient performance in the presence 
of load variations and parametric uncertainties. 


5.1 Introduction 


Nowadays, fuel cells are widely regarded as a potential alternative energy conversion 
technology for stationary and mobile applications since they have high efficiency and 
low environmental pollution [5, 13, 23]. Among various types of fuel cells, polymer 
electrolyte membrane fuel cells (PEMFCs) are particularly suitable for automotive 
applications due to low temperature, quick start-up and favorable power-to-weight 
ratio [11, 34]. However, PEMFCs are multi-input multi-output systems with strong 
nonlinearities, and thus precise control efforts are required to guarantee the optimal 
performance. One of the most challenging tasks is related with the fuel cell air flow 
control, which means to regulate the oxygen quantity in the cathode. On one hand, 
during sudden load increase, it must ensure that the oxygen partial pressure supplied 
by the air compressor does not fall below a critical level (the value of oxygen excess 
ratio is less than one). Otherwise, the so-called oxygen starvation phenomenon occurs 
which will result in hot spots on the surface of the membrane, causing irreversible 
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damages to the polymeric membranes [8]. On the other hand, it must also reduce the 
parasitic losses of the air compressor since it can consume up to 30% of the fuel cell 
power under high load conditions [35]. Therefore, the management of air supplying 
to the fuel cell is one of the key factors for achieving desired fuel cell operation 
conditions. 

Many model based control strategies have been proposed for the control design 
of fuel cell air feed systems, ranging from linear state feedback control [26], optimal 
control [22], model predictive control [9, 32], intelligent control [6, 14] and SMC [8, 
23]. One of the main disadvantages of employing the intelligent control method [6] is 
the requirement of off-line training, which is difficult for real applications because of 
high computational burden. Furthermore, the neural network/fuzzy rule-based model 
does not provide explicit structure property of the physical system. As discussed in 
[14], the critical problem in fuel cell control design is the reconstruction of the oxygen 
excess ratio value. The accurate regulation of oxygen excess ratio can increase the 
system efficiency significantly [28]. Unfortunately, this value depends on the oxygen 
flow entering the cathode, which is unavailable for measurement. Therefore, state 
observers employed as a replacement for physical sensors to estimate the unavailable 
variables, are of great interest. 

Many observation/estimation strategies have been proposed for the fuel cell sys- 
tems in the last few years. Several kinds of Kalman filters (KFs) have been applied 
to the state estimation of the fuel cell systems [26, 33]. It should be noted that these 
methods are sensitive to external disturbances and modeling errors due to model lin- 
earization around pre-defined operating points of the system [2, 15]. Furthermore, 
the calculation of Jacobian matrix of the fuel cell system is time consuming. In [1], 
an adaptive observer was designed to estimate the states of PEMFCs. However, it 
lacks robustness against the fuel cell voltage’s measurement noise and is dependent 
on the system parameters. A finite-time convergent state observer based on HOSM 
was proposed for the PEMFCs, which allows for the finite-time estimation of the 
fuel cell states [23, 27]. Consequently, the separation principle is automatically sat- 
isfied, which means that the controller and the observer can be separately designed. 
This method relies on the calculation of the complex observability matrices which are 
highly nonlinear and state dependent, therefore difficult in real-time implementation. 
More recently, ESO has attracted great attention because it requires the least amoun- 
t of plant information. Unlike traditional observers, such as Luenberger observer, 
high-gain observer and unknown input observer, ESO regards the disturbances of 
the system as a new system state which is conceived to estimate not only the external 
disturbance but also plant dynamics. 

In this chapter, an ESO-based cascade controller is designed for regulating the 
oxygen excess ratio of the PEMFC air compression system to its desired value, using 
sliding mode technique. SMC is known for its finite time convergence, high accuracy 
and insensitiveness to uncertainties and disturbance [3, 12, 19]. These make it an 
effective method to deal with the nonlinear behavior of the considered system. The 
proposed controller is designed in a cascade structure because the relative degree 
between the oxygen excess ratio and the compressor quadrature current is two. It 
consists of two control loops, i.e., ESO-based oxygen excess ratio regulation loop 
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(outer loop) and compressor flow rate control loop (inner loop). The oxygen excess 
ratio is reconstructed by an ESO-based on the measurements of supply manifold 
pressure and compressor flow rate. The outer control loop, which uses the estimated 
oxygen excess ratio, provides the compressor flow rate reference for the inner loop 
based on the STA. The STA is an unique absolutely continuous sliding mode algo- 
rithm among the SOSM algorithms. Therefore, without any introduction of low pass 
filters [19], it is less affected by the chattering phenomenon. For the inner control 
loop, a simple SMC law consists of a linear term and a switching term is adopted, 
ensuring a fast response of the control scheme. The problems related to practical 
implementation of the proposed ESO-based cascade control on an HIL test bench 
has been resolved, and are the main focus of this paper. Experimental results are 
presented and the robustness of the proposed algorithm against parametric variation 
has been validated. 


5.2. Problem Formulation 


A schematic diagram of a PEMFC system is shown in Fig.5.1. The considered 
PEMFC system consists of four major subsystems: the air feed subsystem that feeds 
the cathode by oxygen, the hydrogen supply subsystem that feeds the anode side 
with hydrogen, the humidify subsystem and the cooling subsystem that maintain the 
humidity and the temperature of the fuel cells respectively. The air feed subsystem 
usually consists of an electromechanical air compressor that maintains the required 
oxygen pressure and mass-flow in the cathode of PEMFC. 

There are a number of models of fuel cells and air feed systems in the litera- 
ture, such as the nine-state model proposed by Pukrushpan et al. [26], one dimen- 
sional power density model proposed by Meidanshahi and Karimi [21], and multi- 
dimensional models developed in [4]. Although these works are useful for analytical 
purposes, they require large quantity of calculations. Thus, it is necessary to care- 
fully consider the model since some weakly coupled states can lead to a problem of 
observability of the states [33]. Suh [30] has reduced the nine-state model of Pukrush- 
pan et al. [26] into a four-state model, while preserving the dynamic behavior. The 
advantage of this model is that it has been validated in a wide operating range of a 
fuel cell by experimental results. 

By taking the dynamics of the vapor partial pressure into account, a 5-th order 
model of the PEMFC air feed system is considered in this work. This model facilitates 
both the control and observer design of the fuel cell air feed system in real-time. Some 
assumptions are used to derive a more accurate model of the PEMFC system which 
includes the cathode partial pressure dynamics, the air supply manifold dynamics 
and the compressor dynamics [27, 30]. Mainly, the temperature and the humidity 
of the air at the inlet of the fuel cell stack are assumed to be well controlled. The 
high pressure compressed hydrogen is available and is well controlled to follow 
the cathode pressure. The current dynamics of the motor which directly drive the 
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Fig. 5.1 Fuel cell system scheme 


compressor are negligible due to their small time constant as compared to the me- 
chanical dynamics. 

Some important assumptions have been considered for the system modelling [27, 
30]: 


e The temperature and the humidity of the air at the inlet of the fuel cell stack are 
assumed to be well controlled. This assumption is reasonable due to the fact the 
dynamics of the temperature and the humidity are very slow. 

e The high pressure compressed hydrogen is available and is well controlled to 
follow the cathode pressure. 

e The water inside the cathode is only in vapor phase and any extra water in liquid 
phase is removed from the channels. 

e The spatial variations are neglected, it is assumed that the gases and reactions are 
uniformly distributed in the cell. 

e The current dynamics of the motor which directly drives the compressor is negli- 
gible due to its small time constant as compared to the mechanical dynamics. 


5.2.1 The Vapor Partial Pressures Dynamics 


The dynamics of the vapor partial pressures are described by the following equation 


dp» ca RT rc 
= Wo,ca,in — Wo,ca,ou Wy,ca,gen Wy,mem) » 5.1 
7 Vy, (Wea ,ea,out + Wy,ca,gen + Wo,mem) (5.1) 


where Wy ¢a.in aNd Wy ca.our are the inlet and outlet flow rates of vapor, respectively, 
W,,ca,gen 18 the generated vapor flow and W, mem is the mass flow rate across the 
membrane. The detail calculation of these values can be found in [26]. 
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5.2.2 State Space Representation of PEMFC System 


In view of Eqs. (4.6, 4.9, 4.13, 5.1) and define the state vector, 


T T 
5 oe [x1, X2, X3, X4, x5] = [Weps Psm, POz; PNo> Pv.cal : (5.2) 


Denote hz := Wey and x := Pca = X3 + X4 + X5. Thus, the resulting nonlinear mod- 
el, suitable for control design, can be written in the following form 


. of (xm\" 
Xy = —-c4X1 1} ho + csu, (5.3) 


x] C3 


in = e614 (2) -1]| [is-es (x2 — x) ], (5.4) 


. C10x3 (x) 

X3 = C9 (x2 — x) os c14G, (5.5) 
C11X3 + C12X4 + €13X5 

: Ci0x4 (xX) 

X4 = C15 (X%2 — X) , (5.6) 
C11X3 + C12xX4 + C13X5 

; c10x5p (xX) 

X5 = C16 (x2 — x) ~~ + 2cyaC + €17, (5.7) 


C11X3 + C12xX4 + C13X5 


where u := i, is the control input, ¢ := Js is considered as a measurable disturbance, 
and the function ¢ (x) represents the flow rate at the cathode exit which depends on 
x. The parameters c;,i € {1,..., 21}, are given in Table5.1. 

As mentioned in [25], the vector of measurable outputs is given as follows: 


Y1 hy (x3, x4) 
y=|]yr}=] Ao) |, (5.8) 
¥3 X2 


where y) := hy = Vor, y2 = ho = Wep and y3 := X2 = Psm. It should be noted that 
all physical quantities are measured in the SI system of units. 


5.2.3 Control Objectives 


As the PEMFC system works as an autonomous power plant in automobiles, the com- 
pressor motor is also powered by the PEMFC. Therefore, the net power of the system 
can be expressed as (consumption by the other auxiliary systems is negligible): 


Pret = Pst — Pop, 
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Table 5.1 Parameters 
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where P;, and P,, are the powers produced by the fuel cell and the air-feed system, 
respectively. 

Experimental studies have shown that the air-feed system can consume up to 30% 
of the fuel-cell power under high load conditions [32]. Therefore, it needs to be 
operated at its optimal point, at which it supplies just sufficient oxygen necessary for 
the hydrogen and oxygen reaction. As already discussed in [26], accomplishing the 
power optimization is equivalent to maintain the oxygen excess ratio in an optimal 
value Xj, . 

The oxygen excess ratio Ao, is defined as [26, 28]: 


Wo,in C18 


AO, = = 
Wo, ,react c 196 


(Psm — X) (5.9) 


where Wo, in is the oxygen partial flow in the cathode and Wo, ;eacr 1s the reacted 
oxygen flow. Therefore, the main control objective for the fuel cell air-fed system 
is to force the oxygen excess ratio Ag, track \%,. Given that the relative degree 
between the oxygen excess ratio and control input is 2, an external control loop 
which regulates the compressor flow rate h2(x;) to its desired value At (cy) needs 
to be considered. The desired compressor flow rate is calculated in such a way that 
AO> > 5 


Remark 5.1 High oxygen excess ratio, and thus high oxygen partial pressure, im- 
proves the net power P,,-; and the stack power P;,. However, above an optimum XO, 
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level, further increase will cause excessive increase of compressor power and thus 
deteriorate the system net power. For different load current /,, drawn from the fuel 
cell system, there is an optimum \%, lies between 2 and 2.5 where the maximized 
net power is achieved. 


5.3 Control Design 


5.3.1 Oxygen Excess Ratio Control 


In view of (5.9), the main obstacle of calculating the oxygen excess ratio Ao, is 
the unavailability of the cathode pressure y which is difficult to be measured. Fur- 
thermore, it is not always possible to use sensors for measurements, either due to 
prohibitive costs of the sensing technology or because the quantity is not directly 
measurable, especially in the conditions of humidified gas streams inside the fuel 
cell stack. In these cases, state observers serve as a replacement for physical sensors, 
for obtaining the unavailable quantities, are of great interest [7]. 

Hence, we propose an ESO which provides satisfactory level of accuracy and fast 
convergence for reconstructing the stoichiometry \o,. Then, the estimated value is 
used in the output feedback controller design which takes a cascade control structure 
(see Fig. 5.2). The proposed control structure consists of two control loops: outer loop 
and inner loop. The outer loop generates a reference compressor flow rate based on 
the oxygen excess ratio error. The inner loop computes the control input u based 
on the error information between the reference compressor flow rate and measured 
value. The super-twisting (ST) sliding mode algorithm is employed in both outer and 
inner control loops. 


Extended State Observer Design 


The Eq. (5.4) can be rewritten as follows: 


X2 = fi(%2) (hz — cgx2), (5.10) 
hz pref 
ad “ Eatemtal ae a 82 Inner bal , FC air-feed 
= controller (5.19) ust i, controller (5.24) system 
2 
t a 


aE .—__ 


(5.17)¢ ESO (5.12) DED, 


Fig. 5.2 Block diagram of the oxygen excess ratio control system 
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where fi(x2) = c6 {1 ie [(2) = 1] and d(t) = cg fi (x2) x. 


C3 


Remark 5.2 Given that x2 is a physical variable and the coefficients c3, c4, co and 
C4 
c7 are all positive, the term (2) — 1 > Ois assured in the operation range, which 


means f}(x2) is also positive. Thus, it is reasonable to consider that f|(x2) is a 
bounded invertible function in the operation range x2 € X2. 


Denote z; := x2 and z2 := d(t). Assume that f; (x2) and x are both differentiable, 
it follows z2 = h(t) with h(t) the variation rate of system uncertainty and disturbance 
d(t). Then, (5.10) can be rearranged as, 


Zz = S2(Z1, ha) + Z2, 
Z2=hi(t), (5.11) 


in which f2(z1, 42) = fi (Z1) (hz — €gz1). A linear ESO is designed with the follow- 
ing structure [10, 31]: 

21 = frlzi, fo) + 22 + i (z1 — 21), 

by = Br (zi — 21), (5.12) 


where Z = [Z1, Zo] € R’, the designing gains 3, and (3) are chosen such that the 
polynomial \? + 3, + is Hurwitz stable. Therefore, its natural frequency w, and 


damping ratio € are: 
a ed 
Wn = Bo, f= : 
24/ D2 


Denote the observation errors Z} = z, — Z; and Z2 = z2 — 2, the error dynamics are 
given by, 


Z = —$,Z; +2, (5.13) 
Zo = —Poz1 + h(t). (5.14) 


The system (5.13) can be written in state space form as follows: 


Zi2 = AZ + Bh(t), (5.15) 


where Z)2 = BE — ee al and B = BE 
2 —P2 


Lemma 5.1 Assuming that h(t) is bounded, there exist a constant 6 > 0 and a finite 
time T, > 0 such that the trajectories of the system (5.13) are bounded, ||Z12|| < 6, 
Vt > T, > Oand 6 is dependent on the initial condition of 2,2 and the upper boundary 


of hit). 
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Proof Under the assumption of the boundedness of h(t), the bounded-input- 
bounded-output stability for the system (5.13) is ensured. The detailed proof can 
be obtained directly from the works [31, 36]. | 
From Remark 5.2, f) (x2) is a bounded invertible function in the operation range. 
Then, it follows from Lemma 5.1 that the estimate of y can be constructed as, 


X= cg fi (x2)22, for fi(x2) £0. (5.16) 


Remark 5.3 Since h(t) is unknown in the error dynamics (5.14), it should be noted 
that the estimation of y in (5.16) is not asymptotic stable in the case when h(f) is not 
negligible. For the system having fast varying disturbances/uncertainties, i.e., h(t) 
is large, the performance of ESO presented above does not offer satisfactory results. 
To deal with this issue, the concept of higher order ESO can be employed. The idea 
is simple, if h(t) is not negligible but its derivative h(t) is negligible, then a second 
order ESO can be used. Therefore, in application, one can design a r-th order ESO 
if h—) (t) is negligible [31]. 


Remark 5.4 The observer gains (3; and (3, need to be determined, which decides the 
bandwidth of the ESO. For instance, we can choose 3; = 23 > 0 and 3 = (3, in 
view of (5.15), the larger value ( is, the more accurate estimation is achieved. How- 
ever, this increases the noise sensitivity due to the augment of the bandwidth. More 
generally speaking, the function h(t) represents the rate of change in the physical 
world. If this value is quite large, the physical variables will change very rapidly. In 
this case, the observer bandwidth needs to be sufficiently large for an accurate esti- 
mate of z2. Therefore, the selection of 3 should balance the estimation performance 
and the noise tolerance. 


Outer Loop 


From the result of Lemma 5.1, the oxygen excess ratio can be reconstructed as, 

. C18 - 

No, = —> (x2 — X). (5.17) 
The sliding variable is defined as follows: 

iQ) = dor = Nae 

where AG, is the reference of oxygen excess ratio, whose value is between 2.0 and 
2.5. A variable setting for the reference of Ag, is given in order to not only obtain 
the maximum efficiency from the fuel cell but also avoid the dangers of oxygen 


starvation in the presence of fast load variations [8, 16]. 
Take the first time derivative of the sliding variable sj (t): 
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5n(t) = c2o| fila) (ha — cor + 8x) — 5] 
= gi (x2) ha + yilx2, x, x), (5.18) 


where g1(%2) = gio + Agio, and gio > 0 is a known function based on nom- 
inal parameters while Ago represents parametric uncertainties, yi(x2, xX; Y= 
Pi (x2, X xX; x) + Ay}. It is assumed that 0 < aol <b; <1 and Ay, is assumed 
to be bounded and with first derivative bounded [29]. Both g; and y; are smooth 
enough in the operation range. 

Considering the convergence and stability of the system, the STA can be employed 
for compressor mass flow control in order to force s;(t) and s;(f) to zero in finite 
time. The desired compressor flow rate h} is calculated as follows: 


1 A an 
A= [-o (1) - Giles, %%)], (5.19) 


where a(s,) = Aj |s1|2 sign(s,) +a, fosign(s\)dr is the STA, A; and a, are design 
parameters. Substitute (5.19) into (5.18), the sliding manifold can be rewritten as 


follows: A 
51 (0) --( 20) (si) + Api - 
0 


=-—y10(s1) + 1. 


Agio » 
1 
0 


Through the application of the extended state observer, it can yield mat the term 
go, = Ay, - aoe (1 is bounded and with first derivative bounded, i.e., + 4 \bi| < ®, 
for some positive value ®,. The sufficient conditions for the finite time convergence 
to the sliding manifold s;(t) = s)(t) = 0 are [18]: 


P| Da 49, a, + @, 
t=) 1b) a =O, 


(5.20) 


The boundary values of these functions need to be determined a-priori through exper- 
imentation. Then, the controller gains a; and A; can be selected according to (5.20). 
In this study, we have found the following boundary values as ®; = 4.2 x 10°, 
b, = 0.4. Accordingly, the gains are tuned as a; = 7.5 x 10° and \; = 425. 

The controller (5.19) is followed by a linear first order filter which is given as 


follows [24]: 


Ent = hin, (5.21) 
dt 2 2, 2 


The filter output ae is the command compressor flow rate that provides the inner 
control loop. 


Remark 5.5 It can be seen that (5.19) depends upon x. A real-time robust exact 
differentiator proposed in [18] can be used to estimate the time derivative of ¥ in 
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finite time. The differentiator has the following form 


q= me area LP sign(G, —-XY+, 
Go = —a2Lsign(¢ — x), (5.22) 


where ¢; and ¢) are the real-time estimations of ¥ and x, respectively. The parameters 
of the differentiator 42 = 1, a2 = 1.1 are suggested and L is the only parameter 
needed to be tuned according to the condition L > |x}. 


Inner Loop 


The inner loop will generate the compressor motor quadrature current u based on 
SMC. The sliding variable s2(t) is defined as follows: 


so(t) = hy — hE", 
Take the first time derivative of s(t): 


2 C2 X2 * ref 
$2 (t) = C21 yx] lL} hg+csut —hy 
xX] C3 


= gou t+ Y2 (a1, x2) — AR, (5.23) 


where g2 = 9290 + Agro and Y2 = Yoo + AYr are smooth enough in the operation 
range, (29 > O and gzo are known functions based on nominal parameters while 
Av2 and Agog are parametric uncertainties. It is assumed that 


A 
( 2icy, 
920 


and Ao is assumed to be bounded and with first derivative bounded [29]. 

Note that the sliding variable s2(t) has relative degree one with respect to the 
control input u, because u appears in the s(t) expression explicitly. Thus, the internal 
controller can be designed as follows: 


I 
u = — (—k,sign (s2) — k282 — p20) , (5.24) 
920 
for some positive constants k; and ky. Applying (5.24) to (5.23), the sliding manifold 
can be rewritten as follows: 


$2 = — 2k, sign (s2) = y2k252 =P 2, 
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where 72 = (1 + As ) and ¢2 = Ayo — on yx — hie. Assume that the term ¢ 


is bounded, i.e., |¢2| < ®2, for some positive value 2. The sufficient conditions for 
the finite time convergence to the sliding manifold s7(t) = s(t) = 0 are: 


Py 


es. 
er a 


ki > 


In this study, the following boundary values are obtained as 2 = 150, bz = 0.4. 
Accordingly, the gains are tuned as k; = 180 and ky = 25. 


Remark 5.6 It is clear to see that the control law (5.24) consists of a linear term 
kosy and a switching term k;sign(s2). When the states trajectories are far from the 
surface, the linear term improves the behavior of the controller with a high rate of 
convergence. Conversely, the switching term dominates the the rate of convergence 
when the states trajectories are in the vicinity of the surface. Therefore, kz > 0 can 
be chosen according to the desired rate of convergence towards the sliding surface 
from the initial condition. 


5.3.2 Closed-Loop Stability Analysis 


Theorem 5.1 Consider the system (5.3)-(5.7) in closed loop with the controller 
(5.24), satisfying the conditions in Lemma 5.1 and a smooth reference profile i in 
(5.21). Then, there exists finite time T > 0 such that 


Ao, — a4, | = 10,0), #>T, 


where 6, is a positive function depending on 6 and |. 
Proof The detailed proof can be followed from [24]. | 


Remark 5.7 The reason for introducing the smoothing filter (5.21) is that, the sta- 
bility analysis of the internal compressor flow rate loop requires a bounded first time 
derivative of i [24]. As can be seen from (5.19) that, h> has a discontinuous time 
derivative because of the term o(s,), thus, a first order filter can be used to perform 
the required smoothing action. 


5.4 Experimental Results 


The ESO-based control proposed in Sect. 5.3 is implemented on an HIL test bench 
shown in Fig. 4.5, where the implementable feasibility of the algorithm is evaluated. It 
consists of a physical air-feed system, based on a commercial twin screw compressor 
and a real-time PEMFC emulator derived from the work of [17, 20]. 
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Fig. 5.3. Block diagram of the HIL test bench 


The nominal values of the parameters for the HIL emulator are shown in Table 4.4. 
The PEMFC works under several safety constraints, the most important ones are the 
anode and cathode pressure difference (to be kept minimum) and stack temperature 
(to be regulated). These constraints r equire additional control in real fuel cells, which 
have been replicated in the emulator as well. An internal controller is implemented 
in the emulator to ensure that the anode pressure follows the cathode pressure at all 
times. The temperature is also controlled and its value can be set through software. 
The block diagram of this HIL test bench is shown in Fig. 5.3. 

In order to test the proposed control’s robustness against parametric uncertainties, 
the perturbations of the system parameters are given in Table5.2, i.e., volumes, 
cathode inlet and outlet orifice, ambient temperature and stack temperature, motor 
constant and compressor inertia have been varied around their nominal values (see 
Table 4.4). During the tests, the load (stack current) is varied between 100 and 300 
A as shown in Fig. 5.4. The corresponding compressor flow rate is shown in Fig. 5.5. 

The results of the ESO-STA controller are compared with that of STA, as shown 
in Figs.5.6, 5.7, 5.8, 5.9 and 5.10. The performance of the oxygen excess ratio by 
the proposed controller and STA is shown in Fig. 5.6. It is easy to see that the ESO 
is capable of reconstructing oxygen excess ratio in the presence of parameter uncer- 
tainties, which guarantees that the ESO-STA controller recovers the performance of 
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Table 5.2 Variations of system parameters 


Symbol Parameter Value (%) 
Tst Temperature of the fuel cell, °C +10 

Vea Cathode volume, m? +10 

Vom Supply manifold volume, m3 —10 

Tatm Ambient temperature, °C +10 

Kea,in Cathode inlet orifice constant, kg/(Pa-s) +5 

Kea,out Cathode outlet orifice constant, kg/(Pa-s) +5 

kt Motor constant, N-m/A —10 

Joep Compressor inertia, kg-m? +5 


Stack current (A) 


i) 10 20 30 40 50 60 70 80 90 100 110 


Fig. 5.4 Stack current under load variation 


the state feedback STA controller when the gain / is chosen appropriately (3 = 1000 
is chosen). Little overshoot and fast convergence can also be observed from Fig. 5.6, 
i.e., att = 60s, sudden increase in load current (from 200 to 300 A) causes a sudden 
drop in the oxygen excess ratio (from 2.2 to 1.51). Conversely, at t = 80 s, sudden 
decrease in load current (from 300 to 200 A) causes a sudden increase in the oxygen 
excess ratio (from 2.2 to 2.91). It can be seen from Fig. 5.6a, b that the proposed con- 
troller has a good response and the error remains within an acceptable neighborhood 
of zero. In addition, the settling time of both upward load changes and downward 
load changes is less than 1.5 s for both controllers. 

The performance of the Levant’s differentiator [18] is shown in Fig. 5.7, where fast 
convergence and little overshoot can be observed. The compressor flow rate tracking 
its reference is shown in Fig.5.8. The compressor flow rate reference is generated 
from the external regulation loop, which is based on the error between the desired 
value of oxygen excess ratio and its estimated value. 
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Fig. 5.5 | Measurement of compressor flow rate under load variation 
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Fig. 5.6 Performance of oxygen excess ratio (a) and its regulation error (b) 
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Fig. 5.7 Performance of the differentiator (5.22) 
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Fig. 5.8 Performance of the compressor flow rate 
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Fig. 5.10 Fuel cell net power 
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It should be pointed out that the output of the external regulation loop is compared 
with the experimental measurement. The similar behavior of the two controllers can 
be found, which again guarantees that the proposed ESO-STA controller can recover 
the performance of the state feedback STA controller. The associated stack voltages 
and net power of both controllers are shown in Figs.5.9 and 5.10, respectively. The 
stack voltage varies between 62.5 and 66.5 V, while the net power of the fuel cell 
varies between 6.2 and 18 kW, according to the load current. 


5.5 Conclusion 


In this chapter, an ESO-based cascade controller for regulating the oxygen excess 
ratio of the PEMFC air feed system to its desired value, using sliding mode technique, 
has been presented. The control objective is to avoid oxygen starvation during sudden 
load changes, which is considered as a measurable disturbance in this work. The 
designed cascade controller consists of oxygen excess ratio tracking (outer) loop and 
compressor flow rate regulation (inner) loop. 

Firstly, an ESO is employed to reconstruct the oxygen excess ratio with high accu- 
racy and fast convergence. The outer control loop, which uses the estimated oxygen 
excess ratio, provides the compressor flow rate reference for the inner loop based 
on the STA. Secondly, a simple SMC law consists of a linear term and a switching 
term is designed for the inner loop, which ensures a fast response of the control 
scheme. Finally, the proposed ESO-based cascade control approach was successful- 
ly implemented on an HIL test bench. The experimental results have demonstrated 
the effectiveness, robustness and feasibility of the proposed controller. Furthermore, 
the obtained results are compared to state feedback STA control under the assump- 
tion that oxygen excess ratio is available, which show that the performance of the 
proposed ESO-based STA control can recover that of the state feedback control in 
the presence of parameter uncertainties. In the following chapter, nonlinear observer 
design for PEMFC systems using sliding mode technique will be investigated. 
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Chapter 6 M®) 
Sliding Mode Observer of PEMFC creek 
Systems 


As mentioned in the general introduction, it is not always possible to use sensors for 
measurements, either due to prohibitive costs of the sensing technology or because the 
quantity is not directly measurable. For precise control applications, state observer 
can be used for obtaining unavailable state values instead of sensors. A brief survey 
of existing methods in order to define the context of our work. 

The problems addressed in this chapter are the design of SOSM observer and alge- 
braic observer for PEMFC systems. For the SOSM observer, the goal is to estimate 
the hydrogen partial pressure in the anode channel of the PEMFC, using the mea- 
surements of stack voltage, stack current, anode pressure and anode inlet pressure. 
The proposed observer employs a nonlinear error injection term, where the error is 
obtained from the difference between the system voltage output obtained from an 
experimental validated nonlinear model and estimated voltage output obtained from 
the designed observer. The robustness of this observer against parametric uncertain- 
ties and load variations is studied, and the finite time convergence property is proved 
via Lyapunov analysis. For the algebraic observer, the goal is to estimate oxygen 
and nitrogen partial pressures in the fuel cell cathode side, using measurements of 
supply manifold pressure and compressor mass flow rate. As the proposed technique 
requires the time derivatives of the states, Lyapunov based adaptive HOSM differ- 
entiators are synthesized and implemented for estimating these derivatives without 
a-priori knowledge of the upper bounds of their higher order time derivatives. 


6.1 Introduction 


PEMEFC is an electrochemical device that produces electricity from the chemical 
reaction between hydrogen and oxygen [16, 27]. The electrolyte membrane of a 
PEMFC has a special property that allows only positive protons to pass through 
while blocking electrons. They are under intensive development in the past few 
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years as they are regarded as potential alternative power sources in automotive appli- 
cations due to their relatively small size, low temperature (40-—180°C), quick start up 
and easy manufacturing [41]. Recent developments in PEM and catalyst technology 
have greatly increased the power density of fuel cells, made them one of the most 
prominent technologies for future’s automotive world [25]. 

Automotive fuel cell applications have more rigorous operating requirements 
than stationary applications [7], therefore these applications need precise control 
of performance, in order to guarantee the reliability, health, and safety of both, the 
fuel cell and the user. Along with control, health monitoring and safety systems 
are essential for the application of fuel cells in automotive systems. Both control 
and health-monitoring systems require precise measurements of different physical 
quantities in the fuel cell. However, it is not always possible to use sensors for 
measurements, either due to prohibitive costs of the sensing technology or because 
the quantity is not directly measurable. In both these cases, state observers serve as 
a replacement for physical sensors, for obtaining the unavailable quantities. 

In PEMEFC, hydrogen which is generated from the fuel processing system is fed 
into the anode side of the cell stack, while air is pumped into the cathode side through 
an air compressor. One hand, it is well known that the fuel cell stack life is reduced 
due to the so-called starvation phenomenon of the cell, because of insufficient supply 
of oxygen and hydrogen [49]. On the other hand, excessive supply is also undesirable 
due to the reduction of its efficiency. Pressure regulation remains one of the most 
challenging control problems. The main objective is to ensure minimum pressure 
difference between the anode and cathode side of the PEMFC, avoid the membrane 
from damage and increase the fuel cell stack life [5]. Several model based control 
approaches have been proposed for solving this problem, such as H,, robust control 
[31], feedback linearization plus pole placement [39], proportional integral (PI) plus 
static feed-forward controller [41] and static feedback controller [23]. To evaluate 
the performance of these controllers, a major obstacle is the absence of reliable 
measurements of hydrogen partial pressure, specially in the conditions of humidified 
gas streams inside the cell stack. However, it is not always possible to use sensors for 
measurements, either due to prohibitive costs of the sensing technology or because 
the quantity is not directly measurable. The sensors that do provide satisfactory 
performance usually suffer from slow response times, low accuracy, bulky and high 
cost [2]. Therefore, state observers serve as a replacement for physical sensors, are 
of great interest. 

It is known that low oxygen partial pressure in the cathode reduces the generation 
capacity of fuel cell systems and affects the fuel cell stack life [35]. Therefore, the 
observation of oxygen partial pressure is essential for the feedback control and fault 
detection, in order to ensure the safety and longevity of the fuel cells [21, 44, 50]. 
However, most of the commercially available oxygen sensors do not operate properly 
in presence of humidified gas streams inside the fuel cell stack [2]. The sensors that do 
provide satisfactory performance are usually too big and costly. Therefore, observer 
designs for estimating the unmeasurable partial pressures are of great interest. 

Many research endeavors have been focused during the recent years, on obser- 
vation problems in fuel cell systems [8, 15]. Arcak et al. [2] developed an adaptive 
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observer for hydrogen partial pressure estimation based on the fuel cell voltage. 
Gorgiin et al. [14] developed an algorithm for estimating partial pressures and the 
membrane water content in PEMFCs based on the resistive cell voltage drop. This 
algorithm has incorporated two adaptive observers for hydrogen and oxygen partial 
pressures, adapted from the work of Arcak et al. [2]. However, both of the above 
works lack robustness against the fuel cell voltage’s measurement noise and the 
internal model relies upon unmeasurable values. Several kinds of KFs have been 
applied to the state estimation of fuel cell systems, i.e. classical KF [42], Unscented 
Kalman Filter (UKF) [38] and adaptive UKF [53]. These approaches are based on 
model linearization around pre-defined operating points of the system. Moreover, 
the calculation of Jacobian matrix of complex models like fuel cells are time con- 
suming, and therefore difficult in real-time implementation [4, 22]. Ingimundarson 
et al. [19] proposed a model based estimation approach to hydrogen leak detection 
in PEMFCs. More recently, Linear Parameter Varying (LPV) observer was proposed 
by Lira et al. [32] for the application to fault detection in PEMFC systems, where 
the stack current was taken as the scheduling variable. Kunusch et al. [26] proposed 
an approach based on super-twisting algorithm, in order to estimate the hydrogen 
input flow at the anode of the stack and the water transport across the membrane. 
A Luenberger observer was employed by Thawornkuno and Panjapornpon [51] in 
order to estimate the membrane water content in PEMFCs. The main limitation of 
this method is that it can only converge to a neighborhood of the real system states 
in the presence of disturbances. 

Due to the lack of a straightforward observer design method for a given nonlinear 
system, many observation methods are generally dependent upon state transforma- 
tions, the structure of the system, the form of the nonlinearities and the boundedness 
of the system states [18]. Among the popular strategies, high gain observers [13] 
are usually employed to estimate the system states under the assumption that the 
nonlinearity vector is globally or locally Lipschitz. However, in practice, the Lip- 
schitz constraint is not easy to obtain, which prevents the global convergence of 
the high gain observer. Although the circle-criterion observer design [3, 11] relaxes 
the requirement of Lipschitz constant, it remains limited to systems with positive- 
gradient nonlinearities. 

SMOs have found wide application in the areas of parameter estimation [1], state 
estimation [47] and FDI [54, 55] in recent years. Their well-known advantages are 
robustness and insensitivity to external disturbance. HOSM Observers have better 
performance as compared to classical sliding mode based observers because their 
output is continuous and does not require low pass filters. The STA, which is an 
unique absolutely continuous sliding mode algorithm among the SOSM algorithms, 
therefore it does not suffer from the problem of chattering [28]. The main advantages 
of the STA are: it does not need the evaluation of the time derivative of the sliding 
variable; its continuous nature suppresses arbitrary disturbances with bounded time 
derivatives. 
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6.2 SMO Design for PEMFC Systems 


6.2.1 Nonlinear Dynamic Model of PEMFC 


Figure 6.1 shows a block diagram of a typical PEM fuel cell system [33]. It consists of 
four major subsystems: the hydrogen supply subsystem, the air feed subsystem, the 
humidifier subsystem and the cooling subsystem. As our study is related to the hydro- 
gen supply subsystem, we will restrict the model for observer design specifically to 
the states of the hydrogen supply system and PEMFC anode. Certain assumptions 
have been imposed on the operating conditions, as follows [2, 43]: 


Al. The stack temperature and humidity in the fuel cell cathode and anode are well 
controlled; 

A2. The anode pressure is well controlled to follow the cathode pressure; 

A3. The temperatures inside the anode and the cathode are equal to the stack tem- 
perature; 

A4. The flow channel and cathode backing layer are lumped into one volume which 
assumes uniform conditions inside the anode channel; 

AS. Vapor partial pressure in the stack is considered equal to the saturation pressure. 


Remark 6.1 The first two assumptions are valid from practical point of view as 
the temperature, humidity and anode pressure are usually regulated externally in 
PEMEFC applications. The third assumption is justified as the temperature dynamics 
of a PEMFC stack are slow [43]. Assumption 4 is employed because we are interested 
in estimating the hydrogen partial pressure in the exit, and not its distribution along 
the channel [2]. Assumption 5 means that if the gas humidity drops below 100%, 
liquid water will either evaporate into the cathode gas or it will accumulate in the 
cathode. 
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Fig. 6.1 Fuel cell system scheme 
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6.2.2 Gas Dynamics of Hydrogen Partial Pressure 


Under these considerations, a lumped model is given as follows, 


RT fc 


Hy = ——— (W in — W, out — W. react) >» 6.1 
PH, V,.Ma, | Fb, Hp, out Fb, +) (6.1) 


where M7, is the molar masses of hydrogen Wz, jn is the mass flow rate of hydrogen 
gas entering the anode, Wy, ou; is the mass flow rate of hydrogen gas leaving the 
anode, Wy, react 18 the rate of hydrogen reacted, R is the universal gas constant, and 
Van 1s the anode volume. 

The inlet and outlet mass flow rates of hydrogen are given as follows, 


1 

Winin = Tron Wanin, (6.2) 
1 

Whp,out = Thon Mansour» (6.3) 
an,out 


where Wgn,in ANd Wan,in are the anode inlet, outlet humidity ratios, respectively. Wan, in 
and Wan,our are anode inlet and outlet mass flow rates, respectively. Given that the 
purge of the anode gas is zero, they can be expressed as, 


M, Pan,in Psat (T fc) 
Mu, Pan — Gan,inPsat(T fc)’ 
Warsn = Kania Dania = Pan) 
Wan,out = 0, 


Wanjin = 


(6.4) 


where M, is the molar masses of vapor, @an,in 18 the relative humidity on the anode, 
Psat (Tc) is the saturation pressure at fuel cell stack temperature. Kgn,in is the nozzle 
inlet flow constant for the anode side, Pan,in is the anode inlet pressure, and pay is 
the anode pressure. 

The reacted mass flow rate is defined by, 


ni, 
Wp react = Mine (6.5) 


where, F is Faraday constant, and J,, is the stack current, which is considered as a 
measurable disturbance variable. 
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6.2.3 SOSM Observer Design 


In this Section, we are interested in the estimation of hydrogen partial pressure py,, 
thus we will restrict the model for observer design specifically to the states of the 
PEMFC anode. We assume that the oxygen partial pressure po, is available. This 
assumption is reasonable because we can simply employ a separate observer for the 
cathode channel, e.g. an adaptive numerical differentiation approach is proposed in 
Liu et al. [33] for estimating this value. Unlike partial pressures, the total pressures 
(Pan.ins Pan) are available for measurement. Therefore, an SOSM observer is devel- 
oped for estimating the hydrogen partial pressure from the measurement of output 
stack voltage. 
In view of Eqs. (4.1) and (4.2), it follows that 


RT; 1 ; 
Vfo=n {e + — c (PH) + ne = vi| F (6.6) 


where v(J;,) represents the voltage drop as a function of current. The dynamics of 
Pu, in (6.1) can be rewritten as, 


Puy =Cc (Danis = Das) = c2¢, (6.7) 


Ruy T fe Kan.in = Ruy Ty nMny, 


v. Ta. =  e and ¢ := I,,. Denote py, as the estimate 
of hydrogen partial pressure. Since the voltage V+, is a function of py, and po,, by 
substituting in py, in (6.6), the voltage can be predicted as follows 


where cy = 


4 RT fc ; In(po, 
Vro=n {« zy a c (Bu,) + mee = v0 (6.8) 


The output error term Vie = Vee — Vie can be computed as 


ns RT yf F 
a sp (in (pa) — In (Brn) (6.9) 


Noting from the Mean Value Theorem that, 


In (pu) — In (Py) = & (Pm — Pu)» (6.10) 


for some « > 1. Substituting (6.10) into (6.8), we have, 


Vie = ¢3 (Pi, — Pan) (6.11) 
where c3 = ae The observer for estimating the hydrogen partial pressure py, 
can be designed by adding the correction error term from the measurement of stack 
voltage, 
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Pm =C) (Pang — Pan) + UV fe), (6.12) 


where ju(s) is the output error-injection term generated from STA 


Ls) = Ms|2sign(s) + af sign(s)drT, (6.13) 
0 


X and a are some positive constants. 
The observation error is defined as e := py, — pu,. Subtracting (6.12) from(6.7), 
the error dynamical equation is described by 


é = —p(c3e) — C2¢. (6.14) 


Theorem 6.1 Consider the system (6.1), where the voltage Vr,, current I,,, tem- 
perature Tr,, partial oxygen pressure Po,, total anode pressure Pay and total anode 
inlet pressure Pan,in are available. Then, the SOSM observer (6.12) and (6.13) guar- 
antees the trajectories of the error system (6.14) to zero in finite time, given that 
and a are high enough. 


Proof The detailed proof can be followed from Moreno and Osorio [36] directly. 


Robustness Analysis 


In this part, we will show the robustness of the proposed SOSM observer against 
model uncertainties and measurement disturbances. First, bounded disturbances will 
be introduced in the right-hand sides of (1) and (8), respectively. Thus, we have 


PH, = C1 ( Dancin _ Pan) = co + p(t), 


6.15 
Vee Ve eo + 6(t), 


where Vy, denotes the nominal value in (6), 5(t) incorporates model inaccuracies 
and other disturbances acting on the plant, and e(t) represents the measurement noise 
of stack voltage. 

In view of (6.15) and (6.12), the error dynamical equation is described by 


é = —p(c3e + e€(t)) + p(t), (6.16) 


where p(t) = d(t) — c2¢, the gains and a in the STA algorithm ju(-) are formulated 
as 


A=AVL, a = aL, (6.17) 


where Ao and av are positive constants and L is the only parameter needs to be tuned. 


110 6 Sliding Mode Observer of PEMFC Systems 


Theorem 6.2 Assume that the conditions in Theorem 6.1 hold and the disturbances 
p(t) and measurement noise €(t) satisfy the following conditions 


IPO! Sx, lel so. (6.18) 


Then, the following inequality is established in finite time, 


Ile < — S. (6.19) 
Amin(Q1) — Allgalle — 2\lasllaz 
if 
2xllasil2 
L> = ‘ (6.20) 
Xmin(Q1) _ & — Allgalle 
where 
1 5 fA? +2a —A 
2(0) = [IgillaQ2o)?, Q1 = a] a ae 
x 7 x ve ne 6.21 
gy = (02 +4a),-V], ge = [-A0, 28, en 
gs = [-A,2], A = 3A, = & = C300. 
Proof Denote e; := c3e. The system (6.16) can be rewritten as 
ai = Alex + el? signte +6) +e, eons 
é2 = —asign(e; + €) + p(t), 


where e; = c3e, A= c3A, @ = c3q@ and p(t) = ¢3(t). 
A new State vector is introduced to represent the system (22) in a more convenient 
form for Lyapunov analysis. 


1 eee 
t= le _ pe? (6.23) 
©) 2) 
Denote ; 
A; = lei? sign(e1) — ler + €/? sign(e: + €), (6.24) 
A». = sign(e;) — sign(e; + €). 
In view of (6.18), it follows that 
Ay] < V202, |Ay| < 2. (6.25) 


Then, the system (6.22) can be rewritten as follow 
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: Lf = I A\ 
= —~- = 2 a 6.26 
‘ 2\61| & i - Ba sil oe 


The following Lyapunov function candidate is introduced for the system (6.26) 
et ae ee ee (6.21) 
, 2 -vA 2]? . 


where the matrix P is symmetric positive definite due to the fact that its leading 
principle minors are all positive. Taking the derivative of (6.27) yields 


2A 
V=- ae are ae ais + Lada + est, (6.28) 
where =" _ 
> [7 + 2a -) : ee 
2) 2 , = 2 4 — 2 ; 
QO; = | qi = [Q? +4a), —A’] 650) 
qo. = [—Aa, 24], qx = [-A, 2]. 


It is easy to verify that Q, is a positive definite matrix. Since \min(P)||¢ 13 <V< 
Xmax(P)|I¢ le where || - ||2 denotes the Euclidean norm of a vector or the spectral 
norm of a matrix. Then, Eq. (6.28) can be rewritten as 


: i 
a oe 2 6.30 
< 7G | (Qi) |IClls + sa | 2(o)|ICll2 (6.30) 


+ 2L ||qallallCll2 + xilasllailClle. 


1 
where 21 (0) = IIqill220)?. 
Suppose that there exists a positive constant & such that 


2\(0) < AillCll2- (6.31) 


Therefore, it follows 


V Aes eee 
<- aK (Q1) —&) WClla + 2£L [lgellallClla + xllasilallclle (6.32) 


= =y Clas 


where y = £ (min(Q1) — &) — 2L||qall2 — x|l¢3ll2. Under the condition (6.20), we 
have 


(2 e 


6.33 
~ V Xmax(P) : 
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By the comparison principle [24], the differential inequality (6.32) is finite time 
convergent if 


LQ\(o) 
IWCllo > : (6.34) 
LAmin(Q1) — 4L |lg2ll2 — 2x\laail2 
Therefore, we get 
Clo < te) (6.35) 
2 = . . 
LAmin(Q1) — 4L |lgall2 — 2x\lasll2 
This completes the proof. | 


6.3 Algebraical Observer Design of PEMFC Systems 


The main content in this chapter is the design of an algebraic observer for the PEMFC 
system in automotive applications. This observer is designed for observing the partial 
pressures of oxygen and nitrogen in the cathode of the PEMFC. The necessity of 
this type of observation arises typically in fuel-cell powered vehicles where air is 
used as oxygen source in the cathode of the PEMFC. As air is composed primarily of 
nitrogen and oxygen, the partial pressure of each gas needs to be estimated in order to 
ensure the existence of sufficient quantity of oxygen in the cathode for the chemical 
reaction. Insufficient oxygen quantity can result in low power output and also lead 
to permanent physical damage to the PEMFC [2]. As partial pressures cannot be 
measured from conventional sensors, their observation is essential for the control 
and fault detection, in order to ensure the safety and longevity of the PEMFC [21, 
45, 50]. 

The motivation behind this work is that algebraic observers [10, 17] are precise 
and easily implementable in automotive embedded systems. We first demonstrate 
that the states of the PEMFC air-feed system is algebraically observable, i.e. they 
can be presented in terms of a static diffeomorphism [20] involving the system out- 
puts (compressor flow rate and supply manifold pressure) and their time derivatives, 
respectively. As the algebraical observer requires the time derivatives of the output 
variables, its performance requires robust and exact differentiators [10]. The existing 
robust differentiators, such as Levant [29], require a priori knowledge of the upper 
bound of a higher order time derivative, i.e. the Lipschitz constant. In many practical 
cases, this boundary can not be easily obtained and the results are possibly very con- 
servative. To overcome this problem, we propose new Lyapunov-based adaptive first 
and SOSM differentiators for practical implementation of our algebraic observation 
scheme. The adaptive gains adjust dynamically, therefore the Lipschitz constant is 
not required during the design. 

The performance of the proposed observer is evaluated by implementing on the 
instrumented HIL test bench described in Chap. 4. In our experimental study, the main 
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emphasis has been maintained on the robustness of the proposed observer against 
measurement noise and parameter variations. The use of PEMFC emulation system 
permits to conduct experiments on fuel cell auxiliary systems in real-time, while 
avoiding the risk of accidents (during worst case parametric variations) and cutting 
the consumption of expensive chemical reagents during fuel cell experiments that 
are not linked with fuel cell technology itself [12, 40, 46]. 


6.3.1 Algebraic Observer Design Using Sliding Mode 
Differentiators 


Algebraic observers are ideal for implementation in real-time embedded systems 
because of their low computational requirements. The exact definition of algebraic 
observability is given in Definition 6.1. Our objective is to design an algebraical 
observer for the oxygen and nitrogen partial pressures in the PEMFC air-feed system 
from the available measurements of supply manifold pressure and compressor flow 
rate. 

Let us briefly recall here that these observers are applicable to systems whose states 
can be expressed in terms of input and output variables and their time derivatives 
up to some finite degrees. Further details can be found in Diop et al. [10] and Ljung 
and Glad [34]. In this section, we will first demonstrate the algebraic observability 
of PEMFC air-feed system. Then, we will present Lyapunov-based adaptive HOSM 
differentiators for the implementation of the algebraical observer. Finally, using these 
differentiators, we will present the algebraical observer for this system. 


6.3.1.1 Algebraic Observability 


Definition 6.1 Consider the nonlinear system described by the following dynamic 
equations, 


x(t) = f(x), u(t), (6.36) 
y(t) = A(x(t)), 
where f(-,-) € R” and h(-) € R? are assumed to be continuously differentiable. 
x(t) € R" represents the system state vector, u(t) € R” is the control input vector 
and y(t) € R? is the output. 
The system in (6.36) is said to be algebraically observable if there exist two 
positive integers js and v such that 


x(t) = O(y, 3, ¥.-++ yu, a, ii, --- sul) (6.37) 


where ¢(-) € R” is a differentiable vector valued nonlinearity of the inputs, the 
outputs and their time derivatives [17]. 
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Let us consider the model of PEMFC air-feed system (4.29). Define the cathode 
pressure Peq aS a new variable X, i.e. X = x; + x2 + Cp. It follows from the last 
equation of (4.29) that 


X=y-~ 4 1 = ily. diy2). (6.38) 
“16 c14C16 {1 + cis [hye = 1] 


The time derivative of X is given by 


X= y2 " J 


“16 c14c16 {1 + C15 [() __ 1] 


Yi vega 
eversye?? (Gi) 


= 5 (6.39) 
ci C14C16 {1 +15 [() - 1] 
= p2(¥1, M1, Vi, V2). 
Then, in view of (6.38) and (6.39), the system states can be rewritten as 
¥ 1 c3(Y1 a C2) Wea,out pie - csp 
1= 2C5 6 — C51 
c4 — C5 L(C1 +8) (1 — Yi) — 2 — C7 
= 93 (91, 91, V1, Y2, ¥2), 

X2 = X — x, —c2 = pay, V1, V1, Ys 2), (6.40) 


y2 

x3 = — = 95()2), 
C17 

X4 = yi = Yo). 


It can be seen from (6.40) that all states have been expressed as functions of the 
system outputs and a finite number of their time derivatives, ie. y;, yj and yo. 
Therefore, according to the definition of algebraic observability 6.1, system (4.29) 
is algebraically observable. 


6.3.1.2 Adaptive Sliding Mode Differentiators 


Although system (6.40) is algebraically observable, it is necessary to estimate the 
time derivatives y,, y, and y2 accurately. In this work, y, and ¥, are estimated through 
an adaptive SOSM differentiator and y2 is estimated through an adaptive first order 
differentiator. As discussed in the introduction, gain adaptation have been designed 
to overcome the requirement of knowledge of upper bounds on higher order time 
derivatives. The structure of both these differentiators is presented in the following 
subsections. 
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Adaptive second order differentiator for y,(t): Suppose that the output yj (f) is a 
smooth function, consider the following third order system [30] 


—Ho(t) |zo — yi(t)| sign(zo — yi) + 21, 


in = 
zy = —Ai(t) zo — yi (C1 sign(zo — yi(t)) + 22, (6.41) 
Zo = —Ad(t)sign(zo — yi(t)). 


Denote oo = zo — y(t), 01 = Z1 — y(t) and oy = Z2 — y(t), and then system 
(6.41) can be rewritten as 


Gy = —o(t) lool} sign(oo) + 01, 
6) = —Ai(e) lool} sign(o) + 09, wie? 
0 = —)o(t)sign(oo) — Vit), 


where [1 (t)| < L3 and L3 is an unknown positive constant. The time varying gains 
do(t), Mt (t) and do (t) are designed as follows 


mitt = 2120), HO) = 15030), oO S110: (6.43) 
for some positive time varying scalar L(t) which is adapted according to 


: k, if 
Cee ae i Pol 8, (6.44) 
0, otherwise, 


where k > 0 is an arbitrary positive design constant. 


Proposition 6.1 Consider the error system (6.42). Suppose that the gains o(t), 
A, (t) and A2(t) satisfy (6.43) and (6.44). Then, the states of the error system (6.42) 
converge to zero in finite time, i.e., 0; =0, i € {0, 1, 2}. 


Proof Consider the following Lyapunov function candidate for the error system 
(6.42) 


_ = = Tye = \4 
VELM) = WO+ 7 (LO-L ie (6.45) 
where Vo(€) is based on the Lyapunov function proposed in Moreno [37] 
“ -3%2 0 


WO=ETE, F=|-jfx 2% 37], (6.46) 
0 323 3 
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where £7 = [lool’ sign(dp) o1 3], coefficients (V1, V12> Y2> V13> 23» 3) are chosen 


such that Vo(£) is positive definite and radially unbounded. Since the adaptation law 
(6.44) makes the adaptive law L(t) bounded, it means that there exists a positive 
constant L* such that Lit) < L* fort > 0. 

The derivative of the Lyapunov function candidate (6.45) is given by 


VE, L() = Vo) —k|L@ - L*|’. (6.47) 


Under conditions given in Theorem | [37] and for some positive constant 0, we yaks 


Vo(é) < V4 (€). Apply Jensen’s inequality |x| + |y| > (xl? + lvl)". ¢ = 4 
1, one gets 


VE, L) < —0V; —k|E@) — L*|? < —mino, 480". (6.48) 


According to the Theorem 4.2 in Bhat and Bernstein [6], it follows that € converges 
to zero in finite time, i.e.0; = 0, i € {0, 1, 2}. This completes the proof. | 


Adaptive first order differentiator for y.(t): Let us consider that the output y(t) 
is a smooth function. An adaptive first-order sliding mode differentiator can be con- 
structed for this signal as follows: 


€) = X(t) ler()|? sign(er(t)) — ker (t) + &, 


(6.49) 
&2 = —a(t)sign(e(t)) — kaei@), 


where e;(t) = €, — y2(t). Denote e2(t) = &) — y(t) and its dynamics can be rewrit- 
ten as 
a(t) = —A(Z) ler(1)/? sign(er(t)) — k\Mer(t) + &, 
€(t) = —a(t)sign(e;(t)) — ka(ter(t) — yr(t), 


(6.50) 


where |2(t)| < L2, Lo is an unknown positive constant and the time varying gains 
A(t), a(t), k(t) and k,(t) are formulated as 


Mt) = AVLO, alt) = aL), ex) 
k(t) = kyLQ),  kalt) = koL7(2), , 


for some positive constants Ao, Qo, ky), Ka, and a positive, time-varying, scalar 
function L(t) which is adapted according to 


k, if jei()| £0, 


Ln = . 
0, otherwise, 


(6.52) 


where k is an arbitrary positive constant. 
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Proposition 6.2. Consider the error system (6.50). Suppose that the coefficients Xo, 
ao, ky, and ky, in (6.51) satisfy the following condition: 


Aaigka, > 8ki, a0 + 9AGK;, - (6.53) 


Then, the states of the error system (6.50) converge to zero in finite time, i.e. e\(t) = 
e(t) = 0. 


Proof A new state vector is introduced to represent the system (6.50) in a more 
convenient form for Lyapunov analysis. 


C1 L2()ler(2)|2sign(ei(t)) 
C= |6]| = L(tei(t) (6.54) 
G3 e2(t) 


Thus, the system (6.50) can be rewritten as 


EAU g 
. Lit) 7 
a, t 
C= Tl AiG + L(t)A2¢ + m2 ; (6.55) 
—y2(t) 
where ; 
-¥%¥ 0 } -~ 0 0 
Aj = 0 —-A, 0}, Ar= 0 —k,, 1]. (6.56) 
—ay 0 0 0 —k,, 0 


Then, the following Lyapunov function candidate is introduced for the system (6.55) 


1 4ag + rN Aokyo — Xo 
VOsCR, Fs 5 dokr ky, 2k, —k), |, (6.57) 
— Xo —ky, 2 


where the matrix P is positive definite. Taking the derivative of (6.57) yields 


(cs 


gL it) 
V=-L 
eG 


L(t) 


+ (Tae) —jnaic + PC, (6.58) 


where 
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gn = [-ro ky 2], a2 =[0 G OJ, 


ds rN + 2a0 0 —Xo 
Q, = a O kag + 5k —3kry | 
—Xp —3k), 1 (6.59) 
ao+242 0 0 
Q, = ky, O kag +k, —krg 
0 —ky, 


It is easy to verify that (2; and {22 are positive definite matrices under the condition 
(6.53). 
Since Amin(P)||Cl]? < V < Amax(P)||C||?, Eq. (6.58) can be rewritten as 


. Ainin 2 1 Amin Q 
Vy 2=10)— (20 y} _ py Amin) y (6.60) 
Minox (P) a) 
L Lit 
4 alana ys 3 s ) cTac, 
Anin(P) 
where 
4ag + BY + Aokyo + ao 0 0 
Q= 0 roto 0 (6.61) 
0 0 agk2, + Aokry + 
In view of (6.61), (6.60) can be rewritten as 
1 Lit) 
Vs -“LMOn1 -— 2) V? -— | L@) 3 - UT Vs (6.62) 
where 
_ Amin (921) = Lo\Iqill2 
1 = TW _ a ae, 
Xinax(P) Minin P) (6.63) 
=3 Amin (22) _ Amax(Q) 
we Xmax(P) — DXmin PY: 


Given that L(t) > 0, thus the terms L(t) — y2 and L(t)y3 — inn are positive in 
finite time, it follows from (6.62) that 


V <—c,V? — eV, (6.64) 
where c, and cp are positive constants. By the comparison principle [24], it follows 


that ¢ converge to zero in finite time, i.e. e;(t) = 0 and e2(t) = 0. Thus, Proposition 
6.2 is proven. | 
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Propositions 6.1 and 6.2 imply that the states of the adaptive second order differ- 
entiator (6.41), (6.43), (6.44) and first order differentiator (6.49), (6.51), (6.52) con- 
verge to their real values in finite time, i.e. z} > y(t), Z2 > V(t) and & > y(t), 
respectively. 

The algebraical observer is formulated as the following theorem. 


Theorem 6.3 For any u € U such that y(t) and y(t) are continuously differen- 
tiable, the following dynamic system 


X1 = 93 (y(t), 21, 22, yo(t), &2), 

Xo = va (y(t), Z1, 22, yo(t), £2), 

#3 = gs(y2(0)), xeoe 
X4 = yo(yi(t)), 


where & and 21, Z2 are the estimates of yo(t), yi (t) and y(t) respectively, is a finite 
time observer for system (4.29). 


Proof From the results of Propositions 6.1 and 6.2, it follows that there exists a finite 
time ft, after that £5, z; and z converge to y(t), y; (t) and yj (t) exactly. Consequently, 
limx;=x;, i = {1,2,3,4}. a 
t>ter 

Remark 6.2 In practical implementation, ideal sliding mode is not achievable due to 
measurement noise and numerical approximation errors. Moreover, the adaptation 
laws in (6.44) and (6.52) will result in unbounded L(t) and L (t). A feasible alternative 
to overcome this disadvantage is to modify the adaptive laws by dead zone technique 
[48, 52] as follows 


k, if >e k, if > 
, if lool = €, igs” if je;(t)| >, (6.66) 


L@ = 
0, otherwise, 0, otherwise, 


where € and € are sufficient small positive constants. In this sense, the gains L(t) and 
L(t) will be adapted dynamically according to (6.66) and oo and e; (t) will be forced 
back into a real sliding mode regime in finite time. 


6.4 Simulation Results of SMO for PEMFC Systems 


A completely nonlinear model developed is used in our observer design, comprising 
of the following state vector: 


P 
x= [Po PH: PN; Psm, "sm, Pv,an» Pv,ca, Prms Tt | ’ (6.67) 


where py, is hydrogen pressure in the anode, ms», is mass of air in supply manifold, 
Pv.an 18 Vapor partial pressure in the anode, p,q is vapor partial pressure in the 
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Initial 
Condition 


Fig. 6.2 Schematic diagram of the proposed SOSM observer 
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Fig. 6.3. Stack current under load variation 


cathode, p,m is pressure of return manifold and T,, is the stack temperature. The 
nominal values of the parameters used in the simulation are shown in Table 6.1. 

The dynamic model is based on Pukrushpan et al. [43], with the added temperature 
model described by a lumped thermal model [9]: 


aT, Osc ~~ W.Cp, (Tse ~ Toan) 
dt mst Cp, 


’ 


(6.68) 


: T As 
Osou = Is AF + Vact + Ls: Ronm } 
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Fig. 6.4 Stack voltage response under load variations 


where m,; is the heat mass of the stack, C,,, and C,, are the specific heat, W, is the 
coolant flow rate considered as a control variable, T,.;, is the coolant temperature 
at the stack inlet, vat is the activation loss at catalyst layer and Dig, is the internal 
energy source. The latter is calculated as a function of the stack current, temperature, 
electrical resistance of stack layers Ro~m, Faraday’s number F and the entropy change 
As. The physical parameters were obtained through extensive experimentation. 

In order to obtain the numerical solution, the initial conditions of the state vector 
(6.67) are given as: 


Po,(O) = 0.10bar, py,(0) = 1.1 bar, 
Pn, (0) = 0.68 bar, Psm(O) = 1.48 bar, 
Msm(O) = 0.03kg, Pyan(O) = 0.41 bar, 
Pv.caO) = 0.63 bar, = pyrm(0) = 1.3bar, T;,(0) = 340K. 


(6.69) 


In order to test the proposed observer’s robustness against parametric uncertain- 
ties, the system parameters, i.e. volumes, anode inlet orifice, ambient temperature 
and stack temperature have been varied around their nominal values according to 
the worst-case percentages are given in Table 6.2. The schematic diagram of the 
proposed SOSM observer is shown in Fig. 6.2, where the observer gains are set as 
ro = 5, ao = 10 and L = 25, respectively. The initial errors of the states are set at 
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e1(t) 
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Fig. 6.5 Voltage measurement noise 
Table 6.1 PEMFC system parameters 
Symbol Value Symbol Value 
n 381 R 8.314 J/(mol-K) 
Ra 286.9 J/(kg-K) Patm 1.01325 bar 
Tamb 298.15 K F 96485 C/mol 
Mi 2 g/mol M, 18.02 g/mol 
Mo, 32 g/mol Ma 28.96 g/mol 
My, 28 g/mol Ry 296.8 J/(kg K) 
Ro, 259.8 J/(kg K) Ray 4124.3 J/kg K) 
Ry 461.5 J/(kg K) Psat(T fc = 80°C) 0.4707 bar 
Van 0.005 m3 Vea 0.01 m3 
Kann 2.1x 10-3 kg(Pa-s) | Kea,in 3.629 x 10-6 
kg/(Pa-s) 


20% of maximum deviation from the fuel cell system. For simulation purposes, the 
initial values were chosen as py, (0) = 0.9bar. 

During the tests, the load, i.e. the stack current shown in Fig.6.3. was varied 
between 100A and 200A. The stack voltage shown in Fig. 6.4 varies between 222 V 
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Fig. 6.6 Net power under load variations 


Table 6.2 Parameter variations of the PEMFC system 


Symbol Parameter Value (%) 
Tre Temperature of the fuel cell, +10 
2¢ 
Von Cathode volume, m +10 
Vsm Supply manifold volume, m? | —10 
Tamb Ambient temperature, °C +10 
Kan.in Cathode inlet orifice constant, | +5 
kg/(Pa-s) 


and 248V. Additional noise was added to the measured signal Vr, = Vee, + €1(t), as 
shown in Fig. 6.5. The FC net power varies with the load current (Fig. 6.3), between 
23.3 kW and 41.5 kW as shown in Fig. 6.6. The performance of oxygen excess ratio 
is given in Fig.6.7. It can be easily found that the oxygen excess ratio overshoots 
at t= 15 s with respect to the step decrease of load current (Fig. 6.3). This causes 
an overshoot in the stack voltage (Fig. 6.4). Figure 6.8 shows that the state variable 
(PH,) are well estimated by the proposed observer under step load variations and 
measurement noise. Figure 6.9 shows the response of observation error based on the 
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XO» 


Oxygen excess ratio 


Fig. 6.7 Performance of oxygen excess ratio 


measurements (Pan, Pan,ins PO»> 1s: and V;,). The initial errors of the states are set at 
15% of maximum deviation from the fuel cell system. It is clear from the figures, the 
proposed observer is robust against model uncertainties and the effect of the noise is 
essentially imperceptible. 


6.5 Conclusions 


In this chapter, two nonlinear observers, i.e., An SOSM observer and an algebraic 
observer based on adaptive sliding mode differentiators have been designed for 
PEMFC system. Firstly, with the help of the SOSM observer, the hydrogen partial 
pressure was successfully observed from the measurements of anode inlet pressure 
(Pan,in), anode total pressure (pan), oxygen partial pressure (po, ), load current (Js; ) 
and stack output voltage (V+,). It shows that the proposed observer is insensitive to 
model uncertainty and measurement noise. The effectiveness, robustness and fea- 
sibility of the SOSM observer under load variations, have been validated through 
simulation results. 
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Fig. 6.8 Estimation of py, under load variations and measurement noise 
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Fig. 6.9 Estimation error of hydrogen partial pressure 
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Secondly, the design of an algebraic observer has been presented. The system 
states, oxygen, nitrogen partial pressures, compressor angular speed and supply man- 
ifold pressure, are expressed in terms of input and output variables and a finite number 
of their time derivatives. These variables are estimated via adaptive differentiators 
without any knowledge of the upper bounds of their higher-order time derivatives. 
The oxygen and nitrogen partial pressures have been successfully observed from the 
measurements of stack current, supply manifold pressure and compressor flow rate. 
The effectiveness and feasibility of the proposed observer has been validated. In the 
next chapter, we will focus on the design of fault diagnosis approaches for PEMFC 
systems using sliding mode technique. 
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Chapter 7 ®) 
Sliding-Mode-Observer-Based Fault creek 
Diagnosis of PEMFC Systems 


This chapter presents a fault diagnosis method for PEMFC systems taking into 
account a fault scenario of sudden air leak in the air supply manifold. Based on the 
control-oriented model proposed in the literature, an adaptive-gain SOSM observer 
is developed for observing the system states, where the adaptive law estimates the 
uncertain parameters. 

The residual signal is computed online from comparisons between the oxygen 
excess ratio obtained from the system model and the observer system, respectively. 
Equivalent output error injection using the residual signal is able to reconstruct the 
fault signal, which is critical in both fuel cell control design and fault detection. The 
performance of the proposed observer is validated through an HIL simulator which 
consists of a commercial twin screw compressor and a real-time PEMFC emulation 
system. The experimental results illustrate the feasibility and effectiveness of the 
proposed approach for application to PEMFC systems. 


7.1 Introduction 


Modern complex industrial systems can not be operated safely without reliable fault 
diagnosis and isolation (FDI) schemes in place [55]. Such systems including PEM 
fuel cells are vulnerable to system failures or mechanical faults that can lead to 
catastrophic consequences. PEM fuel cells are electrochemical devices that convert 
the chemical energy of a reaction between hydrogen and oxygen into electricity. 
Among different kinds of fuel cells, PEM fuel cells are suitable for both stationary 
and automobile applications due to the ongoing development of PEM technology 
[4, 6]. However, these energy generation systems based on fuel cells are complex 
and several additional equipments are required to make the fuel cell work at the 
optimal operating point. For this reason, they are vulnerable to system failures or 
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mechanical faults that can cause the shutdown or permanent damage of the fuel cell. 
Thus, reliable FDI schemes for such systems are of great importance. 

FDI is usually achieved by generating residual signals, obtained from the differ- 
ence between the actual system outputs and their estimated values calculated from 
dynamic models. Such approach usually involves two steps: the first step is to decou- 
ple the faults of interest from uncertainties and the second step is to generate residual 
signals and detect faults by decision logic. Several practical techniques for these steps 
have been proposed in contemporary literature, for example geometric approaches 
[31], Ho.-optimization technique [16, 34], observer based approaches (e.g. adaptive 
observers [49, 52, 53], high gain observers [3, 48], unknown input observers [5, 
15]). However, in active Fault Tolerant Control (FTC) systems [55], not only the 
fault needs to be detected and isolated but also needs to be estimated such that its 
effect can be compensated by reconfiguring the controller [1, 8]. Hence, there is need 
of fault reconstruction schemes which estimate the fault’s shape and magnitude. The 
challenge here is that the FDI approach should be able to differentiate between the 
effects of faults of interest and from uncertainties. 

For linear systems, this challenge has been addressed by so called robust FDI 
schemes which are insensitive to model uncertainties [11]. In this respect, signifi- 
cant developments have been made and existing conditions have been given [28]. 
Qiu and Gertler [34] proposed residual generator design schemes based on the 
H,.-optimization technique. Then, the so-called H_/H,. approach was considered 
for the design of residual generator for uncertain systems in [16], where the effect 
of the uncertainty on the residual was taken to be upper-bounded and the effect 
of the fault to be lower bounded. Hou and Muller [15] used UIO to decouple the 
uncertainty from the FDI scheme. Qiu and Gertler [34] proposed residual genera- 
tor design schemes based on the H..-optimization technique. Then, the so-called 
H_/H, approach was considered for the design of residual generator for uncer- 
tain systems in [16], where the effect of the uncertainty on the residual was taken 
to be upper-bounded and the effect of the fault to be lower bounded. Compared 
with linear systems, nonlinear systems are much more complicated, which makes 
their study more difficult. Even though a generic solution to the FDI problem is still 
to be developed, a number of effective observer based FDI approaches have been 
reported. These include techniques based on disturbances decoupling [14], geometric 
approaches [31] and so on. 

Sliding mode technique is known for its insensitivity to external disturbances, 
high accuracy and finite time convergence. SMOs have been widely used for fault 
reconstruction in the past two decades. Edwards et al. [7] proposed a fault recon- 
struction approach based on equivalent output error injection. In this method, the 
resulting residual signal can approximate the actuator fault to any required accuracy. 
Based on the work of Edwards et al. [7], Tan et al. [44] proposed a sensor fault 
reconstruction method for well-modeled linear systems through the LMI technique. 
This approach is of less practical interest, as there is no explicit consideration of 
disturbance or uncertainty. To overcome this, the same authors [45] proposed an FDI 
scheme for a class of linear systems with uncertainty, using LMI for minimizing 
the L» gain between the uncertainty and the fault reconstruction signal. It should be 
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noted that only linear systems are considered in the above works and only few works 
have been reported for nonlinear systems. Jiang et al. [18] proposed an SMO-based 
fault estimation approach for a class of nonlinear systems with uncertainties. Yan 
et al. [50] proposed a precise fault reconstruction scheme, based on equivalent output 
error injection, for a class of nonlinear systems with uncertainty. A sufficient condi- 
tion based on LMI is presented for the existence and stability of a robust SMO. The 
limitation is that requires a strong structural condition of the distribution associated 
with uncertainties. Later, this structural constraint was relaxed by the same authors 
[51], where the fault distribution vector and the structure matrix of the uncertainty 
are allowed to be functions of the system’s output and input. However, these works 
require that the bounds of the uncertainties and/or faults are known. Although in 
the work of Yan and Edwards [51], the requirement on the bound of uncertainty is 
removed, but it still needs to know the bound of the fault signal. 

One of the main problems in the PEM fuel cell operation is the so-called oxygen 
starvation phenomenon during fast load variation. Accurate regulation of the oxygen 
excess ratio is required in order to avoid oxygen starvation [32]. This is a challenging 
task which comes from two aspects, on one hand, it is difficult to measure the oxygen 
excess ratio value, on the other hand, the fuel cell systems suffer from various faults, 
such as sudden air leak in the air supply manifold. Hence, from the point view of 
FTC, only FDI is not enough, the fault signal should be reconstructed and then its 
effect on the system performance can be compensated during active FTC design. 

During the last decades, different kinds of model based techniques have been 
widely studied in the areas of FDI, health monitoring and complex industrial sys- 
tems [38, 47]. Sliding mode based approach is one of the most attractive techniques 
due to its robustness against external disturbances, high accuracy and fast conver- 
gence [12, 23]. Several SMO-based FDI approaches have been proposed for linear 
systems [7, 45], but only few works have been reported for nonlinear systems, espe- 
cially for nonlinear uncertain systems [12, 51]. However, first-order sliding mode 
algorithms are employed in the above works that require low pass filters to generate 
the output injection signals and induce undesirable chattering effects. Furthermore, 
the employment of low pass filters will introduce some delays which results in inac- 
curate estimates or even instability of the system. In recent years, HOSM technique 
has been widely studied due to the reasons that it does not require any low pass filters 
while keeping all the good properties of the standard sliding mode [25, 27]. This 
technique can also be used to alleviate the chattering effect because of its continuous 
output signal. 

From the application point of view, several observer based FDI approaches have 
been studies for the fuel cell systems. Arcak et al. [2] developed an adaptive observer 
for hydrogen partial pressure estimation based on the fuel cell voltage. Gérgiin et 
al. [13] developed an algorithm for estimating partial pressures and the membrane 
water content in PEMFCs based on the resistive cell voltage drop. This algorithm 
has incorporated two adaptive observers for hydrogen and oxygen partial pressures, 
adapted from the work of Arcak et al. [2]. However, both of the above works lack 
robustness against the fuel cell voltage’s measurement noise and the internal model 
relies upon unmeasurable values. Ingimundarson et al. [17] proposed a model based 
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estimation approach to hydrogen leak detection in PEMFCs without the use of relative 
humidity sensors. Escobet et al. [9] designed a fault diagnosis methodology for 
PEMECs where the residuals are generated from the differences between the PEMFC 
simulator included with a set of typical faults and a non-faulty fuel cell model. More 
recently, LPV observer was proposed by Lira et al. [26] for the application to fault 
detection in PEMFCs, where the stack current was taken as the scheduling variable. 
It should be noted that most of the cited papers are based on model linearization 
around pre-defined operating points of the system, depending upon the operating 
conditions such as temperature, humidity and air flow. Moreover, the calculation of 
Jacobian matrix of complex models like fuel cells are time consuming, and therefore 
difficult in real-time implementation. Several efficient control strategies have been 
designed in order to track the optimum operating points [19, 21, 56]. In the aim to 
assure optimal operating conditions, state estimation, parameter identification and 
fault reconstruction are going to play an important role. The state observers serve 
as a replacement of physical sensors for obtaining unmeasurable quantities such as 
oxygen and nitrogen partial pressures in the cathode side. Parameter identification is 
important for online monitoring procedures. The need for complete reconstruction 
of the faults not only facilitates fault diagnosis but also plays an important role in the 
enhancement of system robustness properties. Furthermore, state, parameter and fault 
estimations might be used for control purposes to maintain the fuel cell in an optimal 
operating points [10]. As is known the major obstacle to the implementability of these 
controllers, is the absence of reliable measurements of oxygen and nitrogen partial 
pressures in the cathode side. Therefore, the problems of state estimation and fault 
reconstruction arise due to incomplete knowledge of the parameters and states of the 
system. The auxiliary elements (e.g., valves, compressor, sensors) make the fuel cells 
vulnerable to several faults [9]. Moreover, the models used in [26, 32] are obtained 
through a Jacobian linearization of the PEMFC nonlinear dynamic model around 
its pre-established optimal operating point. As is known that the fuel cell system 
exhibits highly nonlinear dynamics and a fuel cell’s efficiency is highly dependent 
on the various parameters such as operating conditions, such as temperature, humidity 
and air flow. A nonlinear four state dynamic model of the air-feed system proposed by 
Suh and Stefanopoulou [42] is chosen for state estimation and fault reconstruction. 
This model has been validated in the complete operating range of a fuel cell. A set 
of auxiliary elements (valves, compressor, sensors, etc.) are needed to make the fuel 
cell work at the optimal operating point. Thus, the problems of state estimation and 
fault reconstruction arise because of incomplete knowledge of the parameter and 
states of the system. 

In this Chapter, state estimation, parameter identification and fault reconstruction 
are studied simultaneously for a class of nonlinear uncertain systems with Lipschitz 
nonlinearities. This class of systems sufficiently defines the nonlinear dynamics of the 
PEMECs. The approach involves a simple adaptive update law and an adaptive-gain 
SOSM observer. The uncertain parameters are estimated and then injected into an 
adaptive-gain SOSM observer, which maintains a sliding motion even in the presence 
of fault signals. Finally, once the sliding motion is achieved, the equivalent output 
error injection can be obtained directly. Thus, the fault signals are reconstructed based 
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on this information [7, 41]. The proposed SOSM algorithm combines the nonlinear 
term of the STA and a linear term [30]. The behavior of the STA near the origin 
is significantly improved compared with the linear case. Conversely, the additional 
linear term improves the behavior of the STA when the states are far from the origin. 
Therefore, the proposed algorithm inherits the best properties of both the linear and 
the nonlinear terms. 

The proposed fault reconstruction scheme is evaluated by implementing on an 
instrumented HIL test bench, which consists of a real commercial twin screw com- 
pressor and a real-time fuel cell emulation system [29]. This considered model has 
been experimentally validated on a 33-kW PEM fuel cell in a wide operating range 
with less than 5% relative error [43]. A fault scenario, i.e., sudden air leak in the air 
supply manifold is considered. Its effect is simulated with an increment of supply 
manifold outlet flow constant, which presents a mechanical failure in the air circuit 
resulting in an abnormal air flow [9]. The system states (nitrogen partial pressure 
and compressor speed) are estimated successfully. The stack current, considered as 
an uncertain parameter, is estimated through an adaptive update law. This eliminates 
the need of an extra current sensor. The fault scenario considered in this study is 
a suddenly air leak in the air supply manifold, which results a change in the outlet 
air flow in the supply manifold. It is reconstructed faithfully through analyzing the 
information, which is obtained on-line from comparisons between the measurements 
from the sensor installed in the real system and the outputs of the observer system. 
The phenomenon so-called oxygen starvation is monitored through the estimation 
of a performance variable oxygen excess ratio. Since it will cause hot-spots on the 
membrane, resulting in permanent damage to the fuel cell if this ratio is less than 
1 for a long time. The robustness of this observer against measurement noise and 
parameter variations is also validated experimentally. 


7.2 Dynamic Modeling of PEMFC Systems 


A typical PEM fuel cell system is shown in Fig.5.1 which consists of four main 
subsystems, i.e., the air feed subsystem, the hydrogen supply subsystem, the humidify 
subsystem and the cooling subsystem. In order to achieve high efficient operation, 
a set of auxiliary elements (valves, compressor, sensors, etc.) are needed to make 
the fuel cell work at the optimal operating point. In this study, we will focus on the 
controller design of the air compressor for the air feed subsystem. The air compressor 
used to supply the oxygen to the cathode side is the core component and can consumes 
the power generated by the fuel cell up to 30% [46]. Therefore, efficient control of 
the air compressor is critical for the whole system and effects the system’s efficiency 
directly. A typical PEM fuel cell polarization curve is shown in Fig. 7.1. 

As is widely known, the dynamics of the PEM fuel cell system are with highly non- 
linearities. Therefore, suitable control-orient model taking into account the dynamic 
behaviors of the cathode partial pressure dynamics, the air supply manifold dynamics 
and the compressor dynamics is needed for the controller design. Some assumptions 
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Fig. 7.1 Typical fuel cell voltage 


are made to simplify the nonlinear model of the fuel cell system while keeping the 
dynamic behaviors of the air-feed subsystem [33]. Mainly, it is assumed that the 
temperature of the cathode inlet flow is regulated to a constant value through a heat 
exchanger. This is reasonable because the response time of the stack temperature is 
slow [22]. The relative humidity of both anode and cathode sides of the fuel cells are 
regulated to the desired relative humidity through an instantaneous humidifier. The 
hydrogen pressure in the anode side is regulated to follow the cathode pressure by 
the anode valve. Only vapor phase is considered inside the cathode and extra water in 
liquid phase is removed from the channels. The compressor motor current dynamics 
are neglected because the electrical time constant is very small as compared to the 
mechanical dynamics [36]. 


7.2.1 State-Space Representation of PEMFC Modeling 


In view of (4.29), and define the state variables x = [x,, x2, x3, x4]' = 
[P0.» PNr» Weps Psm\" and X = x; +x2+ 2. Then, the nonlinear dynamics of 
the PEMFC system are expressed by the following equations: 


t= F(x) +G-utw-é, (7.1) 


where 


7.2 Dynamic Modeling of PEMFC Systems 135 


by (x4 — X) — x1 fi (X1, 2) 
by (x4 — X) — x2 fi (41, x2) 
= by 
POO = | eyes — bt | (ss) 1] fae o | (7.2) 
falxa)[ favs, x4) — Die (x4 — X) | 


[00b30]", w = [-b,000]', 


G 


b 
where fi (%1,%2) ‘= grFbspthe Wea.ours £203, %4) = Wep and f3(x4) = big 


x4 \P12 
(1+ bis[(55)"" - 1) 
The stack current € := J,, is considered as the external disturbance and the control 


input uv := v,», is the motor’s input voltage. The outputs and performance variables 
of the system are given by: 


= [yi y2 ys] = le Wep Vals l= [z1 za] = (Pret Aa] (7.3) 


where P,,., and Ag, are fuel cell net power and oxygen excess ratio, respectively. 
The fuel cell net power P,,.; is the difference between the power produced by the 
stack P;,; and the power consumed by the compressor. Thus, the net power can be 
expressed as: 

Puer = Pst — Pep, (7.4) 


where Ps; = Is;Vs; and Pep = TcemWcp are the stack power and compressor power, 
respectively. The oxygen excess ratio \o, is defined by the following equation: 


d _ Wo,.in _ big (Ps = Pea) 
OO, = — 7 
7 Wo,,react Diol st 


(7.5) 
The parameters b;, i € {1,..., 19} are defined in Table4.3. More details on this 
model are available in [32, 42]. 


Remark 7.1 As is well known that a set of auxiliary are needed to make the fuel cell 
work at its optimal operating conditions. Particularly, due to it is essential to regulate 
the oxygen excess ratio close to 2 in the presence of fast load variation and fault 
scenario [37, 39]. The level of Ao, is critical because fast load demand will result 
sudden decrease of the oxygen flow rate. On one hand, once the value of oxygen 
excess ratio decreases to a critical value (normally less than 1), oxygen starvation 
phenomenon which leads to the FC degradation occurs. On the other hand, higher 
value improves the fuel cell stack power but also result in higher power consumption 
of the air compressor. Therefore, the problems of oxygen excess ratio estimation and 
fault reconstruction arise due to the reasons of safety and high efficiency. 
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7.3. Nonlinear Observer Design for PEMFC System 


In this Section, we consider the problem of an adaptive SOSM observer for a class 
of nonlinear systems, in which the uncertain parameter is estimated with the help 
of an adaptive law. Then, an SOSM observer with gain adaptation is developed 
using the estimated parameter. Finally, based on the adaptive SOSM observer, a fault 
reconstruction method which can be implemented online is proposed. 


7.3.1 Problem Formulation 


Consider the following nonlinear system 


xX = Ax+ g(x, u)+ o(y, u)O+ w(y, u) f(t), 
y=Cx, (7.6) 


A, Ao 
A3 Ag 
is the control input, y € VY C R? is the output vector, g(x, u) € R" is Lipschitz 
continuous, ¢(y, u) € R’*? and w(y, uv) € R”*” are assumed to be some smooth 
and bounded functions with p > g +r. The unknown parameter vector 0 € R? is 
assumed to be constant and f(t) € R’ is a smooth fault signal vector, which satisfies 


where the matrix A = | x € R” is the system state vector, u(t) € U Cc R” 


IfOl < pr, IFOM < pe, (7.7) 
where p), p2 are some positive constants that might be known or unknown. 
Assume that (A, C) is an observable pair, and there exists a linear coordinate 


I , 
: x= [zt zy’, with z} € R? andz. € 


transformation z = Tx = 
—Ha_—p)xp In—p 


R”~?”, such that 


Ai A12 


TAT oS 
° ie A22 


| where the matrix Ax, = Ay — HA, € R“-?)*-P) is Hur- 


witz stable. 
ecT!= [Zp 0]. where J, € R?*? is an identity matrix. 


Assumption 7.1 There exists a function w;(y, u) such that 
wily, u 
Tw(y,u) = ; (7.8) 


where w)(y, u) € R?*’. 


System (7.6) is described by the following equations in the new coordinate, 
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z= TAT '2+Tg(T'z,u)+To(y,wd+Twly, uw f(t), 
yecrs. (7.9) 


By reordering the state variables, the system (7.9) can be rewritten as 


y = Any t+ Ariza + 91 (Z2, yu) + Oily, WO t+ wily, wf), 
Zo = AnZ2 + Any + g2(Z2, y, u) ae oo, u)o, 
ee (7.10) 


where 


roo. = [20°] rare = [BEB], aay 


dr(-,-) :R? x R™ > R?, do(-,-) 2 R? x R™ > R"?, gi(-,,) RR? x RY? x 
R” > R?, (+): R? x R™? x R" > R™?, 


Remark 7.2 The Assumption 7.1 is a structural constraint on the fault distribution 
w(y, u). It should be noted that there are no such structural constraints on the uncer- 
tain parameter distribution ¢(., -). The structure of (7.10) is more general than that 
proposed in [51], since there exists Lipschitz nonlinearities g;(z2, y, U), g2(Z2, y, u) 
which depend on the states. 


7.3.2 Adaptive SOSM Observer Design 


Now, we will develop an adaptive SOSM observer for the system (7.10), which con- 
sists of an adaptive update law and an adaptive-gain SOSM observer, to estimate the 
uncertain parameter and system state variables, respectively. The basic assumption 
on the system (7.10) is as follows: 


Assumption 7.2 ([51]) There exists a nonsingular matrix T € R?*?, such that 


T [oy u) wily, w)] = hee PG ‘al (7.12) 


where ©; (y, u) € R2*7, o(y, u) € R’*” are both nonsingular matrices and bounded 
in(y,uyeyVxl. 


Remark 7.3 The main limitation in Assumption 7.2 is that the matrices [¢) (y, u), w1 
(y, u)] must be block-diagonalizable by elementary row transformations [51]. For 
the sake of simplicity, the case of only one fault signal and one uncertain parameter 
is considered (¢ = r = 1). 
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Let zy = T y, where T is defined in Assumption 7.2. Then, the system (7.10) can 
be described by 


gy = TAyy + T Ariz. +Tgi(y, 22,4) + las ° a Pe (y, I i 


Z2 = Agoz2 + Aory + go(Z2, y, u) + d2(y, u)O, 

y= T tps 
(7.13) 

where 
_ Aat = Asi Wg, (y, 22, U) 
PA Sl Tea qT. a1 
ul ba 12 bar g(y, 22,4) = be u) 

(7.14) 


Let define z, = ee Z»] > with z,, € R’, zy, € R’. Then, in view of (7.13) and 
(7.14), we can obtain 
= Ayy t+ Ariz2 + Wy, (y, 22,4) + Pi(y, WO, 
Zyy = Any oF Anz ag Waly, 225 u) ae Pr(y, u) f(t), 
Z2 = Aoiy + Ao2Z2 + go(y, 22, U) + Go(y, uO, 
_ T 
=f" [a,enl" 


Nt 
< 
| 


(7.15) 


| 


The adaptive SOSM observer is represented by the following dynamical system 


ex = Auy + An122 + Wa, (y, 22, u) + ®D (y, uo i [ley,), 
a = Apy + A222 + Wy,(y, 22, u) + Mey,), (7.16) 
Any + An2o + go(y, 22, u) + daly, 0), 


22 


where ju(-) is calculated by the SOSM algorithm 


Ls) = \e)|s|2sign(s) + a | sign(s)dt + k)(t)s + kat) | sdT, (7.17) 
0 0 


with the adaptive gains A(t), a(t), k\(t) and k,,(t) which will be determined later on. 
The observation errors are defined as e,, = Zy, Zi Oy = By — ly 2 = 


29, and 6 = 6 — 6. The estimate of 0, denoted by 6, is given by the following adaptive 
law 


6 = —K(y,u) (Any + Ande + Wo, (, 22,0) + i(y, wb ay), (7.18) 


where K (y, u) is a positive design matrix which will be defined later. 


Remark 7.4 It can be seen that the adaptive law (7.18) depends upon z,,. The 
adaptive-gain SOSM algorithm (7.17) can be used to estimate the time derivative 
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of zy, in finite time, without requiring the bound of |Z,,|. The differentiator has the 
following form 


—ov L(t)|Zo0 — zy, [2 sign(zo — Zy,) — ky L(t) (Zo — Zy,) + 21, 
—apL(t)sign(zo — 2y,) — ka? (t)(Z0 — Zy,)s 


- (7.19) 


Z 


where Zo and z, are the real-time estimations of z,, and z,, respectively. The param- 
eters of the differentiator Ap, k),, Qo, ka, and the time varying function L(t) are 
designed as (7.33) and (7.32), respectively. 


Subtracting (7.16) from (7.15), the error dynamical equation is described by 


6. = Aner + gly, 22, 22,u) + d2(y, w)8, (7.20) 

8 = —K(y,u)(Anier + Wa (y. 22. 22.0) + Pi(y, 08), (7.21) 
éy, = —pley,) + Aries + Gi(y, w)O + Wy, (y, 22, 22, W), (7.22) 
éy, = —p(ey,) + Arrer + Po(y, uw) f(t) + Wo, (y, 22, 22,4), (7.23) 


where 


RY, 22> Zo, u) = CAOr 22, u) _ nO, 20, u) 
Wo. (y, Z2, 22, U) = Wo, (y, 22, uw) — Wo, (y, 22, u) 
Wa, 22> 20, u) = Wa, 225 u) ~ WO, 20, u). 
Some assumptions are imposed upon the error systems (7.20)-(7.23). 
Assumption 7.3 The known nonlinear terms g2(y, Z2, u), Wg, (y, Z2, uw) and W,, (y, 


Z2, u) are Lipschitz continuous with respect to z2 


lg2(y, Z2,U) — g2(y, 22, u) || < Qellz2 — Zell, 
I| Wo, (y, Z2, Uw) — Wo, (y, 22, wll < YqllZ2 — Zell, (7.24) 
|W, £25 u) —_ Wo, 20, u)|l < Yq |IZ2 — Zoll, 


where g,, Yg, and y2 are some known Lipschitz constants for W,,(y, Z2,u), 
W,,(y, Z2, u) and go(y, Z2, u) respectively [54]. 


Assumption 7.4 Assume that the Huritwz matrix Az) satisfies the following Riccati 
equation 
ApPi + PiAnt+ BPP t+ 2+e)ln-p = 0, (7.25) 


which has a symmetric positive-definite solution P,; for some ¢ > 0 [35]. 


Assumption 7.5 Assume that the positive design matrix K (y, u) satisfies the fol- 
lowing equation 
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K(y, w®\(y,u) + OF (y,w)K'y,u) — 7%, KO, WK, u) = lg, (7.26) 


for some ¢« > 0. 


Remark 7.5 The Assumption 7.3 means that the nonlinear terms g2(y, Z2, u), Wg, (y, 
22, u) and W,,(y, Z2, u) are Lipschitz with respect to the variable z2. The Assump- 
tion 7.4 means that the Hurwitz matrix A satisfies the Riccati equation (7.25) for 
some small ¢ > 0. The Assumption 7.5 means that the design matrix K (y, u) must 
be designed such that K(y, u)®)(y, uv) + bl iy, u)K"(y, u) is positive definite in 
(YyueVxu. 


Now, we will first consider the stability of the error systems (7.20) and (7.21). 


Theorem 7.1 Consider the systems (7.20) and (7.21) satisfying the Assumptions 
7.3—-7.5. Then, the error systems (7.20) and (7.21) are exponentially stable, if for any 
(y,u) € Y x U, the following matrix 


_ AT pT 
Q1 = late. nT ayy el as _ ° (7.27) 
is positive definite. 
Proof A candidate Lyapunov function is chosen as 
Vi(er,0) = es Pier + 00 (7.28) 


and the time derivative of V; along the solution of the system (7.20) and (7.21) is 
given by 


Vi < ef (ADP + PiAn + PIP) + 2In—p)er 
+ 2e5 (Pida(y, u) — AL, K'(y, u)) 6 — «6'6 
=— [ef 67] O 4 ; (7.29) 
Hence, the conclusion follows given that Q, is positive definite in (y,u) € Y x U. 
a 


Assumption 7.6 It is assumed that ||®;(y, u)||, ||G2Cy, u)|| are bounded in (y, u) € 
Yxu. 


Remark 7.6 Theorem 7.1 shows that lim e2(t)=0 and lim A(t) = 0. Conse- 
too too 


quently, the errors e2, 6 and its derivatives é9, @ are bounded. Under Assumptions 
7.3 and 7.6, the time derivatives of the nonlinear terms in the error dynamics (7.22) 
and (7.23) are also bounded 
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=X, 


@ A T A 
| at (Arre> + Bo(y,u) f(t) + We (y, 22, 22, u)) 


(7.30) 


ad p- are P 
| ry (Aries + @i(y, ujyO+ Wa, (y, 22,22; u)) = x2, 


where x; and x2 are some unknown positive constants. 


In what follows, the objective is to prove the finite time convergence of the system 
(7.22) and (7.23). 


Theorem 7.2 Assume that (7.30) holds and the adaptive gains X(t), a(t), ky(t) and 
k(t) in the SOSM algorithm (7.17) are formulated as 


AG) = AVEO, a) = aL), 
Bl) =] EQ), mt) = kal @, 


(7.31) 


where Xo, Qo, ky, and ka, are positive constants and L(t) is a positive, time-varying, 
scalar function. The adaptive law of the time-varying function L(t) is given by 


k, if |ey,| 40, i = {1, 2}, 


Lit) = 
© 0, else, 


(7.32) 


where k is an arbitrary positive value. Then, the trajectories of the error system (7.22) 
and (7.23) converge to zero in finite time, if Ao, Qo, ky, and Ka, in (7.31) satisfy 


Aagka, > 8ki,a0 + 9AGK},- (7.33) 


Proof We only consider the proof for the system (7.23) and the proof for the system 
(7.22) takes the same way. The system (7.23) can be rewritten as 


is 
éy, = —X(t) ley, |? sign(ey,) — key, + ¥, 
p = —a(t)sign(ey,) — ka(tey, + o(t). 


(7.34) 


Under the condition (7.30), it follows that || o(t)|| < 2, where v2 is an unknown 
positive constant. 

A new state vector is introduced to represent the system (7.34) ina more convenient 
form for Lyapunov analysis. 


G1 L3(t)ley, |2sign(ey,) 
C=1Q]= Lies (7.35) 
G3 y 


Thus, the system (7.34) can be rewritten as 
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» Lit) = Og =) D0 ZHI 
C= a7] 0 -rA0/C+LM)] 0 -k, 1/6+)] 4.6]. (7.36) 
ICi| y KG) 
—ap 0 0 0 —keay 0 o(t) 


Then, the following Lyapunov function candidate is introduced for the system (7.36) 


4ao + AG dok ro —o 
VO =¢'PC, P = 5 Nokyy KX, + 2kag —kro | 5 (7.37) 
—Xo —k, 2 


where the matrix P is symmetric positive definite due to the fact that its leading 
principle minors are all positive under the condition (7.33). 
Taking the derivative of (7.37) yields 


; TQ Lit 
V = —-Li(t) (é is F ca) + ogg + ns ) pe, (7.38) 
ICil L(t) 
where gq = [—Ao —ky, 2). n= [1 e 0] and 
d rN + 2a0 0 —Xo 
Q = = 0 kay + 5k, —3ka, |, 
=\y <3, 1 
5 (7.39) 
ag + 25 0 0) 
2 = ky 0 Kay +2, — ka, 


0 —k, 1 


It is easy to verify that £2; and {22 are positive definite matrices under the condition 
(7.33). Since Amin(P)||Cl]? < V < Amax(P)||C||?, Eq. (7-38) can be rewritten as 


Amin ($22) Xallaille ,2 L(t) 


: Amin 2 1 
(ent ye Dy ge pe 
Nhax(P) ae) 
(7.40) 
where 
4ag ++ Aokyy + 0 0 
o= 0 roto 0 (7.41) 
0 0 2kayk?, + Aokyy + 


Using (7.41), (7.40) becomes 


V <-(LQn -—) V2? - (LO - ge) Vv; (7.42) 
L(t) 
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where 
Amin (2 ) Amin (2 ) Amax ( ) 
poe Ae alla, ney : 4 = ¢ - 743) 
imax (P) Minin (P) max ( ) min ( ) 


Because Lit) > 0 such that the terms L(t)y; — y2 and L(t)y3 — ut are positive 


in finite time. It follows from (7.43) that 


V < -c\V? eV, (7.44) 
where c and cp are positive constants. By the comparison principle [20], it gives that 


¢ converge to zero in finite time, i.e. ey, = 0 and é,, = 0. This completes the proof. 
a 


Remark 7.7 In view of practical implementation, the condition |e,,| = O in (7.32) 
can not be satisfied due to measurement noise and numerical approximations. In 
order to make the adaptive algorithm (7.32) practically implementable, one has to 
modify the condition |e,,| = 0 by dead-zone technique [40] as 


k, if jey|>=T, 


Lit) = (7.45) 


0, else, 


where 7 is a sufficiently small positive value. 


Theorems 7.1 and 7.2 have shown that systems (7.16) and (7.18) are an asymptotic 
state observer and uncertain parameter observer for the system (7.15) respectively. 
In the next Section, we will develop the fault reconstruction approach based on those 
two observers. 


7.4 Fault Reconstruction 


The fault signal f(t) will be reconstructed based on the proposed observer by using 
an equivalent output error injection which can be obtained once the sliding surface 
is reached and maintained on it thereafter. 

It follows from Theorem 7.2 that e,, and éy, in (7.23) are driven to zero in finite 
time. Thus, the equivalent output error injection can be obtained directly 


pley,) = Aner + Boy, u) f(t) + Wo (y, Z2, 22, 4). (7.46) 


From Assumption 7.2, ®2(y, uw) is a bounded nonsingular matrix in (y, uv) € Y x U, 
then the estimate of f(t) can be constructed as 


FQ) = ®7'(y, Wey). (7.47) 
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Theorem 7.3 Suppose that conditions of Theorems 7.1 and 7.2 are satisfied, then 
Ff (@) defined in (7.47) is a precise reconstruction of the fault f (t) since 


A 


lim If) — fOl =O. (7.48) 
Proof It follows from (7.46) and (7.47) that 


IFO — FON = lez, w(Arer + Woe) I 


- = 4 (7.49) 
S ||P, (y, w) Aralllleall + Yq, lO. w)illleall. 
It follows lim || f(t) — fll = 0 given that lim |le2|| = 0. Hence, Theorem 7.3 is 
too t>oo 


proven. | 


7.5 FDI of PEMFC System 


Figure 5.1 shows a block diagram of a typical PEM fuel cell system in automotive 
applications, which consists of four major subsystems: the air feed subsystem, the 
hydrogen supply subsystem, the humidify subsystem and the cooling subsystem. The 
PEMEFCs are supplied at the cathode with compressed air via a supply manifold. The 
compressor draws the air from the atmosphere directly. At the cathode exit, the air 
enters an outlet manifold that is open to the atmosphere. Several assumptions are 
made for the system. Mainly, the anode pressure is assumed to be well controlled 
and equal to the cathode pressure [33]. The input reactant flows are humidified in a 
consistent and rapid way and the high pressure compressed hydrogen is available. 
Vapor partial pressure in the stack is considered equal to the saturation pressure. 

In view of the system (7.1), the state space representation of this model can be 
written is as follows 


: €3(X1 — C2) 
y= (cy - Ca)(X] x4) Wea,out c7€, 
KX] 
; C3X2 
x2= cg (X4 _ x1) = — Wea,out» 
KX] 
ee gee (7.50) 
X3 = —C9X3 1 Wep + ¢13u, 
X3 Ci 
: X4 \ C12 
x4 = C441 +cjs (=) —llyx [ Wer — c16(x4 — x)|, 
ul 

where x := Pcq is the cathode pressure, x7 := py, is the nitrogen partial pressure, 


X3 = Wep is the compressor speed and x4 := Psm is the supply manifold pressure. 
Wea,out 18 the cathode flow rate which is a function of the cathode pressure Wea. our = 
Kea,out/ Pca — Patm and W,, is the compressor flow rate which is a function of the 
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angular speed of the compressor W.p = h3(x3) = Cj7W¢p [29]. The stack current 
€ is considered as an uncertain parameter 6. The control input wu represents the 
motor’s quadratic current component. The cathode and supply manifold pressures are 
assumed to be available for measurement. Thus, the system outputs are y = [x 1 xa)’. 
The system performance variable, oxygen excess ratio Ag, is defined as the ratio 
between the oxygen entering the cathode Wo, ;, and the oxygen reacting in the fuel 
cell stack Wo, react: 


Wi in _ 
Sig es Cig (x4 x) (7.51) 
7 Wo), react C19 


Due to the reasons of safety and high efficiency, it is typical to operate the stacks 
with this value equals 2 during step changes of load current demand [29]. It should 
be noted that positive deviations of Ag, above 2 imply lower efficiency, since excess 
oxygen supplied into the cathode will cause power waste, and negative deviations 
increase the probability of the starvation phenomena. 

We consider a fault scenario: a sudden air leak in the air supply manifold. This 
fault is simulated with an increment Acj¢ in the supply manifold outlet flow constant 
C16 := Ksm.our, Which is translated into a change in the outlet air flow in the supply 
manifold Womour = (C16 + ACi6)(X4 — X1) [9, 26]. We assume that this fault appears 
after time t = 50s, i.e., Acig = 0.2ci¢6. Thus, the fault signal f(t) appears in the 
output channel x4 is defined as 


ine la x (x4 — x1) ke/s, if t > 50s, a 


0, else. 
In order to design the proposed observer for the fuel cell system. Let us define 
Zy, =X, Zyy = X4, Z2 = [x2 x3] and 6 := €. Then, system (7.50) is described as 


the form of (7.15) 


C3 (Zy, = C2) Wea,out 


Zy, = —(c1 + cg) (Zy, Zyy) c4O, 
Ky, 
: Zy, \ 12 
Zy = C14 1+ C15 (- :) -—1 x [#3(D2z2) — (c16 + Ac\6) (Zy, aa ah 
11 


03 Wea,out 
ae —-H 0 ea —Cg Cg as Diza(H ae KZyy ) 
SW, ey 0 0]? — Diz [(2)" - 1 h3(D2Zz2) + ¢13u 
SEE oo —> -—S D222 cS 
———_———_— > 


ci 


g2(Y,Z2,u) 


(7.53) 
where 
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D = [1, 0], Dy = [0, 1], Pi(y, u) = —-C4, Pr(y, u) = C5, 
€3(Zy, — C2) P(Zy,) 


KZy, 


Zy, \ 12 
Woo (ys 200M) = crgcis( =) (Hta(Daz2) — cis(@ — zy); 
11 


Wa 0, Z2, U) = ’ o2(y, u) = 0, 


and the design parameter H is chosen to satisfy the Riccati equation (7.25). The fault 
signal is weighted by cs, modeled as 


oe: {! Teas (2)" = i]}. (7.54) 


The adaptive SOSM observer for the system (7.53) is designed as the form (7.16) 
and (7.18) 


C3 (Zy, = C2) 
KZy, 


zy, —(c =F Cg) (Zy, Zyy) Wea,out c40 =F Ley, ), 


A Zy, \O2 A 
o cu {t+eis| (2) = 1] x hares) — ealen = 25)] + Hen) 


Cil 


22 = AgzZ2 + Ary + go(y, 22, 4), 
(7.55) 
€3(Zy, —€2) Wea, out : 
Rly, + éy) ’ where Cy, = fy 


and @ = —K (cd + (c1 + €8)(Zy, = Zy.) + 
—Zy,, Cy, = Zy, — Zy, and the adaptive-gains of the SOSM algorithm ju(ey,), 1(ey,) 
are designed according to (7.31), (7.32). The value of z,, is obtained from the robust 
exact finite time differentiator (7.19) in [24]. The oxygen excess ratio Ag, and the 
fault signal f(t) are estimated 


te = cise =m) #6 = Hey) (7.56) 
C190 C5 


Remark 7.8 The Assumptions 7.3 and 7.6 are satisfied by the functions W,, (y, Z2, “), 
Wo (y, 22, U), Pi(y, u), O2(y, u) and g2(y, Z2, u). The Riccati equation in Assump- 
tion 7.3 is satisfied by appropriate value of the design gain H = 0.5. The Assumption 
7.4 is also satisfied for some € = 0.02, since the Eq. (7.26) is simplified into a scalar 
equation. 


7.6 Experimental Results 


The structure of the observer based fault reconstruction strategy is shown in Fig. 7.2. 
The nominal values of the parameters for the HIL emulator are shown in Table 4.4. 
In order to test the robustness of the proposed observer based fault reconstruction 
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Motor inverter I 
a 


Twin screw 
compressor 


Exact robust 
differentiator 


Adaptive 
SOSM settsteesscecceeeeeeeeed 
_| observer 


Adaptive real-time observer 


Fig. 7.2. Schematic diagram of the observer based fault reconstruction 


approach against parametric uncertainty, the system parameters, i.e. temperatures, 
volumes, cathode inlet and outlet orifice, motor constant and compressor inertia have 
been varied around their nominal values according to the worst-case percentages, are 
given in Table 5.2. 


0.04 1 1 1 1 1 
Noise e(t) 


Measurement noise (bar) 


0 0.05 0.1 015 O02 025 03 O35 O04 045 0.5 


Time (s) 


Fig. 7.3. System output measurement noise 
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During the tests, the stack current was varied between 100 and 450 A, correspond- 
ing to flow rate variation in the compressor between 7 and 28 g/s. The measurement 


Observation of py, (bar) 
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Time (s) 
(a) Estimate of pn, 


Observation error (bar) 
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Fig. 7.4 Estimate of nitrogen partial pressure and its error 
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noise was also included to test the robustness of the proposed approach, that is, 
Vt = Pea + E(t), Yo = Psm + €(t), Where e(t) represents the measurement noise 
shown in Fig. 7.3. The state observations (py,, W¢p) in response to steps changes in 
the load current are shown in Figs. 7.4 and 7.5. It can be seen from these figures that 
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Fig. 7.5 Estimate of compressor speed and its error 
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the settling time of both upwards load changes (i.e. t = 40s, t = 60s) and down- 
wards load changes (i.e. t = 80s, t = 100s) is less than 2s. It should be pointed out 
that the estimate of compressor speed is compared with experimental measurement. 

The outputs of the adaptive differentiator and Levant’s fixed gain differentiator 
[24] are shown in Fig. 7.6. It is easy to see that the proposed adaptive differentiator 
has faster convergence and less overshoot. Figure 7.7 shows the estimate of the stack 
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Fig. 7.6 Performance of the adaptive differentiator 
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Fig. 7.7 Estimate of stack current and its error 


current, as an unknown parameter @, based on the adaptive law. Both, the result 
obtained from application of Levant’s differentiator and the adaptive differentiator 
are presented. It can be seen that the adaptive differentiator improves the performance 
of the adaptive law and gives a good estimate for the stack current. The fault signal 
is reconstructed faithfully as shown in Fig. 7.8. It is clear that the proposed scheme 
is capable of reconstructing fault signal and state estimation simultaneously in the 
presence of uncertain parameter. The estimate of oxygen excess ratio is shown in 
Fig. 7.9a. It is easy to see that this value decreases after t = 50s. This is due to the 
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Fig. 7.8 Fault reconstruction and its error 
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occurrence of the fault (air leak in the supply manifold) after t = 50s. In this case, 
the air flow supplied by the compressor needs to increase, taking into account the 
effect of the fault, in order to ensure safe operation of the fuel cell. The time history 
of the adaptive gain L,(t) is shown in Fig. 7.9b, where the exponential convergence 
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Fig. 7.9 Estimate of Ao, and gain L(t) versus time (s) 
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of the observer is ensured. The gains of the SOSM algorithm stop increasing when 
the output observation error converges to zero in finite time. 


7.7 Conclusion 


This Chapter has presented a robust fault reconstruction method for a class of non- 
linear uncertain systems with Lipschitz nonlinearities based on an adaptive SOSM 
observer. An adaptive update law has been given to identify the uncertain parame- 
ter. The estimated parameter is then injected into an adaptive-gain SOSM observer, 
which maintains a sliding motion in the presence of the fault signal. The proposed 
fault reconstruction approach has been successfully applied to the PEMFC system 
considering a sudden air leak in the supply manifold. The experimental results have 
shown that the proposed approach is effective and feasible. 
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Chapter 8 M®) 
Sliding Mode Control of DC/DC Power sheet 
Converters 


In the previous chapters, we have successfully designed adaptive HOSM based 
observers for state observation and FDI of the PEM fuel cell systems. We now turn 
our attention towards the power side of the PEMFC-. In fact, the PEMFC itself has 
severe dynamic limitations due to the time response of fuel flow and fuel delivery 
systems (hydrogen and air feed systems) [23]. In order for a PEMFC power system 
to be employed in varying load applications like electrical vehicles, storage elements 
with fast response time need to be integrated into the system. A typical PEMFC 
power system usually relies on rechargeable batteries, super-capacitors or both for 
improved power dynamic characteristics in transient high-power demands. 

Evidently, such a hybrid power system requires power conditioning circuits with 
precise power control algorithms behind them such that the output power is compat- 
ible with the constraints of the power bus. These circuits include DC/DC converters 
which are used to control the voltages of the fuel cell and storage elements to a 
desired value, rectifiers and inverters that are used to convert AC to DC and vice 
versa. Therefore, precise control, and therefore observation, are equally important in 
this area of the PEMFC system as they were in the air-feed system. 

In this chapter, we have studied the control and observation of DC/DC power 
converters that are employed in fuel cell hybrid power systems, using sliding mode 
techniques. The main focus is on the necessary modification and improvement of 
conventional sliding mode methods for their applications to fuel cell power systems. 


8.1 Introduction 


As a ubiquitous direct current (DC) voltage conversion technology in industrial appli- 
cations, DC/DC converters have been widely applied in the uninterruptible power 
supply, telecommunication equipment, DC motor drives, power grids, and clean 
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energy systems [7]. DC/DC converters generally include boost type [26], buck type 
[31], and buck-boost type [18] converters. In DC/DC boost converters, one DC volt- 
age can be converted to another higher DC voltage through a converter [2]. Nowadays, 
this converter is widely used in remote and data communications, computers, office 
automation equipments, industrial instruments, military, aerospace and other fields 
[3, 11]. For an accurate voltage regulation of the DC/DC boost converters, the output 
voltage of the converters to be controlled is always a challenging problem due to the 
existence of various uncertainties. 

To represent the nonlinear systems with complex uncertainties, the interval type-2 
(IT2) fuzzy models [10] were proposed, which are established by applying type-2 
fuzzy sets. They allow uncertain parameters or unknown variables to exist in the 
membership functions of the fuzzy models. This modeling methodology has sub- 
stantially extended its utilization on the uncertain systems, and it has also been 
demonstrated that the IT2 fuzzy model is less conservative than the type-1 fuzzy 
model based on type-1 fuzzy sets [9]. Nevertheless, the defect of this method is the 
lack of self-learning and self-organizing ability. What’s more, fuzzy neural networks 
(FNNs) [1] can make up for the deficiency, and many online learning algorithms 
have been studied for the updating of the IT2FNNs [6, 12, 13]. To mention a few, 
the authors in [6] developed an IT2FNN with support-vector regression to deal with 
noisy regression problems. [13] simplified a type of self-evolving IT2FNN for reduc- 
ing the time-consuming K-M iterative steps, by using a gradient descent algorithm. 
For the problems of noise cancellation and system modeling, a Takagi-Sugeno-Kang- 
type-based self-evolving compensatory IT2FNN was designed in [12]. This approach 
combines the ability of IT2 fuzzy reasoning to manage uncertainties and the capa- 
bility of artificial FNNs to learn from processes. Hence, one of the motivations of 
this work is the utilization of the IT2FNNs to the representation of a more practical 
boost converter model with uncertain parameters. 

For the control of the DC/DC converters, traditional proportional-integral-type 
or proportional-integral-derivative-type controllers are the common control method 
[19] in practice. But the existence of the uncertainties in the DC/DC converters 
remains the difficulty to the stability analysis of the corresponding tracking error 
systems with such a control. Recently, intelligent control methods like fuzzy control, 
adaptive control, backstepping technique control and SMC, have attracted consid- 
erable research interests in dealing with the control problem of the DC/DC boost 
converters. What needs to be emphasized is that, among these control methods, 
SMC technique [15, 16, 25, 30, 32] has been recognized as a very effective and 
popular method due to its capability to achieve wide stability range, fast dynamic 
response, strong robustness, and simple control design procedure. Benefited from the 
SMC strategy, significant robust control methods have been developed for the voltage 
regulation of the DC/DC boost converters over the past few decades. For example, 
a robust adaptive sliding mode controller was designed in [20] using state observers 
for the boost converter with an external input voltage load and unknown resistive. 
The work in [21] presented an adaptive SMC approach to the DC/DC converters, of 
which the state-space model was formulated with some unknown coordinate com- 
ponents of the system equilibrium. The authors in [26] designed an SMC scheme 
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for the robust voltage control of the DC/DC converters. Then, by using the FNN, 
which is a type-1 FNN with Gauss-type membership functions, an integral SMC 
law was designed in [27] for the DC/DC converters. However, it requires a complex 
calculation of the adaptive parameters of the FNN for the updating of the SMC law, 
and the alleviated chattering actually still exists in the designed SMC effort due to 
its discontinuity. 

Inspired by the above discussion, in this paper, we will design a continuous SMC 
law, that is, a dynamic SMC, for the DC/DC boost converters, by utilizing a sim- 
plified IT2FNN. In terms of the representation of the uncertainties, we will fully 
formulate the parameter perturbation of input inductor, output capacitor, load resis- 
tor, input voltages, and noises in the DC/DC converters. The unknown part of the 
formulated uncertainties is estimated by a simplified ellipsoidal-type ITZ2FNN. Then, 
the dynamic SMC law will be designed based on both artificial FNN and adaptive 
IT2FNN. The adaptive rules of the IT2FNN are designed using the gradient descent 
method (GDM). Comparative simulations are detailed to illustrate the effectiveness 
of the designed control scheme. 


8.2 Modeling of Boost Type DC/DC Converters 


A typical boost DC/DC converter is considered, which transforms the DC fuel cell 
stack power to the external power devices with possible voltage-current requirements. 
Boost converters are used for applications in which the required output voltage is 
higher than the source voltage. From a control point of view, the boost converters are 
difficult because its standard model is a non-minimum phase system. Traditionally, 
the control problems of the DC/DC converters are solved by using pulse-width modu- 
lation (PWM) techniques. Sira-Ramirez [22] demonstrated the equivalence between 
SMC and PWM control in the low frequency range for a boost converter. In view of 
practical implementation, SMC is much easier than a PWM control since the maxi- 
mum frequency of commercially available switching elements increases higher and 
higher [24]. The sliding mode approach is expected to become increasingly popular 
in the field of power converter control. 

Figure 8.1 illustrates the circuit topology of a common boost converter, where D 
is the output diode, S is the main switch, and V;, L, R and C are input voltage, input 
inductor, the load resistor and output capacitor, respectively. u denotes the duty cycle, 
which is the percentage of the pulse width d over the total period 7 between an ON- 
and-OFF cycle of the switch. This is to say u = d for some predefined pulse width. 
The switching ON or OFF is determined by the PWM-based duty cycle. Generally, 
the parameters V;, L, R and C are usually treated as ideal constants. But, in practice, 
these parameters should be treated as variables due to the existence of measurement 
errors or changes in actual application requirements. By applying the state-space 
averaging method to the DC/DC boost converter, the average model is formulated as 
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Fig. 8.1 Topology of the 
DC/DC boost converter: a 
Switching OFF state, and b 
Switching ON state 


== Vi, (8.1) 
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where V, denotes the average output (capacitor) voltage, iz, stands for the average 
inductor current, u € (0, 1) is also called the average control effort here, and the 
considered ijq denotes the current source, which is caused by unmodeled dynamics, 
load variations, and unpredictable uncertainties. 

Further, considering the occurrence of the uncertainties caused by parameter per- 
turbations, a practical average model of the studied converter is further formulated as 


4 (1-u) 1 
= V, V; + AV,), 8.3 
a ag Tay G3) 
_, iy 1 iid 


(8.4) 


oO 


= ie V,+——_., 
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where AL, AC, AR and AV; are the perturbations of L, C, R and V;, respectively. 


8.2.1 Representation of the Tracking Error System with 
Uncertain Parameters 


Define the following tracking error variables: 


€i = let — IL, (8.5) 
ev = Veet — Vo, (8.6) 
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where i;e¢ 1s a reference input of the inductor current and Vier is the command of the 
desired output voltage. The dynamics of the tracking error system is thus derived as 


é= Aet+but+O+W, (8.7) 


e=[eev]', b=Bx,x(Q=[ic vo)’, 


_[O-;] g_[ 0 ae | al = [2] 
a(t v]8=[-20 0 | 2 =a] Y=] we | 


LAV; — (Vi — Vo) AL 


t)= 
q(t) L(L+AL) i 
(1) = —L(RAC + CAR + ARAO) Ve ipAC in 
PV RC(RC + RAC+CAR+ARAC) C(C+ AC) C+AC’ 
Viet — V; . Vo lref . 
t) = —— ref > Hn= Sa ref + 
wi(t) 7 + iret, W2(t) RC ct Veer 


In this paper, we assume the allowed parameter perturbations satisfy —aL < AL < 
aL, —bC < AC <bC, and —cR < AR <cR, where a, b and c are prior known 
constants belonging to (0, 1). Considering the uncertainties in B, we represent the 
uncertainty terms by the utilization of the fuzzy model approximation [4, 14] as 
follows: 


2 2 


B=B+AB=Y°Y hyi(Abyhcj(AC)(Bii + Bcj)+ AB, (8.8) 
i=1 j=1 


where —aL < AL = aLrand <aL and —bL < AC = bCrand <bL with 0 < 
a <1 and 0 <b <1 denote the estimates of AL and AC, respectively, used in 
this paper, B is the approximated matrix variable that will be used in further design 
procedure, AB is the approximation error, and 
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pee ae hig =1-hni, 
Fimax — Simin 
fomax — fr 
fei = omax = fomin’ her —— her: 


Remark 8.1 It should be noted that these two membership functions are not the 
component of the fuzzy engine, but only the approximation or estimation of B by 
using A Land AC chosen in this paper. Suppose the approximation error A B bounded 
by || AB|| < @. 6 can be evaluated by repeated simulations, and it has been validated 
in the simulation that || A B|| < 0.5} |B || no matter how the interval uncertainties AL 
and AC change. Based on the approximation (8.8), the error dynamics (8.7) can be 
rewritten as 

@= Ae+(B+ AB)xut+ O+W. (8.9) 


8.2.2 A Simplified Ellipsoidal-Type IT2FNN 


In order to fully formulate the uncertainties existing in the system dynamics, we 
employ an ITZ2FNN model to represent the tracking error system (8.7). For some 
input variable z = [Z, Z2 Z3 J' =[i, V, uJ", the IF-THEN rules are presented as 
follows: 

© Fuzzy Rule i: IF z; is $1;, and z, is S>; and z3 is S3;, THEN q is [g‘, q¥], 
where S$); is an IT2 fuzzy set for the antecedent part, g is the consequent variable 
represented by a weighting interval set [q/, ar |, andi = 1,2,...,r is the number 
of the fuzzy rules. Based on the fuzzy rules above, the ellipsoidal-type IT2FNN is 
detailed as follows. 

> Input layer: The j-th input is z; of the input layer for the IT2FNN (j = 1, 2, 3). 

> Membership function layer: For the j-th input and i-th node of the member- 
ship layer, we adopt the modified ellipsoidal-type of lower and upper membership 
functions [8], which is a kind of exponential-type type-2 membership functions with 
simplified computation complexity and more robustness [29], as follows: 


dé 
_ | xezmij : if |2x; = mij| < Vij, 
on (z;) = vi] (8.10) 
, else, 
i 
icmle 
isc) = — [Par fs Hf by — mul < vp, (8.11) 
, else, 


where 0 < a < 1, mj; and v,;; denote respectively the center and the width of the 
membership function. a, chosen as 0 < a < | in this paper, is the parameter to 
adjust the shape of the footprint of uncertainty (FOU) of the exponential-type type-2 
membership functions. 
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> Rule layer: Based on the lower and upper membership functions in Membership 
function layer, the i-th node in this layer is the product result for the interval degree 
of firing [A es h;| calculated by 


3 - 30k 
n= TL .9;. aS IT _, 5- ore 


> Type-reduction layer: This layer plays a role of reducing the type-2 fuzzy set 
to a type-reduced one for a crisp output [17]. Considering the consequent interval 
weighting g = [q‘, qV], we assign the bounds g/ and q¥ in ascending order, i.e., 
qe <qk <-+- <q andq? < q¥ <--- <q”. Considering the computation com- 
plexity, we use a simplified type of reduction expressed by 
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Via 21 a 
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Further, in order to guarantee h? < h, we update h? := min {h?, hi} and hY := 
max {hi, ee h and h are the firing strength vectors. To avoid singular values in 
actual computation, we here set h’ = Oif )-/_, 4, = 0, andhY = Oif )y_, h; = 0. 

> Output layer: we take the output of the IT2FNN from the average of m,; and 
My: 


hA [nv ny,...,n2] eR, avs 


r 


es a (8.13) 


8.2.3 Control Objective 


For the voltage regulation of the considered boost converter, in the following, we 
aim to firstly design a dynamic SMC against the uncertainties existing in the DC/DC 
model, based on the represented tracking error dynamics (8.9). For a given command 
of the desired output voltage Vier, the capacitor voltage V, will accurately track 
Ver equivalently to minimize the tracking error of the output voltage by properly 
manipulating the switch, which is determined by some designed duty cycle input wu. 
This is to be achieved despite variations in the considered model system. Meanwhile, 
the output voltage is to be controlled to tracking the command as fast and with as little 
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overshoot as possible during transients. Furthermore, the uncertainties, including the 
variations of the load resistance and input voltage, will be eliminated by using the 
introduced IT2FNN, and then an IT2FNN-based SMC law will be designed for the 


compensation of the uncertainties, as well as the robust tracking control of the desired 
inductance current and capacitor voltage of the DC/DC converter. 


8.3 Control Design and Stability Analysis 


8.3.1 Dynamic SMC Design 


8.3.1.1 Establishment of the Sliding Surface 


Recalling the factO < u < | andthe relation V; = V,(1 — u), in terms of the dynam- 
ics (8.9), we propose the following linear sliding surface: 


s=ge+h(u—vo)=gethv (8.14) 
where 7 
vsu—(1- 7") =u-wel-L th, 
and vp = 1 — a € [0,1], g=[91 »] € R'*? is a row vector, and h ER is a 


nonzero constant. Further, we consider V; = p and Veet =o as some known con- 
stants. Considering the sliding motion of s, that is 5s = s = 0, we obtain the following 
dynamic equivalent control law: 


Veq = —h-'g[Ae + (B+ AB)xu + Q+ W]. (8.15) 


Remark 8.2 The sliding surface (8.14) is designed based on the desired final equi- 
librium determined by the desired voltage and the nominal voltage. Besides the 
tracking error variable e, the introduced auxiliary variable v is also to be designed 
to approach to zero. Then, the accurate voltage regulation can be achieved while the 
sliding surface s = 0 is reached. 


Due to the existence of the uncertainties in Veg and the requirement of fast response, 
we introduce the following switching law: 


v, = —nsign(s) — ss — ¢|s|’sign(s), (8.16) 


where 77 > 0 is a gain to be determined, ¢ > 0, and ¢ > 0 and0 < 7 < 1 are prede- 
fined constants. vs is actually an exponential rate plus power rate reaching law [5], 
which is designed to properly and effectively reduce the reaching time of sliding 
mode and weaken the chattering caused by the switching term. 


8.3. Control Design and Stability Analysis 165 


Hence, based on veg and (8.16), we propose firstly the following dynamic SMC 
law: 


b = —h-'g(Ae + Bxv + Bxvup + W) — nsign(s) — ss — Cls|"sign(s), (8.17) 


with 7 = |x|], ¢ > 0, and ¢ > 0 and 0 < 7 < 1 considering || AB|| - ||x|| - |u|] + 
||Q|| < yl|x||, and y can be calculated with the interval uncertainties, and thus we 
omit the expression of 7 here. According to the proposed SMC law for the tracking 
error system (8.9), let us analyze the reachability of the desired sliding mode in the 
following theorem. 


Theorem 8.1 Consider the dynamics (8.9) and (8.17), and the sliding surface (8.14). 
For some given initial conditions u(0) and e(0), and¢, ¢ andT in (8.17), the trajectory 
of e can be forced onto the sliding surface s = 0 in finite time by the dynamic SMC 
law. 


Proof Let us check the reachability condition ss < 0, Vs 4 0, using the SMC law 
(8.17). For all s € 0, from (8.9), (8.14) and (8.17), it follows that 
ss = s(ge+hv) 
= sg| Ae + (B + AB)xu+ (B + AB)xvp + O+ w] 
—sg(Ae+ Bxv+ Bx vo + W) — ns - sign(s) — a ¢s - |s|"sign(s) 


es? tis 7 20, 


IA 


for s 4 0, where ¢ > 0,¢ > O and 7 > 0. Hence, the desired sliding mode is reach- 
able and can be achieved in finite time under the control law (8.17). This completes 
the proof. | 


When the sliding surface s = 0 is reached, the following sliding motion dynamics 
of e and v can be derived from (8.9) and (8.17). 


e=(A- vBye — Bev+ (ABxu+ Q)+ Bx ret (v + vo) + W, 
b= —h-g [oa — upB)e — Bev + Bxyes(v + v9) + w] 


8.3.1.2 The Existence of the Sliding Motion 


On the basis of the resulting dynamic system S, the following theorem states the 
existence of such a sliding motion S. 


Theorem 8.2 Consider the dynamic system S. For some given nonzero constant h 
and row vector g, if there exists matrix 0 < P = PER 0<Fr=F7 «RR 
and scalar } > 0, such that fori, j = 1,2, 
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Ps[A — v9(Bri + Bcj)] + [A — vo(Bii + Bcj)| Ps + <0, (8.18) 


where P; = P + dh~?g"q, then system S is boundedly stable. 


Proof Firstly, from condition (8.18), recalling hp; + hz2 = 1 and hc; +hc2 = 1, 
one can obtain that 


0<F <M& P(A — vB) + (A — vB)" Ps. 
Since s = s = 0 in the sliding stage, we obtain 
v=—h! ge, 
from (8.14). Then, the system S can be further expressed by 


s’: l @é=(A—wB)e+ +[B + AB)xv + ABxvg + Q|+ + (Bx; +W), 
“lv =A 'g[(A- vo Bye + Bxv + (Bx rer Vo + W)]. 


In terms of the dynamic system S, we choose the Lyapunov function candidate 
V =e'Pe + dv’ with matrix P > 0 and scalar 5 > 0. From the expressed é and 
in S’, it follows that 


V = 2e' Pé + 2v60 

= 2e"P {(A — vpB)e + [(B + AB)xv + ABxw + Q] + (Bx er + W)| 
—2e"g"h"5| hog [oa — wB)e + Bxv + (voBxrer + w)| | 

= 2eT P(A — wB)e + 2e™P {[B + AB)xv 
+ ABxuy + Q) + (Bx:eev0 + w)| 4 2e gh" g(A — upB)e 
13a! yh g'g|Bxv + (vo BX ret + w)| 

= Be (Pa oe =agniedt 2e"| P[B + AB)xv + Q] 
+qh~g'gBxv + (P + qh-7g"g) (vo BX ret + w)| 

= 2eT(P + gh-2g'g)(A — Bye + 2e7 {P| a + AB)xv + ABxw + Q 
+qh?P-'vg" g Bx] + (P+ qh-°g"g)(voBXret + w)| 

= 2e'Me + 2e'[Pq+(P+qh’g'g) w| 


—Amin(I") - llell? + 2llell (Pll - llall +P +.4h-7g"gIl- wll) 
—Amin(")|ell (ell — 7). 


IA IA 
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from which one can conclude that V <0 for e€ B°(r), where B‘(r) is the 


complement of the closed ball B°(r) centered at e=O with radius r= 
2(IPll-llg+i(P+94-79"g)Il-Ilwll) 
Amin") 


, where 


q = (B+ AB)xv + ABxvy + O+qh?P |g" gBxv, 
w= Vo BX rer + W. 


Obviously, r is bounded, since the uncertainty in g is bounded and w is known. 
Meanwhile, r can be minimized by increasing I" and h, and decreasing 6 and g. 
Hence, the sliding motion (8.14) exists based on condition (8.18) and is boundedly 
stable. It completes the proof. | 


8.3.2 Ellipsoidal-Type IT2FNN-Based Dynamic SMC Design 


In this subsection, we consider the uncertainties of which the bounds are unknown by 
the representation of the introduced ellipsoidal-type IT2FNN. Define q = [a qQ2 |’. 


Then, based on the considered IT2FNN, g; and q2 represented by gq, = 5qi Nk + 


dx, € R denote the unknown parts in voltage and current dynamics, including the 
; T . -T7T 
noises, where q, = Cine Gx] ER” hy, = [h,. h,| € R7’, and 5, denotes the 


approximation error. Thus, in the following, we will design some proper adaptive laws 
to estimate the unknown parameters g,. Denote the estimation error 4, = q, — Qx- 
Thus, the dynamics in (8.9) can be further written as 


é = (A— wB)e + Bev + (ABxu + O) + Bxrer(v + v0) + W 
= (A—wB)e + Bev + q(x, u) + Bxyepu + W, (8.19) 


where q(x,u) = ABxu + Q. According to this concept, we modify the dynamic 
SMC law from (8.17) in the following theorem, for the reachability of the desired 
sliding mode. 


Theorem 8.3. Consider the dynamic system (8.19) and the IT2FNN in Subsection 
II-B. The trajectory of e and v can be forced onto the sliding surface s = 0 in finite 
time by the following IT2FNN-based dynamic SMC law: 


2 
. is a 1 x : 
v = —h-'g(Ae+ Bxv + Bxuy + W) — 5 Se Gib — 58 — Cls\sign(s), 
k=1 
(8.20) 
with adaptive laws 


A 


m 1 . 
9x = — REG + 59s € R”, (8.21) 


where ¢ > 0, and ¢ > 0 and rT > 0 are some properly selected scalars, q,,(0) = 0, 
Kg > O and matrix E, € R*?". 
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Proof For the adaptive dynamic SMC law (8.20), we consider the Lyapunov function 
ie ee 
L= 55 + 2 sy GE d- 
k=1 


Then, from (8.14), (8.19) and (8.20), it follows that 


2 2 
L=si- Da FK a = s(gé + hi) — 2 EEG Ede + 59x Ehs) 
= sgl (A _ Bye — Bev + O + BX ret (V + uo) + w] _ sg( Ae + Bxv 
: fy ‘ ae 
+up Bx + w) 3 2 saudi — ss” — Cs|s|"sign(s) 

= 
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= yy Gi (— Keg, + 7H EKhks) 
k=l 


2 2 
rT ol i + -TA 
=s[am|[a 2] - 5 So song a — 587 — C181 + eed 
k=1 k=1 
1 2 
5D HG hes 
k=1 
2 1 1 2 2 
=e RRC TLL + 8) — 5 Yo gdp es — 587 — Cis +S) ede 
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1 2 
75 > IG Nes 
k=1 
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= 5(9151 + 922) — ss* — Cs|*7 + Do rede de 
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Meanwhile, for the vector @, = 4, — 4, considering 29,4, < —l@xll? + llaull? 
[28] and the fact that s - (g,d, + g2d2) < 0-87 + 7 - (910, + 9262)” with 0 > 0, we 
can further derive L as follows: 


tN. 
A 


2 2 
1 y willgell? +230 

< Os7 + 9 imo 4+@éa) — es? — Clg? = 5 Kellgll” + 5 Kellgell” 
k=1 k=1 


2 2 
1 . 1 
—(5 — 6)s* — ¢|s|'*7 — 5 ) relldgll” + 5 ) Kelldell’s 
k=1 k=1 


8.3. Control Design and Stability Analysis 169 


Adaptive IT2FNN-based Dynamic SMC 
Se AM Sa ca iehacenl e ear cse E 
l 


Astificial IT2FNN-based Dynamic SMC. |: | V, is 


[| zoe || 

ymca 
: F Shding \.!- 

| Dynamic SMC law urface A 

| on Surlace q References 

Ge 7... | - 

4 Vier |, 


Adaptive rules  -@---- | Initialization ‘| 


Fig. 8.2. Structure of the designed adaptive ITZFNN-based DSMC 


from which one can always find solutions to (¢ — @) > 5 and a > that is 


Sey 
tosay L < —pL + dwith y > Oand d = —<Is|!*7 + } («lai I? + K2llgall?) holds 
for any s 4 0. Since ¢ is bounded, the bounded stability of L can be thus ensured. 
This completes the proof. | 


Remark 8.3 Since the condition of s = 0 (in sliding stage) and q, (0) = O implies the 
dynamic SMC law is same with the one in Theorem 8.2, the analysis of the existence 
of the desired sliding motion is similar to Theorem 8.3, and is omitted here. 


Summarizing from the designed control laws based on the measurements 
{iz, Vo, u} and reference inputs {irer, Vier}, the structure of the designed adaptive 
IT2FNN-based dynamic SMC is illustrated in Fig. 8.2, in which the output uw is the 
desired duty cycle of the converter in Fig. 8.1. Further, in order to design an adaptive 
IT2FNN for a wider utilization on the elimination of various uncertainties, such as 
different unknown perturbations caused by input voltage V;, input inductor L, the 
load resistor R, output capacitor C, and noises, we design some adaptive rules for 
the updating of the parameters m;;, vj; and a;; of the designed type-2 FNN. 


8.3.3 Design of the Parameter Updating Rules 
of the Ellipsoidal-Type IT2F NN 


In terms of the antecedent parameters m;;, vj; and a;; used in each fuzzy rule of 
the IT2FNN, we apply the GDM to update m;;, v;; and a;; for the utilization in the 
flexible estimation of the uncertainties. Let m,;, vj; and a;; be expressed by €. The 
following gradient descent algorithm is designed for the parameter updating rules: 


E(k + 1) = €(k) - (8.22) 


de 
Lb 0&’ 
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where € = se (e here denotes e; or ey), oe denotes a partial derivative of the function 
€ of several variables € (mj,;, vj; and a;;). Meanwhile, we use the following updating 
law, that is, 


de 1 (Om, Oh, 99; Omy Oh; O45; 


86 2 °\ Oh, 3g, GE * Dh, Ty BE 


Oc 
to update ae? where 


Om, Gui — ML (Gui. h;) Omy — Gui — Mu (Gui. hi) 


Oh; 7 Via & "Oh; a ii , 
Oh; 3 Oh; 3 7 
dg. Leroy Gi agi; ere Gil: 


0. 


and 3 and oe are provided in Table 8.1. 


Remark 8.4 Actually, the designed parameter updating rules of the ellipsoidal-type 
IT2FNN can be utilized to a wider ellipsoidal-type IT2FNN with more inputs, 
although there are only three inputs to the IT2FNN. The number of the fuzzy rules 
is related to the number of the fuzzy sets chosen for each input z;. More fuzzy sets 
imply more complex computation, but more precise tracking control of the voltage. 
Therefor, a trade-off between the selection of fuzzy sets and the accurate voltage reg- 
ulation is recommended. In the simulation subsequently, we use two types of type-2 
fuzzy sets for each input to achieve the voltage regulation. 


Og.. ae 
Table 8.1 a and ee of the ellipsoidal-type IT2 membership functions 
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8.4 Simulation Results 


In this study, we numerically simulate the voltage regulation of the considered con- 
verter by using the average model (8.1) and (8.2). The uncertainties in the model 
caused by the perturbations of the input inductor and output capacitor the, as well as 
the variations of the input voltage and load resistance, are fully considered in sim- 
ulations. This section will present simulation results by three cases, i.e., (i) control 
against the input voltage variation, (ii) control against the load resistance variation, 
and (iii) control against both the variations of the input voltage and load resistance. 
For a comparative illustration of the proposed control method, the designed dynamic 
SMC, artificial ITZ2FNN-based dynamic SMC, and adaptive IT2FNN-based dynamic 
SMC are applied for each of the three cases. 


8.4.1 Set Up 


According to the configuration of the DC/DC boost converter, we take the nominal 
values of the model parameters from Table 8.2. 

In simulations, we set the sampling time at T = 50 ws. Besides, the prototype 
of the converter is considered with a nominal output power P, = 480 W, while the 
high-voltage bus is set at V.eg = 120 V for some inverter utilization. Accordingly, the 
reference input of the inductor current is designed as ipep = Kp) + Kj J ey t+ Kae} 
with Kp = 0.12, &; = 3310 and Kg = 0.005 [27]. Based on these setups, we use the 
designed control laws to illustrate the effectiveness and robustness of the proposed 
voltage tracking control method in different cases, subsequently. 


8.4.2 Case 1: Control Against the Input Voltage Variation 


Since the nominal value of the input voltage is Vi, = 60 V, we might consider that 
V; is varied between 48 < V; < 60 and then between 60 < V; < 72 with the load 
resistance R = 30 Q (P, = 480 W) being fixed. A varied V; is depicted in Fig. 8.3. In 
the following, we use respectively the dynamic SMC law (8.17), the IT2FNN-based 


Table 8.2 Nominal values of the parameters of the DC/DC boost converter 


Parameters Values Description 

Ch 860 LF Output capacitor 

fs 20 kHz Switching frequency 
In 860 wH Input inductor 

Rn 30 Q Load resistance 

Vin 60 Vv Input voltage 
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Variation of the input inductor V; 


45 1 1 1 1 
0 0.02 0.04 0.06 0.08 0.1 


Time (s) 


Fig. 8.3. Variation of the input inductor V; 


dynamic SMC law (8.20) with artificial parameters m,;, vj; and a;;, and ITZ2FNN- 
based dynamic SMC law (8.20) with adaptive parameters mj;, vj; and aj; to regulate 
the output voltage of the DC/DC boost converter. 

In this case, specifically, according to Theorems 8.1 and 8.2, the parameters of 
the dynamic SMC law (8.17) are set to 


g= [0.0037 0.0052], h = 850.0, ¢ = 0.0470, ¢ = 0.10, 7=0.10. = (8.23) 


Then, according to Theorem 8.3, the parameters of the ITZ2FNN-based dynamic SMC 
law (8.20) are designed as (8.23) with updated ¢ = 0.0470. Besides, for each input 
z; of the FNN, we use two types of IT2 fuzzy sets, i.e., PS (Positive Small) and 
PB (Positive Big). The corresponding membership functions g;;(zj) with FOUs are 
shown in Fig. 8.4, where the artificial parameters m;;, vj; and a;; for the membership 
functions 9; (z;) and gj;(z;) in (8.10)-(8.11) are chosen as: 


™M({1,2,3,4}1 = 3.9583, M{5,6,7,8}1 = 6.4583, 
M(1,2,5,62 = 78, M3,4,7,8}2 = 102, 
m(1,3,5,7)3 = 9.3, M12,4,6,8}3 = 9.7, 

vi) = 1.8750, vji2 = 18, 

U3 = 0.3, Qit = 0.5, 


a2 = 0.4, aj3 = 0.6, i =1,2,...,8. (8.24) 
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(a) 1 


Fig. 8.4 MFs of IT2FNN for inputs z,, z2 and z3 


Further, we design the adaptive IT2FNN to obtain the ITZ2FNN-based dynamic 
SMC law (8.20). The controller parameters are tuned as (8.23) with updated ¢ = 
0.0528. The parameters m;;, v;; and a;; for the lower and upper membership functions 
are updated according to the adaptive rules in Sect. 8.3.3. The initial conditions of 
the parameters are set as the ones in (8.24). We omit the computational process of the 
adaptive parameters, which are calculated by the gradient descent algorithm (8.22). 

For simplicity, we denote the three control approaches, i.e., the dynamic SMC 
law (8.17), artificial ITZ2FNN-based dynamic SMC law (8.20) and adaptive IT2FNN- 
based dynamic SMC law (8.20), respectively, the Control A, the Control B and the 
Control C. In this case, these resulting duty cycle inputs u are depicted in Fig. 8.5a. For 
the Controls B and C, the parameters of the adaptive laws (8.21) are set to K, = 0.1, 
K2 = 0.2 and FE; = E = I\6. Figure 8.5b shows the trajectories of the corresponding 
sliding surface s (8.14). According to these three control laws under the considered 
case, the regulations of the capacitor voltage V, are shown in Fig. 8.6. Figure 8.7 
illustrates the tracking errors of the capacitor voltage. Besides, the inductor current 
and its tracking error are respectively shown in Figs. 8.8 and 8.9 for the three control 
approaches. From the simulation results, we know that it illustrates the robustness 
of the three control algorithms to the variation of the input voltage, although no 
information of the uncertainties is needed to the IT2FNN-based dynamic SMC law 
(8.20). 
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Fig. 8.5 Duty cycle and sliding surface under Case | by the three controls 
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Fig. 8.6 Capacitor voltage under Case 1 by the three controls 
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Fig. 8.7 Capacitor voltage tracking error under Case 1 by the three controls 
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Fig. 8.8 Inductor current under Case 1 by the three controls 
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Fig. 8.9 Inductor current tracking error under Case | by the three controls 


8.4.3 Case 2: Control Against the Load Resistance Variation 


In this case, the load resistance R is supposed to be varied between 5 < R < 123 Q 
(200 < P, < 800 W), and the input voltage is set at V;, = 60 V. A piecewise variation 
of R is depicted in Fig. 8.10. Also, we use the three control approaches to verify the 
robustness to the load resistance variation. The parameters of the three controllers 


are given as follows: 
g= [0.0037 0.0052], h = 800.0, ¢ = 0.335, ¢ =0.15, 7 =0.10. = (8.25) 


For Control C, we update ¢ = 0.345 and use the same initial conditions as presented 
in (8.24) for the IT2FNN. The resulting duty cycle inputs u are shown in Fig. 8.1 1a. 
Figure 8.11b depicts the corresponding curves of the three sliding surfaces. 

By these three control laws, the regulations of the capacitor voltage V,, in this 
case, are illustrated in Fig. 8.12, and the tracking errors of the capacitor voltage are 
shown in Fig. 8.13. Figures 8.14 and 8.15 show the inductor current and its tracking 
error, respectively. Obviously, the robust control of the DC/DC boost converter with 
the load resistance variation performs well by the designed control approaches. 
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Fig. 8.11 Duty cycle and sliding surface under Case 2 by the three controls 
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Fig. 8.12 Capacitor voltage under Case 2 by the three controls 
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Fig. 8.13 Capacitor voltage tracking error under Case 2 by the three controls 


8.4 Simulation Results 179 


21 
E i, by Control A 
19 F = = =7,, by Control B 
rises d, by Control C 
@ 1) 1 : 1 
oS 
2 
a 
2 11 
= 
oO 
fH 
2 
Q OF 
5 
EI 
5 
3 
i i 


0 0.02 0.04 0.06 0.08 0.1 
Time (s) 


Fig. 8.14 Inductor current under Case 2 by the three controls 
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Fig. 8.15 Inductor current tracking error under Case 2 by the three controls 
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8.4.4 Case 3: Control Against the Variations of the Load 
Resistance and Input Voltage 


Generally, we consider the perturbations from both the load resistance and input 
voltage in this case. The variations of V; and R are used exactly as the ones in Cases 
1 and 2. The parameters of the dynamic SMC law are given as: 


g = [0.0037 0.0052], h = 800.0, s = 0.0468, ¢= 0.10, r= 0.10. (8.26) 


Moreover, we update ¢ = 0.0558 for Controls B and C for good tracking perfor- 
mance, and the initial conditions for the adaptive IT2FNN are the same as the ones 
used in the first two cases. 

Accordingly, the resulting duty cycle inputs u are shown in Fig. 8.16a. Figure 8.16b 
depicts the curves of the sliding surfaces. In this case, the regulations of the capacitor 
voltage V, are shown in Fig. 8.17, while Fig. 8.18 depicts the tracking errors of 
the capacitor voltage. Figures 8.19 and 8.20 illustrate the inductor current and its 
tracking error, respectively. Evidently, the three control scheme is effective to the 
voltage regulation of the DC/DC boost converter. Moreover, from the simulation 
results of the three cases, it can be concluded that the designed adaptive IT2FNN- 
based dynamic SMC scheme is more suitable without chattering in the control efforts 
for the voltage regulation of the converter, despite with unknown information of the 
variation of the system parameters. 
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Fig. 8.16 Duty cycle and sliding surface under Case 3 by the three controls 
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Fig. 8.19 Inductor current under Case 3 by the three controls 
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Fig. 8.20 Inductor current tracking error under Case 3 by the three controls 
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8.5 Conclusion 


The problem of robust voltage regulation of the DC/DC boost converter has been 
solved via the IT2FNN-based SMC in this paper. A modified linear sliding surface has 
been designed, which is depending on the desired voltage and the capacitor voltage. A 
new dynamic SMC against the complex uncertainties has been designed. Moreover, 
artificial and adaptive IT2FNN-based dynamic SMC laws have been synthesized. The 
updating rules of the parameters of simplified ellipsoidal-type membership functions 
have been designed based on the GDM. In the end, simulation results have verified the 
effectiveness of the presented adaptive IT2FNN-based SMC method. To accurately 
estimate the uncertainties existing the dynamics of the converter, some disturbance- 
observer-based control techniques can be considered for the voltage regulation of 
DC/DC boost converters. 
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Chapter 9 M®) 
Sliding Mode Control of Three-Phase sheet 
Power Converters 


This chapter presents an ESO-based SOSM control for three-phase two-level grid- 
connected power converters. The presented control technique forces the input cur- 
rents to track the desired values, which can indirectly regulate the output voltage while 
achieving a user-defined power factor. The presented approach has two control loops. 
A current control loop based on an SOSM and a DC-link voltage regulation loop 
which consists of an ESO plus SOSM. The load connected to the DC-link capacitor 
is considered as an external disturbance. An ESO is used to asymptotically reject this 
external disturbance. Therefore, its design is considered in the control law derivation 
to achieve a high performance. Theoretical analysis is given to show the closed-loop 
behavior of the proposed controller and experimental results are presented to validate 
the control algorithm under a real power converter prototype. 


9.1 Introduction 


Over the last decades, power converters have experienced a tremendous development 
in industrial applications, mainly due to their increased reliability, their advantageous 
properties such as high efficiency and power capacity, the reduced cost of power- 
electronic devices, and the increased application demands in industry [17]. Among 
the power converters, the three-phase PWM power converters are most commonly 
used, whose applications can be found in motor drives, energy storage systems, 
integration of renewable energy sources, etc. [33, 37]. 

Three-phase PWM converters play a key role in industrial applications like inte- 
gration of renewable energy sources, energy storage systems, motor drives, etc [3, 
26]. Particularly, active front ends (AFE) are grid-connected converters that offer 
features as bidirectional power flow, near-sinusoidal currents, and power factor and 
DC-link capacitor voltage regulation capability [35]. For this reason, the control 
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objectives for this application are: (1) maintain the DC-link voltage regulated to a 
certain reference, and (2) supply a desired reactive power and draw grid currents 
with the lowest possible harmonic distortion. 

To make the AFEs achieve good performance, there are many approaches to 
develop the control law for this system. Early solutions were linear regulators like 
the one proposed in [25]. It changes the modulation index slowly to regulate the 
DC-link capacitor voltage. Therefore, its main disadvantage is that it presents a slow 
dynamical response. A faster response can be obtained using a deadbeat current 
control [38]. However, it is well known that it is highly sensitive to the parameter 
uncertainties. In general, these approaches define an operating point and then work 
with a small-signal linearized model around it. This is a drawback because they 
can not guarantee stability against large signal disturbances for the large range of 
operating conditions of the three-phase PWM converter due to the fact that AFE are 
nonlinear systems. Furthermore, the controller implementation requires the system 
parameters which depend on the operating points, otherwise it may result in steady 
state errors in the state variables. On the other hand, in light of the strong nonlinearity 
of AFE, nonlinear control algorithms may be suitable for controlling the power 
converters, which are able to accommodate a wide range of operating conditions. 
The main reason is that there’s no need to have a linear model of the power converter 
for nonlinear controller design. 

Several nonlinear control approaches have been proposed for grid-tied power 
converters, such as input-output linearization [16], nonlinear adaptive control [20], 
passivity based control [5S, 10], feedback linearization [15], differential flatness based 
control [11, 24, 31], model based predictive control [33, 35] and SMC [30]. Among 
these techniques, SMC is suitable for dealing with the nonlinear behavior of the 
considered system due to its characteristics of insensitivity to external disturbances, 
system reduction, high accuracy and finite time convergence [7, 19]. 

SMC has been developed as a new control design method for a wide spectrum of 
systems including nonlinear, time-varying and fault tolerant systems [12]. SMC can 
manage the nonlinear behavior of the three-phase grid-connected power converter. 
Besides, it is characterized to be a robust and effective control strategy. However, 
up to authors knowledge, proposed SMC strategies have only considered the control 
of input current in sliding mode [21, 27]. In general, these works use a proportional 
integral (PI) controller for the DC-link capacitor voltage regulation. This approach 
achieves robustness of input current control, but can’t guarantee robustness of output 
voltage control since they are derived using approximations and linearizations. To 
solve this issue, extended state observer (ESO) based control strategies are proposed 
for voltage control design to reject disturbances and uncertainties [22, 34]. ESO is an 
efficient technique for disturbance estimation, which regards the lumped disturbances 
(such as parametric uncertainty, unmodeled dynamics and load variation) as a new 
state. However, only simulation results are provided in both works, and the stability 
analysis of the closed-loop system with modeling uncertainties is not performed. 

In this chapter, a model-based second-order SMC for three-phase grid-connected 
power converters is proposed. The controller design is based on the system model 
and has a cascaded structure which consists of two control loops. The outer loop 
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regulates the DC-link capacitor voltage and the power factor providing the current 
references for the inner control loop. The current control loop tracks the actual 
currents to their desired values. To design the proposed controller, the load connected 
to the DC-link capacitor is considered as a disturbance, which directly affects the 
performance of the whole system. Thus, a composite control law consisting of SOSM 
based on STA and disturbance compensation via ESO is developed for the voltage 
regulation loop. The current control loop is also designed using the STA. The STA 
is one of the most popular SOSM algorithms and a unique absolutely continuous 
sliding mode algorithm, ensuring all the main properties of first order SMC for 
systems with Lipschitz continuous matched uncertainties/disturbances with bounded 
gradients [8]. The modeling uncertainties are incorporated into the design of ESO and 
SOSM algorithms. The sliding surface is designed ensuring the finite time asymptotic 
convergence of the sliding variables to its desired values in the presence of parametric 
uncertainties. In addition, theoretical study of the stability issue is provided using 
Lyapunov method. 

To show the benefits of the proposed controller it can be compared with some 
conventional control strategies like PI synchronous reference frame (PI-SRF) or 
proportional plus resonant (PR) controllers. The PI-SRF is considered as the stan- 
dard control method in industry for this system and has been finally selected as the 
baseline controller. It should be noticed that the PR controller is also a linear control. 
Besides, when PR is adopted, a simple PI regulator is usually employed for the outer 
control loop. Therefore, it is expected that proposed ESO-SOSM provides the same 
improvements as for the PI-SRF. That is, better transient response and improved 
system performance under a load step. The performance of the proposed control 
is compared with a well-tuned conventional PI-SRF controller. The results show 
that the performance of the designed controller presents a faster dynamic behavior 
while maintaining a lower total harmonic distortion (THD) value in steady state and 
improves system performance under a load step. The validity of the proposed control 
algorithm has been verified by experimental results for a 3.0 kVA insulated gate bipo- 
lar transistor (GBT) PWM power converter using a TMS320F28335 digital signal 
processor (DSP). 


9.2 Problem Formulation 


9.2.1 Dynamic System Model 


The electrical circuit of the three-phase two-level grid-connected power converter 
under consideration is shown in Fig. 9.1. The system is connected to the grid through 
a smoothing inductor L with a parasitic resistance r. It is assumed that an equivalent 
resistive load R; is connected at the DC-link capacitor C. The load is considered as 
an unknown external disturbance. 
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Fig. 9.1 Circuit of the three-phase two-level grid-connected power converter 


The grid current and DC-link capacitor voltage dynamics can be written in a dg 
SRF rotating at the grid frequency [2]. Assuming that grid voltages are balanced then 
the system model is 


dig 
ery =—-rigt+ wig + Ug — bu Vacs (9.1) 
di, . . 
LT STi = tig ty 0, Vie (9.2) 
dVuc : . : 
C a = (Saia + dig) — lLoad» (9.3) 


where dz, 6g are the switching functions, vg, vg are grid voltages, ig, i, are input 
currents, Vz, is DC-link capacitor voltage and w is the angular frequency of the grid. 
From the control point of view, working on the dg SRF has the advantage of reducing 
the current control task into a set-point tracking problem. 

It should be pointed out that the amplitude of the control vector is constrained due 
to the fact that only implementable control vectors are contained in an area limited by 
the well-known hexagon. Therefore, it is necessary to consider the constraint || dag || < 
/2, which implies that the switching function stays unsaturated if its magnitude is 
not larger than ./2. This is a conservative approach but ensures the system operation. 


9.2.2 Modeling Uncertainties 


In practical applications, the system modeling is usually obtained under several 
assumptions, e.g. ignore the switching losses, in order to simplify the control design. 
However, such models are not exact, and may vary from real system behavior under 
some operation conditions. Therefore, a controller should be designed to be robust 
against any modeling parametric uncertainties. For this reason, parametric uncertain- 
ties, i.e. smoothing inductor L, parasitic phase resistance r and the angular frequency 
of the grid w, are fully considered and included into the model such that the control 
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robustness is guaranteed. These parameters are formalized as follows: 
L=Lo+tAL, r=rn+Ar, w= wo + Aw, (9.4) 


where Lo, 7, Wo are the nominal values and AL, Ar, Aw are the parametric uncer- 
tainties which can be considered as unknown slow variant signals. 


9.2.3 Control Objectives 


The control objectives for three-phase two-level grid-connected power converters 
are as follows: 


e The currents iz, i, should track their references i7 and i7, respectively. The refer- 
ence value i7 is calculated in such a way that DC-link capacitor voltage is regulated 
to a certain value. The reference value 7 is set to provide a desired instantaneous 
reactive power. 

ig 17, ig > ee (9.5) 


e The DC component of the DC-link capacitor voltage should be driven to some 
reference value V7... 
Vac > Vj. (9.6) 


9.3 Disturbance Observer Based SMC Design 


In three-phase two-level grid-connected power converters there exist different kinds 
of disturbances, like parameter uncertainties and load variations. If the controller 
does not have enough ability to reject these disturbances then they will degrade the 
performance of closed loop system. A cascade control structure is used to govern the 
system (9.1)—(9.3). The controller consists of a current tracking loop, inner loop, and 
an ESO-based voltage regulation loop, outer loop. For the current tracking loop, an 
STA controller is designed which ensures fast convergence of the currents ig andi, to 
their references i7 and i; , respectively. For the voltage regulation loop, an ESO is used 
to estimate the load power which is considered as an external disturbance. Besides, an 
STA controller is implemented in parallel to regulate the DC-link capacitor voltage 
to its desired value using the estimate of the disturbance. In this section, firstly, the 
basics of STA will be briefly shown. Secondly, the design of both control loops will 
be presented. 
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9.3.1 Recalling of Super-Twisting Algorithm 


SMC has been developed as a new control design method for a wide spectrum 
of systems including nonlinear, time-varying and fault tolerant systems [12, 14]. 
The STA is one of the most popular SOSM algorithms and a unique absolutely 
continuous sliding mode algorithm, ensuring all the main properties of first order 
SMC for systems with Lipschitz continuous matched uncertainties/disturbances with 
bounded gradients [8]. This part discusses the STA in a general case for a single input 
nonlinear system. 

The sliding mode design approach consists of two steps. The first step considers the 
choice of sliding manifold which provides desired performance in the sliding mode. 
The second step concerns the design of a control law which will force the system 
states to reach the sliding manifold in finite time, thus the desired performance is 
attained and maintained. This part discusses the STA in a general case for a single 
input nonlinear system. 

Consider a nonlinear system 


x = a(x)+bd(x,u), 


; oes (9.7) 


where x € X C R" is the state vector, u € U C Ris the input, s(t,x): Rt! > R 
is the sliding variable and a(x) and b(x, u) are smooth uncertain functions. 

The control objective is to force s and its time derivative s to zero. By differenti- 
ating the sliding variable s(t, x) twice, the following relations are derived: 


: O 0 
s= a (t,x) + aa (t,x) [a (x) + b(x,u)], 


i= o (xt os (t,x, u) [a (x) +b (x, w)] + a (t,x,u)a 08) 


= g(t, x,u)+ (t,x, u)u. 
Assuming that the sliding variable s has relative degree one with respect to the 


control input u, i.e. 23 (t,x, u) 4 0, there exist positive constant values ®, Ij, and 
I°y such that the following conditions are satisfied, 


0<In <y@t,x,u) < Ivy, (9.9) 
—-<y(t,x,u)<@. (9.10) 


Under the conditions (9.9) and (9.10), the following differential inclusion can be 
obtained: 
s €[-@,4+46]4+ Un, Mula. (9.11) 


In the sequel, a control law based on STA is designed. It consists of two terms, one 
is the integral of its discontinuous time derivative while the other is a continuous 
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function of the available sliding variable s. 

u=uU; +2, & = —asign(s), wu. = —A|s|? sign(s), (9.12) 
where a and A are design parameters that can be determined from the boundary 


conditions (9.9) and (9.10). The sufficient conditions for finite time convergence to 
the sliding manifold s = s = 0 are [18]: 


ce) , 4PIuat+? 
- ~ 12 T,a-@ 


m 


: (9.13) 


Remark 9.1 In the case of relative degree one systems, traditional first-order sliding 
mode could also be applied. However, motivated by the chattering elimination aim, 
super-twisting based control is employed which means that the control signal u is 
continuous and chattering is avoided. 


9.3.2 Extended State Observer Design 


According to the control objective (9.6), the outer control loop is designed to regulate 
the output capacitor voltage to its reference value V7... For the system (9.1)-(9.3), 
the current dynamics are much faster than the DC-link capacitor voltage dynamics 
[32]. Under this condition, it can be considered that ig ~ i7 and i, ~ i7. Based on 
the singular perturbation theory [13], if the fast dynamics are stable then (9.3) will 


be reduced to Pr i 
dc 
C = 
dt Vue ( 


Dp Pload) , (9.14) 


where p* = vgi7 + Ugig and Pjoad = Vaclload- 

It should be noted from (9.14) that pjoaq can be considered as an external distur- 
bance. This paper proposes to design an ESO to estimate the disturbance asymptot- 
ically. Then, it will be injected into the control design. Unlike traditional observers, 
such as Luenberger observer [23], high-gain observer [1] and UIO [4], ESO regards 
the disturbances of the system as new system states which are conceived to estimate 
not only the external disturbances but also plant dynamics [9, 28, 29]. 

To design the ESO, the new variable z = V2. /2 is introduced in (9.14), yielding 


Ci = p*—d(t), (9.15) 


with d(t) = Pioaa. A linear ESO is given by 
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C2 = p*-d(t)+ Ac), (9.16) 
d(t) = -h(-2, (9.17) 

where the positive gains 3, and (3 are chosen such that the polynomial 


fen Bo 
b aa aay er oe 9.18 
+ c + Fe (9.18) 


is Hurwitz stable. Therefore, its natural frequency w, and damping ratio € are: 


_ JB -«. Pr jt 
rene c= 2V BC (9.19) 


Denote the observation errors €, = z — Z, €g = d(t) — d (t), the error dynamics are 
given by, 


Ce S66 =e; (9.20) 
€q = foez +h(t), (9.21) 


where h(t) = d (t) is the variation rate of load power. The system (9.20) and (9.21) 
can be written as follows: 
€= Ae+y, (9.22) 


fa 
where ¢ =[e,, eg]", A= BC f| and a = [0 A(t)]'. 


Lemma 9.1 Suppose that h (t) is bounded, there exist a constant 6 > 0 and a finite 
time T, > O such that the trajectories of the system (9.22) are bounded, |\e|| < 6, 
Vt = T, > 0. 


Proof Solving (5.15), one has 
t 
E(t) = ee (to) + i ef a) de, (9.23) 
i) 


where fy is the initial time. Using the bound jet | < ke“), with k > 0, 
(6 > 0, then (9.23) can be estimated as follows: 


t 
lle (t)I| < ke Je (to) I + . ke PO IIb (1) || dr 


to 


Sone k 
< ke“? |e (to) || + ees IIb (71. (9.24) 
to<T<t 
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It follows from (9.24) that |le|| < 6, Vt => 7, > 0, where 6 is a positive constant that 
depends on k, ( and the upper bound of ||~ (7) ||. | 


Remark 9.1 In view of (9.20) and (9.21), the parameters 3, and (3, determine the 


bandwidth of the ESO, i.e. 
Wh = Wry (4é4 — 4€? + 2). (9.25) 


Generally speaking, the larger the ESO bandwidth is, the more accurate estimation 
will be achieved. However, this increases the noise sensitivity due to the augment of 
the bandwidth. The function h(t) represents the rate of change of load power. If this 
value is quite large then it means that the load power changes very rapidly. In this 
case, the observer bandwidth needs to be sufficiently large for an accurate estimate 
of d(t). Therefore, the selection of (3, and G2 should balance between the estimation 
performance and the noise tolerance. 


9.3.3 Capacitor Voltage Regulation Loop 


Define the regulation error zZ = z* — z with z* = (vi) /2, it follows that, 
Cz = —p* +d(t). (9.26) 


Notice that the perturbation d(t) is an unknown time varying variable. The proposed 
ESO-based STA controller for the voltage control loop is given by 


p= 1,2) +d(0), (9.27) 


in which jige(Z) is the STA which takes the following form, 
1 t 
Mac(Z) = Ade|Z|? sign(Z) + au f sign(z)dT, (9.28) 
0 
with some positive constants Ay, and ag,. Substituting (9.27) into (9.26), yields 
CZ = —HMac(Z) + €a. (9.29) 


It can be easily obtained from Lemma 9.1 that 


lleall < |All 6+ sup |lb(7)Il = Fa, (9.30) 
to<T<t 


with ||A|| = Amax(A™A) and Fy is a positive value. The sufficient conditions for 
the finite time convergence to the sliding manifold z = z = 0 are [18]: 
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Fig. 9.2 Block diagram of the STA + ESO controller for the voltage loop 


Qdc + Fg 


Qc > CFy, . > 4C7F, . 
dc d ee ges Ey 


(9.31) 


The structure of the ESO-based voltage regulation loop is shown in Fig. 9.2. 


9.3.4. Current Tracking Loop 


The objective of the inner control loop is to force iz and i, to their references i7 and 
i; respectively. The reference current i 7 is computed from the output of the outer 
control loop and is calculated to achieve dynamic voltage regulation. The reference 
current i7 is set to provide a desired instantaneous reactive power q*. Therefore, the 
current references are calculated as 


f@a tt, ita, (9.32) 


Sq =i5 — ia, Sq = 17 — ig. (9.33) 


Taking the first time derivative of su, = [sa 5q]" yields, 


Ae | cs) er cep te a Nae i (9.34) 
Sq i, lg = pele L |g 


In order to satisfy the saturation constraint, the controllers dg and 6, are designed as 
follows: 


6g =a (mg), (9.35) 
0g ="a° (mq) ‘ (9.36) 
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where 
Lo Ve To. : . 
ee val — pala) + = Pia = i + wig |, (9.37) 
Lo u’, YO. : ; 
ig = al = Ug (Sq) + 7m Lyi! be woiu], (9.38) 
[a(Sa) and jg (Sq) are STAs which are in the form 
; t 
yta(sa) = ralsal !sign(sa) + og / sign(sg)dr, (9.39) 
0 
t 
i - 
fg (Sq) = Ag|Sq|? Sign(sg) + ay | sign(s,)dT, (9.40) 
0 


with some positive constants 4;, a;, i € {d,q}, and o(x) is a saturation function. 


Theorem 9.1 Consider the system (9.1)-(9.3) in closed loop with the saturated 
controller (9.35) and (9.36). This yields 


: Lo 
Sa = 5 Ma) + galt), (9.41) 
: Lo 
Sq = = HG) + Gq (t), (9.42) 


where 


Ar AL. woly -— wl. 


galt) = Z ig + i ig + : igs (9.43) 
Ar. AL. woly -— wl. 
q(t) = Tia + 7 is : ig. (9.44) 


The state trajectories of the system (9.41) and (9.42) converge to the origin sq = 
0, Sg = O in finite time if the gains of (sq), [4(Sq) and Vj. are chosen such that the 
following conditions are satisfied, 


Qa > ; a2 > aa, 
: ; (9.45) 
Qg >5 : ra Qg; 
and 
2 
Ve 2 (Lolli, Il) 4 3E?, (9.46) 


where Yo, Ya and Yq are positive constants and E is the amplitude of the grid source. 
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Fig. 9.3. Block diagram of the current tracking regulation loop 


Proof The proof is divided into two parts. Firstly, if the control vector ||duq || < wf, 
then Jaq = [ma, mq]'. Secondly, in the case when the control vector ||uaq || > af, 
thus the control vector is saturated to that value. 

Case 1: The controller is given by 6g = mg and 6, = mg in (9.35) and (9.36). 
Considering closed loop behavior (9.41) and (9.42) and given that (9.1)-(9.3) is a 
physical system, thus it is reasonable to assume that the variables yy(t) and ¢, (t) 
and its time derivatives are bounded functions: 


Iga < ya. Ilee@Il < Ya, (9.47) 


with some positive constants yg and y,. Taking into account that 
A-wL<los d+ wk, (9.48) 


for some scalar 0 < yo < 1. Then, according to [18], the trajectories of the system 
(9.41)-(9.42) will converge to sy = Sg = 0 and s, = 5, = O in finite time when the 
controller gains satisfy conditions (9.45). 

Case 2: According to [6], the sufficient condition for the control vector da, to enter 
into the circle of radium /2, (i.e., []dagl| < /2) is that V7. satisfies the condition 
(9.46). Thus, Theorem 9.1 is proven. a 


Remark 9.2 It should be noted that the cross-coupling terms and the source volt- 
age are compensated by (9.35) and (9.38). However, (9.37) and (9.38) require the 
information of the derivative of i; and 7. Usually, these values are zero during the 
steady state and only affect slightly the system performance during the transient state. 
Therefore, these terms are neglected in the final control law. Similarly, parameter ro 
is very small in order to reduce system losses. Moreover, this parameter is usually 
unknown and its use it is also avoided. Taking into account these considerations, 
the block diagram for the current tracking control loop is shown in Fig.9.3, where 
Udg = Vacddq: 
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9.4 Experimental Results 


In order to demonstrate the feasibility of the proposed control algorithm, practical 
results have been performed, comparing the proposed ESO-based SOSM control to 
a well tuned linear conventional PI-SRF regulator. The PI-SRF is considered as the 
standard control method for this system and has been selected as a baseline controller. 
The electrical parameters of the power converter, DC-link voltage reference, and 
switching and sampling frequencies for the experimental setup are summarized in 
Table 9.1. 

The power converter prototype used for the experiment is shown in Fig.9.4. A 
digital implementation of both current and DC-link voltage control algorithms is 
executed in a TMS320F28335 floating point digital signal processor board. Two sets 
of experiments are done. The first one consists of a load step at DC-link from no load 
to full load (3.125 kW). To perform the experiment, the capacitor voltage reference 
is set to 750 V and a 180 £2 resistive load is suddenly connected to the DC-link. 
The second test focuses on the reactive power tracking ability. For this purpose, 
an instantaneous reactive power command step is done. Measurements of DC-link 
voltage, phase voltages and currents, harmonic contents of currents, active power, 
reactive power, and power factor have been taken. 

The parameters of the PI-SRF and ESO-SOSM controllers are given in Table 9.2. 
They are chosen so that the current dynamics are much faster than the capacitor 
output voltage dynamics. It should be noticed that a Phase Locked Loop (PLL) is 
also implemented in the digital platform in order to work in the SRF [36]. 


Table 9.1 Electrical system 


Parameter Description 
parameters 
Phase-to-neutral voltage (RMS) 230 V 
Grid frequency wo 50Hz 
Filter inductor Lo 15 mH 
DC-link capacitor C 2800 uF 
DC-link voltage reference V7, 750 V 
Sampling frequency f; 10 kHz 
Switching frequency fs, 10 kHz 
Table 9.2 Controller design PI-SRF Values 
parameters ; 
Voltage regulation loop (Kpdc, kidc) 0.04, 0.5 
Current tracking loop (Kpa, kKpg, kia, Kig) | 75, 400, 75, 400 
ESO-SOSM Values 
Voltage regulation loop (Age, Ade, 1, 2)| 3, 750, 3, 300 
Current tracking loop (Ag, Qa, Ag, q) 85, 20000, 85, 20000 
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Fig. 9.4. Power converter prototype: a Power converter, b DC load 
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9.4.1 Load Step at DC-Link 


The first test consists of evaluating the proposed controller performance under a 
load step at the DC-link capacitor. Figure 9.5 shows the transient response of the 
DC-link capacitor voltage for a load step from no load to a load composed by a 
resistor of 180 92. Three waveforms are presented. Figure 9.5a corresponds with the 
conventional approach with a PI-SRF controller. Figure 9.5b is the result achieved 
with the proposed ESO-SOSM algorithm and Fig. 9.5c is associated to the PI-SRF 
but increasing the bandwidth of the controller. 

The results show that, both control laws can achieve the DC-link capacitor voltage 
regulation. Comparing Fig. 9.5a with Fig. 9.5b, the settling time is roughly the same 
when PI-SRF and ESO-SOSM are used but the last one reduces the DC-link capacitor 
voltage drop. The PI-SRF results in a voltage drop of 45 V while in the case of the 
ESO-SOSM approach this is only 22 V. Therefore, the proposed controller reduces 
the DC-link capacitor voltage drop under a load step in 48.9%. Clearly the ESO- 
SOSM requires less energy from the DC-link capacitors compared with the PI-SRF. 
Same DC-link voltage drop can be achieved with the PI-SRF but it is necessary 
to increase the outer controller bandwidth. It is possible to do this, but it is well 
known that increasing the bandwidth of the outer control loop in the conventional 
PI-SREF affects the grid current harmonic content. Table 9.3 collects the current THD 
computed up to 50th harmonics for the three situations. Clearly the ESO-SOSM 
allows to achieve a better transient response while maintaining a low THD in steady 
state. 


(a) (b) (c) 
Tek JU @ Acq Complete M Pos: 300.0ms Tek JU @ Acq Complete M Pos: 3000s Tek JUL @ Acq Complete _M Pos: 300.0ms 
* + + 
/ — ’ fcnastiets 
| gees re 3 
« 
/ | rT ] « 
\e ft 
Chl 100¥ M 1000s On J Chl 10.0¥ M 100ms CH J CHT Wa¥ M 100ens ral 


Fig. 9.5 Transient response of DC-link capacitor voltage under a load step: a PI-SRF algorithm, 
b ESO-SOSM strategy, ¢ PI-SRF algorithm with increased bandwidth 


Table 9.3. Grid current THD for the load step experiment 
Controller 

PI-SRF 

ESO-SOSM 

PI-SRF ((Kpde, Kide) * 3) 
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Fig. 9.6 Transient response of grid currents in the dq frame under a load step. Top: PI-SRF algo- 
rithm. Bottom: ESO-SOSM strategy 


Figure 9.6 plots the grid currents for the load step at the DC-link capacitor tests. 
The waveforms are represented in the SRF. Only the currents for Fig. 9.5a and b are 
plotted. Clearly the PI-SRF presents a slow transient response and needs some time 
to reach the steady state value. On the other hand, the ESO-SOSM approach quickly 
changes the iz current component value. It should be noticed that 77, comes from the 
outer control loop. Therefore, the ESO introduced in the control law provides a clear 
improvement compared with the conventional approach that only considers a PI for 
the DC-link voltage regulation. 


9.4.2 Reactive Power Command Step 


The second experiment assesses the control law capability to provide a desired 
amount of reactive power. To do so, it is introduced a instantaneous reactive power 
command step from 0 to 3 kVAr. The system response is evaluated for both the 
PI-SRF and ESO-SOSM controllers. 

Figure 9.7 presents the transient state for the experimental results achieved by the 
PI-SRF and Fig. 9.8 shows the corresponding one to the ESO-SOSM. As expected, 
both controllers provide the desired reactive power. However, the dynamic response 
for the ESO-SOSM is faster than the PI-SRF. The performance in most interest is the 
steady-state behavior. Figure 9.9a—c show the steady state features for the PI-SRF 
controller. Figure 9.9a—c are the grid currents, the current spectrum and information 
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Fig. 9.7 Transient response of grid currents in the dq frame under a load step for the PI-SRF 
algorithm. Top: Grid currents in the SRF. Bottom: Instantaneous active and reactive power 
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Fig. 9.8 Transient response of grid currents in the dq frame under a load step for the ESO-SOSM 
strategy. Top: Grid currents in the SRF. Bottom: Instantaneous active and reactive power 
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Fig. 9.9 Steady state currents for a instantaneous reactive power command of 3 kVAr: a Current 
waveforms for PI-SRF algorithm, b Current spectrum for PI-SRF algorithm, ¢ Power factor for PI- 
SRF algorithm, d Current waveforms for ESO-SOSM strategy, e Current spectrum for ESO-SOSM 
algorithm, f Power factor for ESO-SOSM strategy 


about the active power, reactive power values and the power factor respectively. On 
the other hand Fig. 9.9e-f present the same information in the case of the ESO-SOSM 
approach. Comparing the grid spectrum in Fig. 9.9b and e the ESO-SOSM controller 
has better performance. Particularly, the Sth and 7th harmonics are reduced compared 


9.4 Experimental Results 203 


with the PI-SRF results. Furthermore, the current THD is 1.2% for the ESO-SOSM 
and 1.7% for the PI-SRF. This supposes a reduction of 29.4% when the ESO-SOSM 
approach is used in the same power converter prototype. 


9.5 Conclusions 


In this chapter, the problem of disturbance observer based SMC of three-phase two- 
level grid-connected power converters has been investigated. SOSM technique has 
been adopted to solve its control problem, which is considered as a promising alter- 
native due to its features of robustness and effectiveness for nonlinear systems. 

Considering the dynamics of the output voltage, the load connected to the DC-link 
capacitor is regarded as a disturbance for the voltage control loop, which directly 
affects the performance of the whole system. To improve the disturbance rejection 
ability, a composite control law consisting of SOSM based on STA and disturbance 
compensation via ESO is proposed for the voltage regulation loop. With the dis- 
turbance compensation using ESO, the gains of STA for the outer control loop can 
be reduced without decreasing the settling time and reducing the DC-link capacitor 
voltage drop. The current loop has been designed to track the current references in 
the presence of system parameter uncertainties using an STA. Experimental results 
have demonstrated that the proposed ESO-SOSM controller performs better than that 
of the conventional PI SRF control. The power converter operated with the proposed 
algorithm has achieved less DC-link voltage drop under a sudden load step, shorter 
settling time and better grid current quality in terms of lower THD and reduced values 
of low-order harmonics content. 
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